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1. Einleitung 


Die Regulierung ist ein Produkt des technischen Zeitalters. Sie wurde aber 
erst entwickelt, als es notwendig wurde, den Menschen fiir die Uberwachung 
eines maschinellen Vorganges zu ersetzen. So entstand zuerst die Tourenzahl- 
regulierung der Kolbendampfmaschine mittels des Fliehkraftpendels in der 
ersten Halfte des 18. Jahrhunderts. Die Starkstromtechnik bedurfte weiterer 
selbsttatiger Uberwachungseinrichtungen, wie Spannungsregler, Leistungsreg- 
- ler und andere mehr. In der Folge wurden immer mehr Vorgange vollautoma- 
tisiert, zum Beispiel Gliihbehandlung nach einem vorgegebenen Programm im 
Fabrikationsprozess, selbsttatige Kurshaltung eines Schiffes oder Flugzeuges 
(Autopilot), automatische Dampfungsregulierung eines Ubertragungssystems. 
Das sind nur drei von vielen Beispielen. In jedem Gebiet waren Spezialisten mit 
der Berechnung und Entwicklung der besondern Regulierungsprobleme be- 
schaftigt. Stabilitatsfragen mussten im Zusammenhang mit der Differential- 
gleichung des vollstandigen Systems diskutiert werden. Dies erforderte die 
Kenntnis der Lésungen der charakteristischen Gleichung. Von STopora, der 
das Regulierungsproblem bei der in Entwicklung begriffenen Dampfturbine 
bearbeitete, angeregt, gelang es Hurwitz Kriterien aufzustellen, mit denen 
entschieden werden konnte, ob eine charakteristische Gleichung nur Lésungen 
mit negativen Realteilen enthalte oder nicht. Im ersten Fall ist das System 
stabil. Die Berechnung der Lésungen ist aber schon bei einfachen Systemen 


: 
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langwierig und undurchsichtig. Erst mit der Abstraktion, dass jeder zeitliche 
Vorgang durch eine Superposition unendlich vieler und unendlich kleiner expo- 
nentiell zu- oder abnehmender Sinusschwingungen dargestellt werden kann, 
gelang eine Algebraisierung der Berechnung von solchen Problemen. Die Fern- 
meldetechnik war schon von Anfang an genétigt, mit Frequenzspektren zu 
rechnen. Es ist deshalb gut zu verstehen, warum die bedeutendsten Beitrage 
gerade von dieser Seite kamen (HEAVISIDE, CARSON, K. W. WAGNER, CAMP- 
BELL, Nyquist, BopE, KUPFMULLER). Dank der heute vorhandenen allgemei- 
nen Methodik sollte es jedem Ingenieur méglich sein, seine speziellen Regulier- 
probleme lésen zu k6nnen. 


2. Begriffe 


Unter einer Regulierung (in Deutschland Regelung) versteht man eine Ein- 
richtung, welche selbsttatig dafiir sorgt, dass eine gewiinschte Grésse (Span- 
nung, Temperatur, Druck, Drehzahl, Winkel, Drehmoment, Menge usf.) fort- 
laufend mit einem Sollwert verglichen wird und bei einer Stérung einen solchen 
Vorgang auslést, dass die Abweichung in vorgeschriebene Grenzen zuriickge- 
fiihrt wird. Aus dieser Definition ergibt sich das allgemeine Prinzipschema einer 
Regulerung (Figur 1). 

Meist handelt es sich darum, dass die Regulierung in einem bestimmten 
Bereich der regulierten Grésse wirksam sein soll. Der Zusammenhang zwischen 
dem Istwert und dem Fehler ist in sehr vielen Fallen nichtlinear. Um trotzdem 


Storung 
| 
1 
: 
i oo i: Stérungsumformer 
| | 
Vergleichs - tL. 4-4 
mittel | 


Sollwert Istwert 


Sollwert - 
Umformer 


Regler - Reguliertes 
element System 


Reguiierte 
Grosse 


Messwertumformer 


Fig. 1 
Allgemeines Schema einer Regulierung. 
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die Regulierung mathematisch behandeln zu kénnen, werden nur kleine Sté- 
rungen betrachtet. Selbstverstandlich ist die Betrachtung iiber den ganzen 
Regulierbereich auszudehnen. Der Reguliervorgang selbst ist ein zeitlicher Ab- 
lauf, der bei linearen Systemen aus der Frequenzcharakteristik und der gege- 
benen Stérfunktion berechnet werden kann. Sind verschiedene Stérfunktionen 
(Eingange) vorhanden, so addieren sich die Wirkungen auf den Istwert. Fiir 
jeden Eingang (Sollwert, Stérung) lasst sich das Schema der Figur 1 auf das- 
jenige der Figur 2 reduzieren, wobei die komplexen Ubertragungsfunktionen 


Vergleichs - Istwert 
Wert 
Messwert 
Fig. 2 Fig. 3 
Vereinfachtes Schema einer Regulierung. Einfachste Form einer Regulierung. 


R Vergleichswert (Reference Input); E Fehler (Error); C Messwert (Controlled Variable). 


H, und G, verschieden ausfallen, je nachdem der eine oder andere Eingang 
betrachtet wird. Das Produkt H, G, wird aber in allen Fallen die gleiche Funk- 
tion ergeben. Betrachtet man den Messwert C als Ausgangsgrosse, so entsteht 
schliesslich das Schema der Figur 3, wobei G = H, G,. Damit aber die regulierte 
Grésse den gewiinschten Verlauf nimmt, muss der Messwertumformer genau 
proportional sein, wobei das Auftreten bestimmter Zeitkonstanten im allgemei- 
nen nicht vermieden werden kann. 

Die drei Gréssen R, E und C sind die Zeiger von sinusoidalen Zeitfunktionen 
mit der Frequenz w. G ist die komplexe Ubertragungscharakteristik des auf- 
geschnittenen Kreises: 


C(b) = E@) Gib) (fir p=jo) 
E() = Rb) — Cl), (2) 
Cb) = ey RP) = ME) RE), 3) 
MO) = 7 e5 () 


M ist die komplexe Ubertragungsfunktion des geschlossenen Kreises. Sie kann 
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aber auch als Frequenzspektrum einer zeitlichen Funktion m/(t) aufgefasst wer- 
den, die sich dann ergabe, wenn als Zeitfunktion des Sollwertes die Impulsfunk- 
tion 6(t) gewahlt wiirde. Bei einer beliebigen, vorgegebenen Zeitfunktion des 
Sollwertes 7(¢) kann der zeitliche Ablauf des Messwertes c(¢) und damit verkniipft 
der Istwert der regulierten Grésse berechnet werden aus dem Faltungsintegral : 


t 


c(t) = m(t) * r(t) = [me r(t— 1) dr, (5) 
m(t) = af | Me) ertap. (6) 


m/(t) darf nur abklingende exponentielle Glieder aufweisen (Stabilitatsbedin- 

gung). 

Beim Entwurf einer Regulierung sind folgende Punkte zu beachten [1]?): 

a) Vollstandige Darstellung der gewiinschten Bedingungen, wie Regulierungs- 
bereich, maximal zulassige Abweichung des Istwertes vom Sollwert bei allen 
méglichen Stérwerten, zeitlicher Ablauf des Regulierungsvorganges bei 
plétzlicher Anderung des Sollwertes und von Stérwerten. 

b) Die Bedingungen unter a) miissen in einer Fréquenzcharakteristik M des 
geschlossenen Kreises ausgedriickt werden. 

c) Aus b) soll eine approximative Frequenzcharakteristik G des aufgeschnitte- 
nen Kreises abgeleitet werden kénnen. 

d) Bestimmung der statischen und komplexen Ubertragungscharakteristiken 
der einzelnen im Regulierkreis vorkommenden Organe, wie Motoren, Ver- 
starker, Ventile, Messfiihler usf., die bereits vorhanden sind. 

e) Auswahl der geeigneten, zusatzlich notwendigen Elemente, wie Verstarker, 
stabilisierende Netzwerke, Einfiihrung weiterer Gegenkopplungen usw., 
damit die unter c) berechnete Frequenzcharakteristik G zustande kommt. 

f) Vereinfachen mit dem Ziel einer 6konomischen Lésung des Problems. 

g) Experimentelle Priifung, insbesondere des zeitlichen Ablaufs der Regulie- 
rung bei plétzlichen Anderungen der verschiedenen Einflussgréssen. 

Aus der modernen Anwendung der Regulierung hat sich eine gewisse Typi- 
sierung herausgeschalt. 

Typ 0: Eine konstante Abweichung E hat einen konstanten Istwert der regu- 

lierten Grdésse zur Folge. 

Typ 1: Eine konstante Abweichung EF hat einen konstanten Wert der ersten 

zeitlichen Ableitung des Istwertes der regulierten Grésse zur Folge. 

Typ 2: Eine konstante Abweichung E hat einen konstanten Wert der zweiten 

zeitlichen Ableitung des Istwertes der regulierten Grésse zur Folge. 


1) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis auf Seite 259. 


a 
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Da der Messwert proportional zum Istwert sein muss, kann diese Typisie- 
rung auch formelmassig zum Ausdruck gebracht werden: 


: SOP) 
Le a 


Typ 02 G(0) = (K = Konstante) , 


Typ 1: G(p) = 


, 


K 
a 
p 
rs 
UL SOS oe 


Das statische Verhalten (p> 0) der Regulierung 


G 1 


Pomel Go eae 


ist infolgedessen bei Typ 0 charakterisiert durch eine dauernde Abweichung des 
Istwertes vom Sollwert, die um so kleiner wird, je grésser die Verstarkung K 
ist. Dagegen wird bei den tibrigen Typen M = 1, da 1/G gegen Null strebt. Die 
Abweichung wird dann ebenfalls Null. Fiir Positionsservosysteme werden des- 
halb nur die Typen 1 und 2 verwendet. 

Ferner wird zwischen kontinuierlicher und diskontinuierlicher Regulierung 
unterschieden. Die letztere kann in vielen Fallen durch Mittelwertbildung auf 
die erste zuriickgefiihrt werden (zum Beispiel Tourenzahlregulierung an Klein- 
motoren durch Fliehkraftunterbrecher). Besondere Vorsicht ist bei der experi- 
mentellen Ausmessung geboten. Um mit der Theorie vergleichbare Resultate 
zu erhalten, diirfen im allgemeinen nur kleine plétzliche Anderungen des Soll- 
wertes vorgenommen werden. Ahnlich verhalt es sich mit plétzlichen Anderungen 
der Stérwerte, da im geschlossenen Kreis an einzelnen Organen Sattigungs- 
erscheinungen infolge von Ubersteuerungen eintreten kénnen. In der folgenden 
methodischen Betrachtung werden solche Sattigungserscheinungen nicht be- 
rticksichtigt. 


3. Stabilitatsbetrachtung 


Ausgehend von der Frequenzcharakteristik M(p) des geschlossenen Regu- 
lierkreises bei gentigend kleinen Anderungen der verschiedenen Gréssen ist das 
Kriterium fiir Stabilitat leicht anzugeben. Fiir jeden Pol von M muss Re (f)<0 
sein. Nyquist [2] hat daraus ein allgemeines Kriterium fiir die Frequenzcharak- 
teristik G(p) des aufgeschnittenen Kreises abgeleitet, damit der geschlossene 
Kreis stabil sei. Es lautet folgendermassen: 

Hat die Funktion G(p) Pole in der rechten Halbebene von p, so muss die 
Abbildung der (j w)-Achse in der G-Ebene den Punkt —1-+ 70 n-mal im 
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Gegenuhrzeigersinn umschlingen, wenn alle Frequenzen von —@ bis +@ durch- 
laufen werden, damit das System, charakterisiert durch M =G/(1+ G), stabil — 
ist. Da fiir alle physikalisch realisierbaren Ubertragungsfunktionen gilt: 


Re G(—j w) = ReG(+j7@) und ImG(—j)=ImG(+j@), 


so liegen die Polstellen symmetrisch zur reellen p-Achse. Es gentigt die Abbil- 
dung der positiven (j w)-Achse in der G-Ebene. Der Polstrahl vom Punkt 
—1-+ 70 zu dieser Ortskurve muss von w= 0 bis @ = +co den Winkel 
+ (n/2) 2m = +2 tiberstreichen, damit das geschlossene System M stabil ist. 


Fig. 4 
Beispiel einer stabilen Regulierung, deren offener Kreis unstabil ist. 
K pK G 
Ge) ==, M(p) 


p?—2Ept+ 1’ — p+ (K—28p+1 1+G 


fiir > 0 offener Kreis instabil, fiir 4 — 2§> 0 geschlossener Kreis stabil 


Es kann also vorkommen, dass der aufgeschnittene Kreis instabil, der ge- 
schlossene Kreis dagegen stabil ist. Ein Beispiel dafiir ist eine Anordnung nach 
Figur 4. Fiir Regulierungen kommen praktisch nur stabile offene Kreise in 
Frage, so dass das Stabilitatskriterium wesentlich einfacher gefasst werden 
kann. Weist G(p) keine Pole mit positivem Realteil von # auf, so geniigt die Be- 
dingung, dass die Abbildung der (7 w)-Achse auf die G-Ebene den Punkt 
—1-£70 nicht umschlingen darf, damit der geschlossene Kreis stabil sei, oder 
noch praziser, dass der Polstrahl vom Punkt —1-+70 nach der Ortskurve 
beim Durchlaufen von —co < @< +oo eine totale Drehung von 0° erfahrt 
(siehe Figur 5). 

Die Abbildung der (j w)-Achse in die G-Ebene wird Nyquist-Diagramm 
genannt. Figur 56 zeigt einen Fall, der bei Verringerung der Verstarkung un- 
stabil wird. Auch dies ist bei Regulierungen méglichst zu vermeiden. 

Das Stabilitatskriterium vereinfacht sich dadurch auf folgende Forderung: 

Der geschlossene Kreis ist stabil, wenn die Frequenzcharakteristik G(j w) 
des aufgeschnittenen Kreises bei allen Absolutwerten von G>1 eine kleinere 
Phasendrehung als +- 180° aufweist. 
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wy, 


b) 


a Gegengekoppelter Wechselstromverstarker, stabil; 
b stabile (bedingte) Regulierung, Typ 0; 
ce stabile Regulierung, Typ 1. 


Fig. 5 
Verschiedene Formen des Nyquist-Diagramms. 


Der Zusammenhang zwischen dem Verlauf des Absolutwertes und der 
Phase von G(j w) fiir realisierbare Frequenzgange lieferte neuartige Berech- 
nungsmethoden, welche das Stabilitatsproblem auf einfache Weise zu iiber- 
blicken gestattet. BopEs grundlegende Arbeiten [3], [4] zeigten dazu den Weg. 

Die komplexe Ubertragungsfunktion G kann durch ein komplexes Uber- 
tragungsmass dargestellt werden, das ebenfalls von der Frequenz abhangig ist. 
‘Das vereinfachte Stabilitatskriterium im komplexen Ubertragungsmass ausge- 
driickt, lautet dann 

A=In|G| 20, —180°<B<+180° 


(8) 
A<0O, B keine Beschrankung. | 


Es hat sich als zweckmassig erwiesen, die Winkeldifferenz von B gegen 
+ 180° als Phasenspielraum y einzufiihren. 
Stabilitat ist dann vorhanden, wenn gilt: 


Ave Oley 10. Ax<0 | oy -behiebigs (9) 


in Figur 5a und 5c erfiillt, in Figur 50 nicht erfiillt. 
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BopE untersuchte die Bedingungen, unter welchen sich eine eindeutige Zu- 

ordnung zwischen den reellen Funktionen A(@) und B(w) ergibt. A heisst die 

Verstairkungsfunktion und B die Phasenfunktion. Bei gegebener Verstaérkungs- 

funktion A kann die Phase B unter folgenden Voraussetzungen berechnet 

werden: 

a) F(j w) = A(w) +7 B(w). F(p) darf in der rechten p-Halbebene keine Pole 
enthalten. 

b) Auf der (j@)-Achse diirfen nur solche Singularitaten vorkommen, dass 
lim (p — po) F(p) zu Null wird. 


b—Po 
c) F(p) soll im Unendlichen analytisch sein, das heisst 
lim 7?) 9. 
poo P 


d) Die Funktion F(f) soll im Nullpunkt und im Unendlichen durch folgende 
asymptotische Reihen darstellbar sein: 


$>0 F=A,+joB,+4,09+ 7B, 0+ -- 


A‘ Bi 
ree 21 


proo F=A,+j 224 + 


7) w? ws 


Diesen Bedingungen entspricht beinahe jede Ubertragungsfunktion, insbe- 
sondere auch der Logarithmus davon, welche physikalisch realisierbaren Ele- 
menten zugeordnet sind. Gleichbedeutend mit den obigen Bedingungen ist die 
Forderung, dass sich das dynamische Verhalten des offenen Kreises in einer 
linearen Differentialgleichung endlichen Grades mit konstanten Koeffizienten 
ausdriicken lasst. Elemente, deren Verhalten mit partiellen Differentialglei- 
chungen beschrieben werden kénnen, entsprechen nicht den genannten Voraus- 
setzungen. Man kann ihr Verhalten naherungsweise beriicksichtigen, wenn es 
moglich ist, die partielle Differentialgleichung in eine gewohnliche zu verwan- 
deln. Physikalisch heisst das, Konzentrieren von homogen verteilten Eigen- 
schaften. 7 

Es gilt, wenn obige Voraussetzungen erfiillt sind, 


Blo,) = -2f. [A ( (w) — ae Wy dw 


— a 


(10) 


Bei Einfiithrung einer logarithmischen Frequenzskala u = In (w/w ) kann diese 
Beziehung umgeformt werden auf 


Bli\ ed Hi in (ctgh 41) du (11) 
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oder 
+00 
pio = (44), + 2 [LGA — (28), Jo (onl) ow. a 


Da die Funktion In (ctgh||/2) symmetrisch zu «= 0 mit wachsendem wu 
rasch klein wird, liefert das erste Glied den Hauptbeitrag zur Grésse der Phase B 
bei der Frequenz wy. Geniigt die Ubertragungsfunktion G(p) den vorstehenden 
Bedingungen, so ist sie eine rationale Funktion von der Form wie zum Beispiel 


7 K+ 7) (1+?) 
Spit r hia pimiaprs2ete SF SOS 


Die Verstarkungsfunktion 


Ajo) =InK+n|j1+joT,|+n|1+j7o0%|—-In|jo| vee 
—In|1+jT,|—Inj1+jo7,|—In|1—@?T?+2jT¢|. | 
Die Phase 
B(w) = arctgw T, + arctgw T, — s 
(15) 
T; t Tea tg ee 
—arctgw J; —arctgw I, Sr 


Wird der FrequenzmaBstab logarithmisch gewahlt und sowohl A(w) wie B(w) 
aufgezeichnet, so erkennt man sogleich die Bedeutung des ersten Gliedes von 
Formel (12). Die Glieder A, = In|1+ 7 T,| sehen alle gleich aus (Figur 6): 


10° ie: 


Fig. 6 
Verstarkungs- und Phasenmass der Ubertragungsfunktion 1 + 7@ T. 


F(j 0) = A(o) +7 Blo) =n (1 +7 @7,) 


ZAMP IV/16 
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Oly | IE ih | oT,>1 


A, 0 In V2 In @ T, 
dA 1 

E 0 = 1 
du 2 


A und u miissen im gleichen logarithmischen Mass ausgedriickt sein, zum 
Beispiel 


A = 5 In (1 + ® T2) in Neper (N) 
oder 

A = 10 “log (1 + w? T?) in Dezibel (db) , 
dann 

nm =In rh 
oder 


oO @ 
n= 20 Plog = —; 
0 


= 10 (1 Dekade), «=20. 
Mo 
Die Asymptote von A hat demnach die Steilheit von 20 db pro Frequenz- 
dekade, und die dazugehérige Phase ist 90° fiir die Faktoren der einfachen 
Form (1+ # 7,). 
In erster Annaherung kann mit den Streckenziigen die Verstarkungsfunk- 
tion aufgezeichnet werden. 
Die Beitrage 
yA ao Bs) 
1 =i 


In|l1—@w?T?+72¢€Tq@| und arctg 


k6onnen der Figur 7 entnommen werden. 

In Figur 8 wird eine Ubertragungsfunktion G(p) durch die Verstarkungs- 
funktion A und die Phasenfunktion B vollstandig dargestellt. 

Bei der kritischen Frequenz wird A = 0 oder B == -+- 180°. Bei A = 0 muss 
B weniger als 180° betragen oder bei B = + 180° muss A negativ sein. Damit 
im Beispiel (Figur 8) der geschlossene Kreis stabil bleibt, darf die Erhéhung 
der Verstarkung nur wenige Dezibel betragen. Der Einfluss der Verstaérkungs- 
anderung ist aus dem Diagramm A (@) durch Verschieben der Nulline ersichtlich. 

Aus der Bodeschen Beziehung [Formel (12)] lasst sich die Stabilitatsbedin- 
gung grob auch so ausdriicken, dass im Bereich A =0 dA/du< 2, das heisst 
die Steilheit der Verstarkung kleiner als 40 db pro Frequenzdekade sein soll. 

Die beschriebene Darstellungsart bietet gegeniiber der herk6mmlichen Art 
des Nyquist-Diagramms wesentliche Vorteile, besonders wenn es sich nun 
darum handelt, die Stabilitat oder das dynamische Verhalten einer Regulierung 
mit zusatzlichen Mitteln zu verbessern. 
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Akh ge 


oT, 


€=005 015 05 


= E =e a8 
Verstarkungs- und Phasenmass der Ubertragungsfunktion 1 — (w T)? + j2w T ¢. 


FG o) =In [1—(@T)*+ 7°27, 8) 
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Fig. 8 
Verstarkungs- und Phasenmass der Ubertragungsfunktion G. 


( 
p(1+ PTs) (1+ PT%) + 
mit T,=1, T,;=10, 7,=100; T=0,5 


4. Massnahmen zur Verbesserung der Ubertragungsfunktion 


Durch vorhandene, in der Regulierung zu verwendende Organe ist die Uber- 
tragungsfunktion des offenen Kreises schon weitgehend bestimmt. In vielen 
Fallen k6énnte ohne zusatzliche Massnahmen die Stabilitat nicht gewahrleistet 
oder manchen Bedingungen des Pflichtenheftes nicht entsprochen werden. Die 
folgenden Methoden, einzeln oder in Kombination angewandt, erlauben ent- 
scheidende Verbesserungen. 


a) Korrektur durch Zuschaltung von Gliedern H zur Kette G, (Figur 9) 


Die gesamte Ubertragungsfunktion G des offenen Kreises ist das Produkt 
G, H. In den Figuren 10 bis 12 sind verschiedene Méglichkeiten zusammen- 
gestellt, wie eine gegebene Funktion G,, welche allein eine unstabile Regulierung 
ergabe, so verbessert werden kann, dass der geschlossene Kreis stabil wird. 

Am giinstigsten wirken Glieder mit kombinierter Vor- und Nacheilung, die 
auch keine weitern Verstarker nétig machen. 


Fig. 9 
Stabilitatsverbesserung durch Kettenschaltung von Korrekturgliedern zum G-Kreis. 
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1000 14+ 6p 
G, = pa See EN 
‘> ( + 100 4) (1 = fp)?’ “a 1+ 100p 
preys 1+ pif, T, = Rk, C, 
2 Lot Pt; T= (A, + R,) 
Fig. 10 


Wirkung eines Korrekturgliedes mit Phasennacheilung. 


Fig. 11 


Wirkung eines Korrekturgliedes mit Phasenvoreilung. 
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(1 + 3,16 p) (1+ P) 
(1 + 31,6 p) (1+ 0,1 p) 


(1+ p 7,) (1+ pT) 


7 T, 7, p? + (+ I+ Tye) P+ 1 
Cy | T, = RG, T, = R, Cy, T2= RO, 
eee ae eet PE ake Bag iy poke pe whos 
Fig. 12 


Wirkung eines Korrekturgliedes mit kombinierter Phasennach- und -voreilung. 


b) Beeinflussung der Zeitkonstanten bet vorhandenen Elementen durch Einfiihrung 
von Teilgegenkopplungen 


Ein Element habe die Ubertragungsfunktion 
= I 
G, =, ee we . 
Diesem Element wird eine Gegenkopplung mit dem konstanten Faktor H, zu- 
gefiigt. 

Die neue Ubertragungsfunktion G{ des Elementes lautet dann 


Go een me 
I 1+p7T, 1+8,G, 1+7T7,+H,8, | 


(16) 


ist, 2, > 0,so'wird' 7, <= Tpmndu, < hy. 
Die Verkleinerung des Faktors K, muss unter Umstanden in einem andern 
Element kompensiert werden kénnen, damit die Regulierung dem Pflichtenheft 
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entspricht. Die Verkleinerung einer der grossen Zeitkonstanten im Nenner der 
gesamten Ubertragungsfunktion des offenen Kreises erlaubt das Verschieben 
eines Gefallsbruches des Streckenzuges nach héhern Frequenzen und damit die 
Verkleinerung der Phase B bei der kritischen Frequenz. Diese selbst nimmt 
ebenfalls zu. Unter Umstanden wird mit dieser Massnahme bereits eine stabile 
Regulierung erreicht. 


c) Korrektur durch eine zusatzliche Gegenkopplung tiber das ganze System 


In Figur 3 bedeute die Ubertragungsfunktion G selbst ein geschlossenes 
System entsprechend der Figur 2 (siehe Figur 13). Sie ist dann zugleich die 
Ubertragungsfunktion des offenen Kreises des vollstandigen Systems. Dann ist 


G 1 1 
amie ca, 41a Oe, Gh’ ca 
; 1 1 
\., G, | S15 das heisst Ic eS are und Cire 
|H,G,|<1, das heisst fen ee und GwG. 
1 


Die Ubertragungsfunktion H, G, soll ihrerseits die Stabilitaétsbedingungen fiir 
den geschlossenen Kreis G erfiillen. In Figur 14 ist ein Beispiel einer solchen 
Stabilisierung dargestellt. In der Praxis wird von dieser Art viel Gebrauch 
gemacht. 


S$ G (p) 


Fig. 13 
Stabilitatsverbesserung durch zusiatzliche Gegenkopplung. 


5. Frequenzgang des geschlossenen Systems 


Nach Bestimmung einer Funktion G(p), welche unsern Stabilitatsbedingun- 
gen fiir den geschlossenen Kreis entspricht, lasst sich nach Ausdruck (4) M(A) 
berechnen. G(f) ist eine rationale Funktion 


Gtp) = ane) (18) 
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3,16 p? ; 
= 4 F 9. 
A, (1+ 3,166) (1+ 0,1p)’ G, siehe Figur 
WT KT Tp . 
* 14 7(4 4+ fd +"? FG 
Kt, tp" 


(1 + 1 p) (1 + TP) 

T=R,C, T= R,C, Ty = RC, 
Fig. 14 
Wirkung der zusatzlichen Gegenkopplung. 


wobei P,(p) und Q,,(f) Polynome vom Grade » und m in # bedeuten. » < m 


mee | 
MAND) = Bat Oa Bl (19) 


P,(p) und Q,,(p) liegen meist als Produktsummen vor in der Form 


P,(p) = KOT +1) (6 +1) ++ (6 Tr + VY 


und 
Qm(P) a (p Ty a 1) (p Te + 1) ne (p Tm =e 1) , 


wobei Paare von konjugiert komplexen Zeitkonstanten vor allem im Polynom 
des Nenners auftreten kénnen. Will man den Frequenzgang des geschlossenen 
Systems genau berechnen, so muss man die Lésungen des neuen Polynoms vom 
m-ten Grad P,(p) + Q,,(f) = 0 bestimmen, was unter Umstanden eine grosse 
Rechenarbeit erfordert. Sofern |G(j w)|>>1, dann ist M(j @) = 1, das heisst, 
Gefallsbriiche des Streckenzuges A = 1In|G(j @)|, welche in diesem Gebiet 
liegen, haben auf den geschlossenen Kreis einen sehr kleinen Einfluss, der im 
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allgemeinen vernachlassigt werden kann. Zeitkonstanten T des Polynoms P, 
die in diesem Gebiete liegen, miissen Zeitkonstanten 1’ des Polynoms P + Q =0 
entsprechen (t’ ~ T). Im Gebiet, wo | G(j@)|<1, ist M(j w) » G(j w), das 
heisst, diejenigen Faktoren von G, deren Zeitkonstanten in diesem Gebiete 
liegen, treten auch in der Funktion M des geschlossenen Systems auf. Damit 
k6énnen in beiden Fallen einige der Lésungen des Polynoms P + Q = 0 ange- 
nahert angegeben werden, worauf mit bekannten Verfahren die genauen Lé- 
sungen eruiert und dann das Polynom im Grade stark reduziert werden kann. 
Die restlichen und wichtigsten Lésungen lassen sich nun leicht gewinnen. Wich- 
tig sind sie deshalb, weil sie den Frequenzgang des geschlossenen Systems in 
der Umgebung der kritischen Frequenz bestimmen (A ~ 0 und y ~ 0). Auf 
graphischem Wege bedient man sich der Nichols Charts[5], aus welcher sich 
bei gegebener Verstarkungsfunktion A(w) = In| G(j@) | und des zugehérigen 
Phasenspielraumes y die Werte fiir In| M(j@) | und der Winkel von M(j @) 
bestimmen lassen. 


6. Einschwingvorgang des geschlossenen Kreises 


Aus dem Frequenzgang M(p) lasst sich der Einschwingvorgang bei einer 
plétzlichen Anderung des Sollwertes R durch Anwendung der inversen Laplace- 


Transformation berechnen. ° 
r(t) = Ry), RO) =>". (20) 
C(p) = R(p) M(p) ;_ c(t) = L-*[R(p) M(A)], (21) 
Ry M(p) _ Ry F,(P) 
O Raman p ite On) ee, ON oa 


Aus 0+ P=S und der Lésungen p, = —1/t, fiir S(p) =0 findet man 
mit der Partialbruchzerlegung und Riicktransformation 


Ry P(0) 7 Ry P(p,) =t/t’, 
basa tay (3) Oats (23) 
dp /p, 


vorausgesetzt, dass das System zur Zeit t = 0 in Ruhe ist. Die wichtigsten Bei- 
trage zum Einschwingvorgang liefern die Lésungen /, nahe an der kritischen 
Frequenz. Es sind dies aperiodische und unter Umstanden periodisch exponen- 
tiell abklingende Funktionen. Sofern der Frequenzgang M(j w) eine Uberhé- 
hung in der Nahe der kritischen Frequenz aufweist, so hat der Einschwingvor- 
gang sicher periodischen Charakter, und die regulierte Grosse zeigt deshalb ein 
Uberschwingen. Sehr oft werden in dieser Hinsicht Anforderungen gestellt in 
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Verbindung mit kurzer Einstellzeit und Genauigkeit, welche sich mit den gege- 
benen Elementen nicht erfiillen lassen. Da sich der Frequenzgang eines offenen 
Kreises relativ leicht iiberblicken asst, sollten die Forderungen an den Ein- 
schwingvorgang eines geschlossenen Systems direkt in Forderungen an den 
Frequenzgang des offenen Kreises iibersetzt werden kénnen. CHESTNUT und 
Maver haben dies im letzten Kapitel ihres Buches [1] mit Erfolg durchgefiihrt. 
Allerdings ist das Verfahren nicht sehr genau, da es sich nur einer graphischen 
Darstellung vieler berechneter Zuordnungen bedient. Fiir den Berechner eines 
Reguliersystems ist es aber schon sehr wertvoll, einen angenaherten Frequenz- 
gang des offenen Kreises als vorlaufiges Ziel angeben zu kénnen, wozu die 
genannten Kurvenkarten eine wertvolle Hilfe bedeuten. 


7. Beriicksichtigung der Stérungen 


Der Einfluss der Storungen auf die regulierte Grésse (Istwert) soll méglichst 
klein sein. Die Rechnung folgt dem bisher gezeigten Gang, wobei jedoch diejeni- 
gen Teile, die im Gegenkopplungsweg H, des geschlossenen Kreises legen, 
meist das Messwerk, den Verstarker und einen Teil des regulierten Systems um- 
fassen (S = Stérungsgrésse, C = regulierte Grdsse) 


G(p) 
S(O) = SF) Se) ee (249 


Ist keine Regulierung, also auch kein geschlossener Kreis vorhanden, so ware 
C(p) = S(P) Gil(6) . (25) 


Bei geschlossenem Kreis und unterhalb der kritischen Frequenz ist |G, H,|>1. 
Insbesondere bei p > 0, das heisst bei einer plétzlichen Anderung der Stor- 
grosse nach Abklingen des Einschwingvorganges, wird der Faktor 


1 


i+ G0) Ho) | <1- 


Beim Typ 0 bleibt eine Verschiebung der regulierten Grésse bestehen, bei 
den Typen 1 und 2 gibt es keine dauernde Veranderung von C, weil G, H, bei 
p > 0 gegen Unendlich strebt. 

(24) kann auch geschrieben werden: 

eee H,G, M 
On Es eRe as Pa, (24a) 


Da sowohl M (siehe Abschnitt 4) als auch H, rationale Funktionen von p 
sind, ist es ohne weiteres méglich, den Frequenzgang anzugeben und den Ein- 
schwingvorgang bei vorgegebener zeitlicher Stérungsfunktion s(¢) entweder 
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durch das analoge Faltungsintegral zu (5) oder via inverse Laplace-Transfor- 

mation, wenn S(p) bekannt ist, zu berechnen. 
Als Stérungen treten insbesondere auf: 

a) Lastschwankungen; zum Beispiel Veranderung des Lastdrehmomentes bei 
regulierter Tourenzahl. 

b) Schwankungen von Gréssen des Energiereservoirs; zum Beispiel bei elek- 
trischen Antrieben, Netzspannungsschwankungen ; bei thermischen Motoren 
Druckschwankungen usf. 

c) Veranderung des Verstaérkungsgrades eines Teils des geschlossenen Kreises 
infolge Speisespannungsschwankungen, Alterung usf. 

In den ersten beiden Fallen kénnen es sehr rasch veranderliche Vorginge 
sein, wahrend es sich im Fall c) meist um zeitlich langsame Ablaufe handelt. 

Veranderungen der Ubertragungsfaktoren nach c) werden wie folgt be- 
handelt. Nach Figur 2 ist 


E Gy 
Ce Re (25) 
Veranderungen von G, und H, haben zur Folge: 
om 1 Shree ol tN gaeF 
C nie ess © Gy t+ G) A, Bak (26) 


Fir £ > 0 wird 
1 ebnes 
pate ies + x1. 
Ton <1, dagegen it G 


Da sich der zweite Summand 6H,/H, voll auswirkt, miissen an das Messwerk 
H, hohe Anforderungen in bezug auf seine zeitliche Konstanz gestellt werden. 
Dagegen tiben blosse Amplitudenanderungen von G, kaum einen Einfluss auf 
das statische, aber unter Umstanden einen erheblichen auf das dynamische 
Verhalten der Regulierung aus. 


8. Zusammengesetzte Reguliersysteme 


Ein Reguliersystem kann aus verschiedenen Teilreguliersystemen zusam- 
mengesetzt sein, wie zum Beispiel in Figur 15a gezeigt wird. Die Berechnung 
erfolgt dann sukzessive von innen nach aussen. 


O= Rape RG, (27) 
CBG RCS (28) 
C =Ravat = Gt, (29) 
ee ew abelg C= GG (30) 
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Eine genaue Berechnung ware ausserordentlich umstandlich, deshalb macht 
man zunachst Gebrauch von den gebrochenen Streckenziigen in der Bodeschen 
Darstellung und der vereinfachten Bestimmung des Frequenzganges der ge- 
schlossenen Teilsysteme nach Abschnitt 4. Das Verhalten der beiden Schal- 
tungen von Figur 15a und d ist nach aussen identisch, da eine Vertauschung der 


Fig. 15 


Schema einer komplizierteren Regulierung. 


Reihenfolge von G; mit G, fiir die Berechnung keine Rolle spielt. G{ ist bereits 
ein Servosystem mit stark verbesserten Eigenschaften. Durch Hinzufiigen des 
Verstarkers G, mit einer angepassten Ubertragungscharakteristik kénnen Eigen- 
schaften verwirklicht werden, die ohne eine solche Schaltung nicht méglich 
waren. Die Approximation in der Berechnung kann dadurch verbessert werden, 
dass zum Beispiel Uberhéhungen in der Nahe der kritischen Frequenz von Gj 
und G{ durch quadratische Formen 


1+2&pT+T?p? 


im Nenner beriicksichtigt werden. Eine Kontrolle fiir einige Frequenzwerte um 
die kritische Frequenz muss unter Verwendung der Nichols-Diagramme durch- 
gefiihrt werden. Zur Bestimmung der Einschwingvorgange geniigt meist ein 
angenadherter Ausdruck fiir das gesamte System. 
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Z 
Cas. De lee 
6—— 6 Se Oa NN \ 
107! ; SN Se 
| SQ AK 
(e) 
N 
ASS Ne 
| IS Ne 
ee 
Xe 


Fig. 16 


Abschatzung der Wirkungsweise der Regulierung von Figur 15. 


teal Gs 
1 mite 1000 OO SESS 
H, 1 “(1+ 100 p) (1 + p)? * 1+ HG, 
eae 15 (1 353,16 p) (14-)0)1 P) cn. as 
. E Gs H, 3,16 p2 1 1+ Gi G, 
> aS Ose OTE z 10 (1 + 100 p) 
2 Gy Hs + sL(lb 1000 p) 
—oO—o— Gj 
1 


1. Approximation Gj ae 6) (1 0,316 p)2 


: ” 1 9 J 
2. Approximation Gj G+ Pp) Ut 2eT per py 1 == OSM, <== 0) 1 


4 ; - (1 + 3,16 p) ‘ | . 
3. Approximation Gy, (enteral £5) er 26, Ip + 13 P®) 


T,= 1,86, £= 0,88, T,= 0,301, & =°0,306 


I 
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9. Bestimmung der Ubertragungsfunktion eines Elementes 


Beispiel. Tourenzahlregulierung eines Gleichstromnebenschlussmotors mit 
Hilfe des Felderregerstromes (Figur 17a). Als Stérungsfunktionen treten auf: 
a) das Lastdrehmoment M,; b) die Batteriespannung am Anker U9. 


U, M,; 


Fig. 17 


Ersatzschema eines Nebenschluss-Gleichstrommotors. 


Statische Verhdltnisse : 
Fluss der Erregerwicklung : 


@ = K, 1, (Ankerriickwirkung vernachlissigt), 
Motordrehmoment M, (= Lastdrehmoment M,): 
M,=—=M,=-K,1,@, 
Klemmenspannung am Anker: 
U,=1,R+K, On, 
daraus die Tourenzahl 


al 1 U, M,R 


K, eRihyta) Aya erat 

Die Konstanten K,, Ky und K; sind von J,, Jy, 7 und M, abhangig und kén- 
nen nur bereichweise als konstant angenommen werden. Die Tourenzahl ist 
eine nichtlineare Funktion der Eingangsgrésse J,. Auch hier ist nur eine be- 
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schrankte Variation des Erregerstromes J, zuldssig, damit das Problem lineari- 
siert werden kann. 


Abkiirzung: 
R 
Ke Keke Ss 
nis U, M, . . . 
n Rat K P Gleichgewichtsbedingung. (31) 


Kleine Stérungen aus der Gleichgewichtslage: 


Tie aah b18) |e n=n, [1+ ¢,(é)], | 
U, = Ug [1+ ego(é)], Ig = Ig [14+ &0()], (32) 
M, = My [1+ eyo(t)], My = Myo [1 + e,,(8)]. | 
Gleichungen: 
}=K,1,. (33) 
Elektrischer Kreis: 
Cee Ll. Ke re (34) 
Bewegung der Achse: 
22nj,=M,—M,, (35) 


J, totales polares Tragheitsmoment bezogen auf die Motorachse, 
M, SKyd, Oe hy Ky dls (36) 
Werden die e(t) als rein sinusoidale Schwingungen der Frequenz w = —]7p 


gewahlt und bedeuten E;,4-, die zugehérigen Zeiger, so ergeben sich unter 
Beriicksichtigung der Gleichgewichtsbedingungen die Gleichungen 


Uy9 Eye = Lo Eig (R+ 6 LYSE Ke Ky hig) (Bey Es (34a) 
M4 Ey, =2n? J, U) E,, + Mg Ems» (35a) 
My Ey = K, Kg Lo9 Lio (Ei + E;9) - (36a) 


Die Doppelprodukte £;, E,, und E£;, Ej, werden vernachlassigt gegentiber den 
linearen Gliedern. Man erhalt als Ausgangsgrésse bzw. deren relativer Wert E,, 
in Funktion der relativen Werte des Erregerstromes £;. und der Storungen Ey2 
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und Ey. den Ausdruck (37) 


E,, K’ (p+ 1) — Ema Moy (1 + Ty) — Evo K" 


ee K" pT, (PT, +) +1) 


. 


rr ree (Eppa oe gor 1Epdarg eel a 
K' = Ky Kal le aie oie ae Be 
Wa LC PAM feed a 
Ae ee wal oe Steet te ma JD 
Jae —- K ) o R , iP I, ’ 
id ifs Ky Ky Tyo Ls ; 
TR Ky Ky Tyo L290 — (Tin %o/®) 
In der Figur 17c bedeutet nun: 
ECS Kk" 
G=zr Gutl), R= zr, | 
M. 1 as 
Ase, cot: d . 
= K" (1+ pT) , Gy [p Ts (p T, + 1) 2b 1] | 


Alle diese Ubertragungsfunktionen sind dimensionslos und beziehen sich nur 
auf den Zusammenhang der relativen Anderungen der Eingangsgrésse J,, der 
Stérungsgréssen U, und M, mit derjenigen der Ausgangsgrésse 1. 

Daneben muss die statische Gleichgewichtsbedingung (31) erfiillt sein. Soll 
die Regulierung ein breites Gebiet umfassen, dann ist es erforderlich, die obigen 
Funktionen fiir alle in Frage kommenden Gleichgewichtslagen neu zu bestim- 
men. Die zusatzlichen Verstairker, Tachometerdynamo usf. miissen dann im 
ganzen verlangten Regulierbereich solche Ubertragungsfunktionen ergeben, 
dass der geschlossene Kreis in allen Fallen stabil ist. 


10. Einfluss einer endlichen Laufzeit im offenen Kreis 


Eine endliche Laufzeit t kann in einer Ubertragungsfunktion durch den 
Faktor exp(—p t) beriicksichtigt werden. Die bisher beschriebenen Methoden 
miissen mit Vorsicht angewandt werden, da die im Abschnitt 2 vorausgesetzten 
Bedingungen nicht mehr erfiillt sind. KGPFMULLER behandelt in verschiedenen 
Ver6ffentlichungen [6] diesen Sonderfall, der sehr oft in der Praxis auftritt 
(zum Beispiel Schalldruckregelung in einem freien Schallfeld; Leistungsregu- 
lierung von Generatoren auf Grund ferniibertragener Messwerte). Im Gegensatz 
zu dem in den vorigen Abschnitten Gesagten kann sich eine Stérung auf die 
Ausgangsgrésse zunachst voll auswirken. Sie muss dann durch die Regulierung 
sukzessive auf den urspriinglichen Wert zuriickgefiihrt werden, wenn sie stabil 
arbeitet. 

In vielen Fallen aber wird der gemeinsame Faktor im Nenner der Ubertra- 
gungsfunktion (Ausgangsgrésse dividiert durch Fehlergrésse) fiir irgendeinen 
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Eingang mindestens eine Zeitkonstante von erheblicher Grésse aufweisen (Trag- 
heit des zu regulierenden Systems). Mit andern Worten, die Ausgangsgrésse 
kann sich nicht sprunghaft andern. Ist das zu regulierende System tragheitslos, 
so muss der Messkreis eine Zeitkonstante aufweisen, damit eine erhebliche Re- 
gulierwirkung (Verminderung der Stéreinfliisse) iiberhaupt méglich ist. An 
Stelle der Betrachtung des Frequenzganges wird anschaulicher mit einer ver- 


4 T 
Fig. 18 Fig. 19 
Beispiel einer Sprungtibergangsfunktion Regulierung, bei der die Ubertragungs- 
eines Systems mit endlicher Laufzeit. funktion A des Leistungsflusses eine 


Funktion der Abweichung ist. 


einfachten Sprungiibergangsfunktion oder deren zeitlicher Ableitung, der Im- 
pulsiibergangsfunktion, gearbeitet. Als einfacher Fall fiir die analytische Be- 
handlung untersucht KipFrMULLER die normierte Sprungiibergangsfunktion des 
offenen Regelkreises 

ei) = 1 =e 9, thy ebi= 0. (39) 


Das Reguliersystem sei nach Figur 19 aufgebaut. A sei eine Funktion von 
E—R—H,C 
C = A S(S Quelle, Energiereservoir) 


66 = 5A, S; Agios 100 = (FF), dE 


fir R konstant: 


dA 
eH, 06 on, OCS Solly fae 8C + A, dS 
und daraus 
oC OSs 1 Os 1 (40) 
pm ee dA So aie Carte aa 
eee eee as ala 


ZAMP IV/17 
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C, (aA bil 
G= 4° (ae), a— KP Ue = K 57> se 
C, (aA 
K = 7° Cal H,(0) . (42) 


In Anwendung fritherer Prinzipien erkennt man leicht die Grenze der Stabilitat 
fiir 


bG | Sil are (G) ae —70.. 
Ki V1 + wt? und arctg(m@pt) + wot, S72. (43) 


Bei K > 1, was beieiner effektiven Regulierung immer der Fall ist, wird wyt » K 
und arctg@,)t ~ 7/2. Es muss dann sein 
Kt, 
t 


VAS 


3° (44) 


damit die Regulierung stabil ist. 
KUPFMULLER lést das Problem mit Hilfe der Integralgleichung (45) und 
diskutiert die verschiedenen Stabilitatsbereiche. 
Se 
nit) = f(t) — K [ elt 0) WC) a (45) 


0 


h(t) gesuchte Funktion, 

/(t) Ausgangsfunktion einer gegebenen Stérungsfunktion bei aufgeschnittenem, 
aber in gleichem Zustand sich befindendem Regelkreis [1 (¢ — 4) gesetzt]. 

Mit (31) und einem Ansatz fiir /(¢) 


h(t) = hy + hy e~*' sin (wp t + @) (46) 
ergeben sich zwei Beziehungen: 


1 Mgt , Tt 
Toe CROCE rr 


e~ HIT — Si(qg ty) ePetil(t ooh) | (47) 


Entsprechend den periodischen Funktionen si(@, ¢;) und tg(@, ¢,), gibt es viele 
Wertepaare « und @ , von denen aber nur das erste fiir den Einschwingvorgang 
wichtig ist. Die gleichen Beziehungen kénnen auch leicht gewonnen werden aus 
der Lésung der charakteristischen Gleichung 1 + G(p) = 0, welche unter Ver- 
wendung von (41) lautet: 


pt+1+Ke*?=0. (48) 
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Mit den Lésungen 4,,, = —a + 7 w, eingesetzt, erhalt man 


tg wg =——°"_ und Ke%*=Vagr? + (1—at)?, (49) 


t— 1 
welche, etwas umgeformt, die Beziehungen (47) ergeben. 
Der aperiodische Grenzfall ist bei w, = 0, das heisst 


Tt 
Ki, 


K e**=at—1 oder aus (47) e Hlt — (50) 


wobei infolge ¢,/7 <1 


t 


—t,/t 
é = 
1 und Ki, 


=€. 


Dies in (50) eingesetzt und unter Beriicksichtigung, dass K e*“ 1, also 
auch « t >1, erhalt man e** = « ¢, e und daraus 


Org = las 


Trotz der sehr grossen Zeitkonstante t in der Ubertragungsfunktion wird die 
Ausgangsgrosse rasch auf die bleibende Abweichung zuriickgefiihrt. 
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Résumé 


Se basant sur de nombreuses publications, on expose et on démontre les 
méthodes les mieux appropriées 4 l’analyse des problemes de réglage. L’article 
traite les différents moyens pour améliorer lastabilité, l’influence de perturbations, 
les effets transitoires, les syst¢mes de réglage composés, l’influence d’une vitesse 
de propagation finie dans le systéme et la détermination de la fonction de transfert 
a l’exemple d’un moteur en dérivation. Les méthodes décrites ne sont valables 
que pour des systémes linéaires. Des systémes non linéaires peuvent aussi étre 
étudiés par la méme méthode, pour autant qu’il s’agisse de petites perturbations 


autour d’une position d’équilibre. 
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Summary 


Based on the numerous publications dealing with feedback-control systems, 
the essential features of the phase-attenuation method most suitable for the 
design of regulators and servo-mechanisms are shown. Measures to increase 
stability, series and parallel compensation, influence of disturbances, transient 
behaviour and the effect of a finite time delay within the system are discussed. 
As an example, the fransfer function of a DC shunt motor is derived. The method 
is suitable for linear systems only, nonlinear systems may be analyzed by assum- 
ing small deviations from a steady state value. 


(Eingegangen: 9. Mai 1953.) 


Characteristic Surfaces of the Equations of Motion 
for Non-Newtonian Fluids’) 


By JERALD LAVERNE ERICKSEN, Bloomington, Indiana, U.S.A.?) 


1. Introduction 


In phenomenological theories of continuum mechanics, characteristic curves 
or surfaces, when they exist, are usually of some importance. In the classical 
linear theory of elasticity, waves of compression and distortion are such sur- 
faces [1]*). In gas dynamics and in plasticity, one frequently makes use of the 
properties of characteristics in solving problems, and physical interpretations 
for them have been suggested [2], [3]. There is thus some reason to believe that 
such surfaces may be of importance in modern theories of fluid dynamics such 
as those proposed by REINER([4], RIVLIN[5], and TRUESDELL[6]. We consider 
RIVLIN’s equation of motion for non-Newtonian fluids. These are 


Hitep=o(S +ui, v'), (1.1) 
ti = —p 6; + G,(II, II) di + G,(II, IT) di, d?, (1. 2) 
v,=di=0, (la) 


where ¢ is the stress tensor, /¢ the body force per unit mass, v? the velocity 
vector, d,; = (v;,; + v;,;)/2 the rate of deformation tensor, II = — (di d3)/2, 
III = det. dj, p is an arbitrary hydrostatic pressure, and 9 is the density, as- 
sumed constant. It is assumed that, in the range of values of II and III con- 


1) Prepared under Army Contract DA-33-008 ORD-454 with Indiana University. 


*) Graduate Institute for Applied Mathematics, Indiana University. On leave of absence from 
the Naval Research Laboratory, Washington, D.C. 


3) Numbers in brackets refer to the Bibliography on page 266. 
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sidered, G, and G, are single valued, continuously differentiable, and that not 
both of these functions vanish identically. For these equations, the characteristic 
condition is obtained. For a plane motion, it is found that, at a point where 


G=11G,=0, 
or 
II 0G, 
fete cay eres 


every direction parallel to the plane of the flow is a characteristic direction. At 
a point where neither of these conditions is satisfied, there are at most six 
distinct characteristic directions parallel to the plane of the flow, and explicit 
determinations for these are given. Certain inequalities must be satisfied in 
order that these exist. Similar results are obtained for a slightly more general 
class of motions. 

Setting G, = k (—II)~ 1”, G, = 0, where & is a positive constant, one obtains 
the von Mises equations of plasticity. This case was discussed recently by 
Tuomas [7]. 


2. An Alternative Form of the Equations of Motion 


Using (1.2) to eliminate ¢*/ from (1.1), one obtains 


—pit eft Gait G (diat) + (at 2 + ah ahi 32) am, wi 
o foes or + 484" ur) dn, a (FF + vij0/). 
Using (1. 3) and the relations 
one di, ie he aie 
one can put this equation in the form 
—pitot+ Tho" —0 (Se +0404), eis 
where 
Ti Gy, s On ig aee (Of Un + a dj, + Om Fe) 
= ayy (0°) Soe + 4” arr) 2.2) 
+ in Gy (a sr arr tant eok 
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In all of which follows, it is assumed that o, /;, G, and G, have been specified. 
Then (1. 3) and (2. 1) yield four equations for determining the four unknowns # 
and v'. Given # and v’ satisfying these equations, one can obtain #4 from (1. 2). 
The stresses thus determined satisfy (1. 1). Therefore we need consider only 
(1. 3) and (2. 1). 


3. The Characteristic Condition 


In obtaining the characteristic condition, we consider only stationary flows. 
However, it is easily verified that the same condition holds for nonstationary 
flow. 

In this section, Latin indices take on the values 1, 2, 3, Greek the values 
2,3. Let S be a surface defined by an equation of the form /(x1, x?, x8) = 0, 
where f is continuously differentiable and f:* /,; > 0 at each point of S. If P is 
an arbitrary point on S, one can chose rectangular Cartesian coordinates such 
that, in the neighborhood of P, S is represented by an equation of the form 


a) = p(x", x*) , (3. 1) 


where g , = 0 at P. In this coordinate system, the unit normal vy, to S has com- 
ponents v? = + df at P. 

Consider the function # on S. Using (3. 1), one can express # as a function 
of the surface coordinates x*. Then 


da 
a 3 8,0» (3. 2) 


where dp/dx* denotes surface differentiation of p. Similarly, 


dv; 


Ve = “axe Vib ae» (3. 3) 
A dv;,; 
Oa aa i Cae (3. 4) 
It follows from (3. 4) that?) 
= dvjy 
Opal ae oe ol Gos (3. 5) 
dv; p dv; : 
UY Pode pos os per ae oO Hg ee (3. 6) 


Differentiating (1. 3) with respect to x1 and using (3. 5), one obtains 


0, 4, = v* Mhad Sm, 
1,11 10 8 gt yee (3. 7) 


1) A derivation of these equations and a general discussion of characteristics is given in 
HADAMARD [1]. 
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If the functions , v,, and v,, ; are assigned on S subject to (1.3) and (3. 3), 
equations (2. 1), (3. 2), (3. 5), (3. 6) and (3. 7) yield 21 linear equations in the 
21 unknowns #,; and v;,;,. At a point of S where the determinant D of the 
coefficients in these equations is different from zero, p, ; and v;, ;, are determined 
uniquely. If D = 0, they are not. If, at each point of S, D -= 0, S is said to be 
characteristic relative to the assigned data. The equation D = 0 will be called the 
characteristic condition. 

Consider the 21 equations at P, where g , = 0. The derivatives P ,, De b> 
U;,41 and v,,,, are then given directly by (3.2), (3.5), 3.6), and (3.7). It then 
follows from (2. 2) that the quantities 


21 21 31 31 11 rb 
TH Usa St san, FSi Py + 134 3,41, —P + TH e+ TH 3,1. 


can be determined at P. The derivatives p,; and v,,;; are thus determined 


uniquely unless 
TH TH — TH TB =O. (3. 8) 


Clearly, (3. 8) is the characteristic condition, expressed in this special coordinate 
system. To express (3.8) in invariant form, it is convenient to introduce a scalar 


U=6,+G,"', 


and a tensor 


: : . i i m 0G 0G 
Us = Ge (di — vw) + 2 (ui = 9 oy") (x dni eT or) 


+ 1m i m n 0G 0G 
+ 2 (4, we — ¥ me”) (1 dy Sate — Pi are) 


where p= divi. We note that », Ui = 0. It follows from (2. 2) that, at P, 
273, = U 65+ Uj. 
Thus (3. 8) can be written as 
U2 + U U% + U2 U3 — UZ UZ =0. (3. 9) 


Furthermore U! = + », Uj = 0, so that (3. 9) can be written as 


(Uj)? — UZ Ui 
Z 


U2+U Ui+ 0) (3. 10) 
which is an invariant form of the characteristic condition. A characteristic 
direction is the direction of a unit vector y, satisfying (3. 10). The number of 
characteristic directions corresponding to a given set of functions Gy, G,, and 
d,; may be finite or infinite, depending on the nature of these functions. 
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In certain theories of plasticity and fluid dynamics, it is assumed that 
G, > 0,G,=0 Then U= G,, Uj = A‘ B;, where 

0G, 0G, 
ms“orir 4 O11’ 


At=2(ut—vi mr), Bo=pnd 
thus 
(Ui)? — Ui Ui = (A’ B)?-— AB A! B.=0. 
Equation (3. 10) then becomes 
G,+ At B,=0. 


Now consider a flow in which there exists a surface across which # and v,, ; 
are continuous, but on which there is a finite discontinuity in at least on of the 
quantities #, ; and v,,;;. Such a surface must be characteristic. However, a char- 
acteristic surface need not be a surface of discontinuity. 


4. Plane Flow 


Introduce an orthogonal coordinate system in which the surfaces x3 = const 
are parallel planes, and consider the problem of determining the characteristic 
direction v,; for the case where 


Ut, Shan 8; Soar ao ark a 
at each point of the flow. Then 
= 1 2 42 
d;, = d, &; €; + dy ee; , 


where e} and e? are mutually orthogonal unit eigenvectors of di and d,, dy are 
the corresponding eigenvalues. Using (1. 3), one obtains 


d;, = d, (e; & — ee). 
Thus, if m is the angle determined by e} and »,, 
d,,v' = d, [(et v')? — (ef v')"] = d, (cos*g — sin?y) = d,cos2p. (4.1) 
Since vy; Uj = 0, the determinant | U}| must vanish. Thus 
(U)2= Uj U=2|Uj1=0. 
From (3. 10), either 


U+vUi=0, (4, 2) 
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or 
UE 0: (4. 3) 


Let P be a point where (4. 2) holds. Introduce rectangular Cartesian coordinates 
such that, at P, »; = 6%. Then 


dy; uw a Ay 5 di = dye (dy, + don) = 0. 


Using this result it is easy to show that 


U+Uui=U+ UZ = G, — 2 (dy)® =F = 0- 
But 
— (yp)? = dy dog — (dy)? + (dep)? = dy dg + (d;; v¢ v4)? = —d? + (d;, vé v*)?, 
and II = —d?. Thus (4. 2) can be written in invariant form as 


Gy + 2F dy 9* 4)? + 11] 2 <0, 


or, using (4. 1), 
0G, 
oll 


Greil sin 29 = 0. (4. 4) 
We exclude the exceptional case G, = I] 0G,/0II = 0, which is mentioned in 
the introduction. Then (4. 4) can be satisfied if and only if 
0G,\-1 

oe say on (2 I a) Ee (4. 5) 
When (4. 5) is satisfied, (4. 4) determines four characteristic directions. These 
are all distinct if the strict inequalities hold in (4. 5). If w is an angle satisfying 
(4. 4), two make angles m and —q with e} while the other two make angles » 
and — with e?. A necessary condition that (4.5) hold for «<IL< 6 < 0, 
where G,(f, 0) + 0, is that 


Glee B \1/2 
7 eye.) = (ar) 
fone. ll p. 
Suppose now that (4. 3) holds. Then 
G,+ G,d,,v'vi=0, 


or, using (4. 1), 
G,+ G,d,cos2p=0. (4. 6) 


Excluding the case G, = II G, = 0, one must have 


G 
gia] si (4.7) 
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in order that (4.6) define real characteristic directions. When the strict inequal- 
ity holds in (4.7), (4.6) defines two distinct characteristic directions. The 
vectors e} and e? bisect the angles formed by these two. 


5. Locally Plane Flow 
In any flow one has 


Let 2 2 3 8 
d,, = d, €; €; + dy & & + dye; e , 


where e}, e?, and e? are mutually orthogonal unit eigenvectors of dj and d,, 
dy, dg are the corresponding eigenvalues. A flow for which III = d, d, d; =0 
will be called locally plane. One can, without loss of generality, take ds = 0. 
Then, using (1. 3), 

d;,= d, (ee — ra) e;) ‘ 


For the special case of the von Mises equations of plasticity, it turns out that, 
for real characteristic directions to exist, the motion must be locally plane. 
Furthermore, 


ye=0. (5. 1) 


Consider the problem of determining the characteristic directions for such a 
flow for the case where (5. 1) holds. If one introduces a rectangular Cartesian 
coordinate system such that, at an arbitrarily preassigned point P in the flow, 
e? = 03, one sees that this problem is formally identical with that just solved. 
These characteristic directions are thus given by 


: bGs irs 
G,+ 21I 5 sin’29=0, G+ G,d,cos2p=0, 
where cosy = e} »*. In such flows, the characteristic surface elements will not 
necessarily unite to form surfaces [7]. 

RIVLIN’s general solution for Poiseuille flow provides an example of a flow 
which is locally plane but not plane [8]. This flow is axially symmetric and 


(5. 1) is the condition that the characteristic surfaces be surfaces of revolution 
about the axis of symmetry. 
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Résumé 


Aprés avoir dérivé la condition caractéristique pour les équations du mouve- 
ment des fluides «non newtoniens», nous avons déterminé les directions caracté- 
ristiques pour le cas du mouvement plan. 


(Received: February 2, 1953.) 


Luftkrafte an einem schwingenden Schaufelkranz 
kleiner Teilung 


Von HEINz SOHNGEN, Darmstadt?) 


1. Anlass der Untersuchung und Zusammenfassung 


Die vorliegende Arbeit befasst sich mit rein theoretischen Uberlegungen 
zum Problem des Schwingens von Schaufelkranzen in axial durchstrémten 
Maschinen. Zur Vereinfachung nehmen wir dabei an, dass der stehende oder 
umlaufende Kranz sich in einem Ringkanal geringer Hohe (ebenes Problem) 
befinde, der sich nach beiden Seiten ins Unendliche erstreckt. Das str6mende 
Medium mége inkompressibel sein, Stérkérper mégen sich in der Nahe des 
Kranzes nicht befinden. 

Um das Schwingungsverhalten eines derartigen Kranzes beurteilen zu k6n- 
nen, ist die Kenntnis der Luftkrafte, die an den schwingenden Schaufeln an- 
greifen, unerlasslich. Aussagen tiber Luftkrafte sind aber bisher nur fiir den 
Einzelfliigel bekannt, und zwar unter den folgenden Voraussetzungen: 


a) diinne und schwachgewélbte Profile (Streckenprofile), 

b) kleine Auftriebsbeiwerte, 

c) kleine Schwingungsamplituden. 

Das strémende Medium wird dabei als reibungsfrei (siehe zum Beispiel [1], [2] 
und [3])2) und im allgemeinen auch inkompressibel [4] angenommen. Die Er- 


1) Technische Hochschule Darmstadt. 
2) Die Ziffernin eckigen Klammern verweisen auf das Literaturverzeichnis auf Seite 297. 
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gebnisse dieser Theorie wird man anwenden kénnen, wenn das Teilungsverhalt- 
nis a/c (siehe Figur 1) gentigend gross ist. Ist das Teilungsverhaltnis aber klein, 
so muss erwartet werden, dass infolge der gegenseitigen Induktion der Schaufeln 
Abweichungen eintreten, und zwar um so starker, je kleiner a/c und je grosser 


ft 


Fig. 1 
Gitter in der Abwicklungsebene. 


der Auftriebsbeiwert c, ist. Aus diesem Grunde sollen die Luftkrafte an einem 
schwingenden Gitter unter den folgenden Voraussetzungen berechnet werden: 
a) ebenes Streckengitter kleiner Teilung, das heisst lineare Theorie in a/c, 

b) inkompressibles und reibungsfreies Medium, 

c) kleine Schwingungsamplituden. 


Der Tatsache, dass es sich um ein Gitter aus einem Schaufelkranz handeln soll, 
tragen wir dadurch Rechnung, dass Schaufeln des Gitters mit Nummern, die 
sich nur um Vielfache der Blattzahl 3 unterscheiden, gleiches Verhalten auf- 
weisen sollen. 

Unter diesen Voraussetzungen wird ein Verfahren angegeben, das bei be- 
liebig vorgegebener Schwingungsform des Schaufelkranzes sowohl in Tiefen- 
richtung als auch tiber dem Umfang die Berechnung der zugehérigen Luftkrafte 
gestattet. Explizit wird die Rechnung aber nur fiir die praktisch hauptsachlich 
interessierenden Falle der Drehung und Biegung durchgefiihrt, und zwar wer- 
den hierfiir die an den Schaufeln angreifenden Krafte und das Moment um die 
Vorderkante angeschrieben. Diese reichen aus, um das Eigenschwingungsver- 
halten, das heisst die kritische Geschwindigkeit eines Schaufelkranzes, zu be- 
stimmen. Um auch das Verhalten eines Schaufelkranzes bei Fremderregung 
beurteilen zu kénnen, bendtigt man ausser den Kriaften fiir die genannten Frei- 
heitsgrade auch noch diejenigen, die sich an den Schaufeln bei periodischer 
Stérung des Zustroms einstellen. Diese sollen spiter angegeben werden. 

Es kann nicht erwartet werden, dass die Luftkrafte an einem schwingenden 
Gitter eine einfache Form haben, dazu ist die Anzahl der in das Problem ein- 
gehenden Parameter zu gross. 
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_ Abschliessend wird an einem Beispiel gezeigt, wie mit den angegebenen 
Kraften die kritische Geschwindigkeit berechnet werden kann, insbesondere 
wird fiir den Fall der Biegeschwingung eines Kranzes mit dem Staffelungs- 
winkel B = 0 eine explizite Formel dafiir abgeleitet. 


2. Liste der Bezeichnungen 


= Schaufeltiefe, 

= Teilung, 

Staffelungswinkel, 

= Blattzahl des Schaufelkranzes, 

= 22C/a 3, 

= Q etB 

= Blattnummer, 

dimensionslose Tiefenkoordinate (€ = 1 entspricht Profiltiefe), 
= dimensionslose Umfangskoordinate (7 = 1 entspricht Umfang des 
Kranzes), 

= m/3 Umfangskoordinate der m-ten Schaufel, 

rechtwinklge Koordinaten (dimensionsbehaftet), 

= stationdrer Anstellwinkel (weit vor dem Gitter), 

= stationare Anstré6mgeschwindigkeit vor dem Gitter, 
stationdre Abstromgeschwindigkeit, 

= mittlere stationdare Anstr6mgeschwindigkeit der Vorderkante, 
mittlere instationare Anstro6mgeschwindigkeit der Vorderkante, 
=) Luttdichte, 

Druck, 

= Auftriebsbeiwert, 

2 a sina/c cos(a+ f) modifizierter Auftriebsbeiwert, 

6 = stationdre Zirkulation an der Vorderkante, 

W, 2(€, 1, s) = instationare Wirbeldichte, 

W,cG(n,s) = instationare Zirkulation an der Vorderkante, 

v = Frequenz der Schwingung (Hertz), 

t = Zeit, 

S = WK We. 

a) = 22 c/W,, reduzierte Frequenz, 

AA EH nS) Auslenkung des Schaufelpunktes &, 7,, in *-Richtung, 

AY(E, Nm, S) = Auslenkung des Schaufelpunktes €, 7,, in y-Richtung, 

AME, Nm 8) = AXE, tm 5)14 3, 

A,(§, Nim s) = Ay(é, Nim s)/a 3, 


HY BS iv} 
I 


3 


HSl S wy 
= 
I I 


SS REN 
55 
Il 


S 


g 
S 
l| 


i) 


ae Beau) 
I 


ries 
I 


n = Ordnung der Fourier-Transformation (siehe Abschnitt 3), 
Wes ges a, all OY Me Ae, ony 

() = Fourier-Transformierte u-ter Ordnung, 

On = (n/|n}) ©, 

o z= |w| 2. 


3. Vorbemerkungen 


Den Schaufelkranz mit den 3 Schaufeln denken wir uns in eine Ebene ab- 
gewickelt und erhalten so ein Gitter, dessen Schaufeln wir numerieren ..., — 2, 
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—1, 0, 1, 2, .... In dieser Ebene fiihren wir ein rechtwinkliges Koordinaten- 
system x, y ein, dessen x-Achse den Schaufeln parallel ist und dessen Ursprung 
in der Vorderkante der nullten Schaufel liegt. Dieses System, das noch dimen- 
sionsbehaftet sein mége, dient zur Orientierung der Geschwindigkeiten. Mit 
u bezeichnen wir stets eine Geschwindigkeit in x-Richtung und mit v eine 
solche in y-Richtung. 

Zur Orientierung auf den Schaufeln ist es zweckmissig, noch ein zweites 
— bei umlaufenden Kranzen schiefwinkliges — Koordinatensystem &, 7 (Gitter- 
koordinaten) einzufiihren, dessen &-Achse in Richtung der x-Achse weist und 
dessen 7-Achse mit der Gitterachse zusammenfallt. €, 7 seien dimensionslos, 
und zwar entspreche € = 1 der Gittertiefe (c) und 7 = 1 dem Umfang (a 3) 
des Kranzes. Es ist dann 


z=x+tiy=cé&+te azn, (3,1) 
bzw. 


1 ¥ 
=—(x—ytgf), n= a3cosp ° 


Ist f eine Funktion von x, y, so schreiben wir auch in nicht misszuverstehen- 
der Weise (x, v) = f(é, 7). 

Hat das Gitter Streckenprofile, so sind die Koordinaten der einzelnen 
Schaufeln durch 


? 
n=, _ ——, 0, =; oe Rey Oe eee 


gegeben, und wir setzen noch 7, = m/3. 

Da sich das Verhalten der Schaufeln mit der Blattzahl 3 periodisch wieder- 
holt, so miissen sémtliche in dem Problem wesentlichen Funktionen in » die 
Periode 1 haben. Man wird daher Fourier-Reihen bzw. -Polynome einfiihren. 
Wir bedienen uns dabei der komplexen Darstellung. Ist /(7) eine Funktion der 
Periode 1, so schreiben wir 


ee) 


f(n) = py) fe pene (3, 2a) 


-—c 


mit 
1 
~™ =| ta) ei2™nn dn ; 
0 


Handelt es sich um eine Funktion, die nur fiir die diskreten 7,,-Werte definiert 
ist, zum Beispiel um die Schwingungsamplituden der einzelnen Schaufeln, und 
ist f, der Wert, den diese Funktion auf der k-ten Schaufel, also fiir Hp an- 
nimmt, so interpolieren wir diese Funktion durch ein endliches Fourier-Poly- 
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nom und setzen dazu 


n 1 : 1200 m4 
Ee per PIN Ss (3, 2b) 
m=1 
und 
iu 4, wenn 3 gerade, 
3 — 
3 == il 
2» Wenn 3 ungerade. 
Dann ist 
Kh) =e De pee (3, 2c) 
=%+3%4+1 


ein trigonometrisches Polynom, das an den Stellen 17 = 7, die Werte f, an- 
nimmt. 
In der vorliegenden Untersuchung werden des 6fteren Integrale der Form 


F(n) = | Hn) h(n — 7’) dy’ (3,3) 


auftreten, wobei /(7) und g(7) stetige Funktionen der Periode 1 sind. Fiir die 
Fourier-Transformierten von F (7) gilt dann der sogenannte Faltungssatz 


FO — f) 7 | (3,4) 


An Stelle der Zeit fiihren wir die dimensionslose Grésse s = W, t/c ein; 
dabei bedeutet W, die stationare Abstro6mgeschwindigkeit hinter dem Gitter 
und c die Schaufeltiefe. Wir beziehen hier im Gegensatz zur Tragfliigeltheorie 
auf die Geschwindigkeit W, anstatt W, oder W,,, da die freien Wirbel mit der 
Geschwindigkeit W, abschwimmen. Diese drei Geschwindigkeiten sind nicht 
gleich, da wir nicht nur den Fall schwacher Belastung betrachten. Ist y die 
Frequenz in Hertz, mit der die Profile schwingen, so fiihren wir noch die so- 
genannte reduzierte Frequenz 


Qe vy GC 
W, (So) 


ein. 

Die Schwingungen der einzelnen Schaufeln seien vorgegeben und mogen 
mit der reduzierten Frequenz w erfolgen. Die Auslenkung des Schaufelpunktes 
E, Mm in Richtung der x-Achse kann dann in der Form 


AX(E, Nims 8) = 4m(€) cos[@ s + Ym(E)] (3, 6a) 


bzw. die Auslenkung in y-Richtung 
Ay(E, Nm, 8) = bm(E) coslo s + Xm(E)] (3, 6b) 
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geschrieben werden, wobei @,,(&), Om(§), Ym(&) und 7,,(&) gegeben sind. Setzt man 


if 


txm (6) 
2a3 bm(&) € 


Armlé) = aq enlé) f°", AanlE) = 


und kennzeichnet den Ubergang zur konjugiert komplexen Grosse durch einen 
Querstrich, so erhalt man gemass den Formeln (3, 2b) und (3, 2c) in 


3 “6 
Ax(é, n; s) = ase % D A” (é) e7t2amn +a 3 e-ims P A” (E) eee 


—4+3+1 —3§t+3jrl 


bp 


(3,7) 


as 7 2, bs ee ae 
Ayéé, n, s) = “se py: AW (é) e7t2aNn ae axere* py: AWE) e?2? ” 


—$+ 3+1 —$+3t+1 


zwei Fourier-Polynome, welche die nur fiir die diskreten 7,, definierten Aus- 
schlage der Schaufeln interpolieren und als Superpositionen von elementaren 
Schwingungsformen iiber dem Umfang des Schaufelkranzes darstellen. 

Da wir uns auf eine lineare Theorie in den Schwingungsamplituden be- 
schranken, in der sich die Krafte entsprechend den Ausschlagen tiberlagern, 
so reicht es aus, die aerodynamischen Reaktionen zu berechnen, die zu den 
elementaren Schwingungsformen 


AxE, n, 8) = ag Fo AME), Ay(E, 9, 8) = a3 e8*-™™7 AME) (3,8) 


gehoren. Dabei soll w eine beliebige reelle und 1 eine beliebige ganze Zahl sein. 
Die Grosse |”| gibt die Anzahl der Harmonischen auf dem Umfang an und 
beschreibt damit die Schwingungsform tiber dem Umfang. Von den beiden 
Funktionen A” (€), A$? (€), die die Schwingungsform der Profile iiber die Tiefe 
beschreiben, wollen wir nur voraussetzen, dass sie in dem Bereich 0 < € < 1 
stetig differenzierbar sind. 

In den technisch interessierenden Fallen haben diese beiden Funktionen 
einen sehr einfachen Aufbau. Ist das Profil starr, so kann es nur eine kombi- 
nierte Biege-Dreh-Schwingung ausfiihren. Diese kann als Uberlagerung der 
folgenden drei Grundformen dargestellt werden: 

1, Schlagschwingung normal zur Schaufel, Schlagamplitude A\” 


(n) 
AME =0, AME) = = = const; (3,9) 
2. Schwingung tangential zur Schaufel, Amplitude A!” 
n Ay") n 
AM(é) = mae AWE) = 0; (3, 10) 


3. Drehschwingung um die Vorderkante, Winkelamplitude A“) 


c A) 


AME) =0, APE) = = 


é. (3,11) 
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Im allgemeinsten Falle eines starren Profils ist also 


( 
A” () = Ay” ; A” (é) = Ay”) c Ay’ 


eo oan) vs Oe 


é. (3, 12) 


Entsprechend iiberlagern sich die Krafte. Es reicht also aus, die zu den oben 
angefiihrten drei Grundformen gehérenden Krafte zu: berechnen. Dement- 
sprechend kénnte man sich auf die Betrachtung linearer Funktionen von é 
fur A, und A, beschranken. Da dies aber bei der Herleitung keinen wesentlichen 
Vorteil bringt, so soll in der allgemeinen Theorie davon kein Gebrauch ge- 
macht, sondern ganz allgemein eine Methode angegeben werden, die es ge- 
stattet, die zu den Ausschlagen (3, 8) gehérigen Krafte zu berechnen. Die 
explizite Bestimmung der Krafte fiihren wir dagegen nur fiir die soeben ge- 
nannten drei Grundschwingungsformen durch. 


4. Gang der Untersuchung 


Wird ein ruhendes Streckengitter mit der Geschwindigkeit W, angeblasen, 
so legt sich, wie wir spater sehen werden, unter der Voraussetzung kleiner 
Teilung in die Vorderkanten der einzelnen Schaufeln ein Wirbel der Zirkulation 
I, der die Umlenkung besorgt. Zwischen den Schaufeln und dahinter herrscht 
die Abstrémgeschwindigkeit W,. 

Schwingt das Gitter, so treten zu diesem stationaéren Wirbel noch instatio- 
nare hinzu, die nicht nur in der Vorderkante liegen, sondern sich tiber das ganze 
Profil verteilen. Unter einem instationaéren Wirbel verstehen wir dabei das 
Feld eines Wirbels mit zeitlich veranderlicher Zirkulation einschliesslich des 
Feldes seiner Wirbelschleppe. Wird namlich ein Wirbel veranderlicher Zir- 
kulation angeblasen, so schwimmen nach dem Kelvinschen Satz in dem Masse, 
wie sich die Zirkulation andert, freie Wirbel ab. Diese unterliegen dem Anstrém- 
feld und dem von ihnen selbst und dem Ausgangswirbel induzierten Feld. Da 
wir uns auf eine lineare Theorie in den Schwingungsamplituden beschranken, 
so diirfen wir annehmen, dass die freien Wirbel nur dem Anstrémfeld unter- 
worfen sind, in unserem Falle also mit der Geschwindigkeit W, abschwimmen. 
Anders ist es bei einem Gitter mit gewdlbten Schaufeln; dort variiert die Ge- 
schwindigkeit langs der Schaufeln und ist erst von der Hinterkante ab kon- 
stant. Entsprechend schwimmen die freien Wirbel ab. Um das Problem nicht 
auch dadurch noch zu komplizieren, beschranken wir uns auf schwachgewdlbte 
Schaufeln, das heisst Streckengitter. 

Definiert man das Feld eines instationaren Wirbels wie oben, so ist die 
Zirkulation um den Ausgangswirbel identisch mit derjenigen, die man tiblicher- 
weise als gebundene oder tragende Zirkulation bezeichnet, und die Kraft, die 
auf diesen Wirbel ausgeiibt wird, ist wie im stationaren Fall gleich 9a VI’ 
und steht auf der Richtung von V senkrecht. Dabei ist V der Mittelwert der 
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Geschwindigkeit itber eine volle Umgebung des betrachteten Wirbelpunktes 
und J’ die gebundene Zirkulation. 

Wir bezeichnen mit W, ¢(&, 1m, s) fiir 0 << € <1 die Dichte der gebunde- 
nen Wirbel auf der m-ten Schaufel. Ausserdem miissen wir analog zum statio- 
naren Fall in die Vorderkanten isolierte instationére Wirbel legen. Deren Zir- 
kulation werde mit W,c G(j», 8) bezeichnet. In dem Punkte &, 77,, mit € > 0 
liegen dann ausser dem gebundenen Wirbel W, g(&, 7,,,5) noch diejenigen 
freien Wirbel, die sich von den davorliegenden gebundenen Wirbeln abgelést 
haben. Ihr Einfluss braucht aber nicht mehr gesondert betrachtet zu werden, 
da er nach der obigen Definition in dem Feld des zugehérigen gebundenen 
Wirbels enthalten ist. 

Die Wirbel liegen eigentlich auf den schwingenden Schaufeln, und ihre 
Starke ist so zu bestimmen, dass das von ihnen induzierte Feld die Normal- 
komponenten kompensiert, die sich durch die Anstrémung und das Schwingen 
ergeben. Die Anordnung der Wirbel auf den schwingenden Schaufeln braucht 
aber nur bei den stationaren Wirbeln mit der Zirkulation J’) zu erfolgen. Die 
instationdren Wirbel kénnen wir uns auf den stehenden Schaufeln angeordnet 
denken, da der hierdurch begangene Fehler in den Schwingungsamplituden 
quadratisch klein ist. 

Fiihrt man die hier skizzierten Rechnungen durch, so ergibt sich bei Be- 
schrankung auf eine lineare Theorie in der Teilung a/c fiir die beiden gesuchten 
Groéssen G(7) und g(&, 7) eine einzige Integralgleichung zweier Veranderlicher. 
Diese kann nach Einfithrung von Fourier-Entwicklungen tiber dem Umfang 
des Schaufelkranzes gelést werden. Abschatzungen dariiber, bis zu welchen 
Werten von a/c eine derartige Theorie sinnvoll ist, lassen sich im allgemeinen 
Fall nur schwer durchfiihren. Im stationaren Fall, der ja auch exakt behandelt 
werden kann, zeigt sich, dass fiir a/c < 0,8 die linearisierte Theorie den Ge- 
samtauftrieb mit hinreichender Genauigkeit liefert [5). 

Hat man die Gréssen Jy, G(7), g(€, 7) bestimmt, so kann damit die Kraft- 
verteilung auf den einzelnen Schaufeln nach dem Kutta-Joukowskyschen Satz 
berechnet werden, falls man die an dem betrachteten Punkt herrschende mitt- 
lere Geschwindigkeit kennt. Im Inneren des Gitters kann hierfiir W, genommen 
werden, da die Abweichungen hiervon in den Schwingungsamplituden linear 
sind und dies auch fiir die Dichten der dort liegenden Wirbel gilt. An der Vor- 
derkante ist im stationaren Fall die mittlere Geschwindigkeit gleich W,, (Mit- 
telwert von W, und W,), durch die instationiren Belegungen treten Abwei- 
chungen davon auf, die in den Schwingungsamplituden linear sind. Diese 
diirfen nicht vernachlassigt werden, da dort ausser der instationaren Zirkulation 
W,¢G(7) auch noch der stationére Wirbel J’, liegt, der nicht klein zu sein 
braucht, da wir uns nicht nur auf den Fall der schwachen Belastung beschrianken. 

Ganz entsprechend erfolgt die Rechnung, wenn ein Laufgitter durch ein 
inhomogenes Feld hindurchfahrt. 
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5. Geschwindigkeitsfeld einer schwingenden Wirbelreihe 
zeitlich konstanter Zirkulation 


Legt man durch die mit Wirbeln belegten Schaufeln einen Schnitt parallel 
zur Gitterachse, so erhalt man eine Wirbelreihe (siehe Figur 2). Die Wirbel 
dieser Reihe haben aber bei dem hier zu behandelnden Problem nicht alle die 
gleiche Zirkulation, sondern diese variiert mit der Blattnummer m und wieder- 
holt sich mit der Blattzahl 3 periodisch. Dasselbe gilt fiir den Abstand der 


pete 


Nien 
Tigi" 
Tp,)o 
Thy.) ox 
io 
Fig. 2 
Wirbelreihe aus einem Schaufelkranz. 


einzelnen Wirbel. Aus der Periodizitat ergibt sich, dass man das Feld der Aus- 
gangsreihe durch Uberlagerung von 3 Wirbelreihen mit der Gitterkonstanten 
a3 aufbauen kann. 

In diesem Abschnitt wollen wir nur den Fall betrachten, dass die Zirkula- 
tion der Wirbel zeitlich konstant ist. Dagegen lassen wir zu, dass sich die Lage 
der Wirbel mit der Zeit andert. 

In den Punkten P= miae-'? (m=0, +1, +2, ...) mégen Wirbel der 
Zirkulation I'(y,) mit I'(jm-+ 1) = I'(jm) liegen. Gesucht ist das von diesen 
Wirbeln erzeugte Feld. 

Die Wirbel in den Punkten P= (m+ k3)iae~'® (k=0, +1, +2, ...), die 
alle die gleiche Zirkulation ['(7,,) und den Abstand a3 haben, erzeugen das 
Feld 


sa IEG) a a ) } 
mz) Pee e ctg faa FE rl) eee ( 


das weit vor dem Gitter verschwindet. Das gesuchte Feld ergibt sich durch 
Uberlagerung der Felder fiir m = 0, 1, ..., 3 — 1, so dass also 
zeib 


: eb 3-1 It . 
Ut Vie 203 ¥ Tom) { ote( 275 — Nm) +? 
= m=0 


ist. Geht man zu dem Fall kleiner Teilung a/c, das heisst grosser Blattzahl 
iiber, so kann man die Summe naherungsweise durch ein Integral ersetzen und 
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erhalt, falls der Aufpunkt nicht in einem Wirbel liegt, also € + 0 ist, 


eiBb | eB ; ; ! 
UF mn 5, Vi eae fre ctg(— —an')+t;¢ dn’. (5,1) 


a3 


Dies gilt, falls die Wirbel ruhen. Schwingen die Wirbel und ist 
Az(Hm, 8) = 43[Ay(Hm, 8) + t A(m; $)] 


die Auslenkung des m-ten Wirbels, wobei wir vorlaufig noch annehmen, dass 
A, und A, reell sind, so ergibt sich aus (5, 1) 


it 
: eB . a e'B 
v-iv=< fre)| otal 2 
0 


: [z — Az(n’, s)]—2 n' +4 dy 


und fiir kleine Schwingungsamplituden 


e'B vn, + 68 — AI(n’) Az(n’, s) 
ee ca {Ten +8 On (5,2) 
{ cto( ze"? 
x | ota : z—n1) | an’, | 
falls € + 0. Versteht man unter 
k,(§, 0) ==| * tel E+ 29] 
+e ctg[— 3 E+ 07] — 2sing |, 
5,3) 
Ay, n) = |e ctg[ A é +29] ) 
y\S> 4 | Ge 1) 
—e"ctg[— E+ n| + 2icosp | 
sowie 
1 Og 
hiel€, m) = G | de? * ote[ A € + 279] 
—te-*Fotg|— 7 E+ xm] — 2c0s2 6}, 
bul, 2) = {| — e%? ctg[ eg se 
tele 5 7]) "a ee cte| fy E+2 1] 
—e*Petg|— E+ an] + 2sin2 6 
mit 
Pepe els |B 0 js 2 
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so kann damit das Feld der schwingenden Wirbelreihe zeitlich konstanter Zir- 
kulation in der Form 


OI’) A,(n/, 
= Ef (Te) bles 9) + AOA a ty 


UE bln —10') | an, ae 
OI’) Ay(n’, 
yaa / (i Ay(E, 9 — 1) + EE) a, 9 = 71) 
OI’) Ae(n’ : 
aa ) Bene ) | ay 


geschrieben werden. 
Wir wollen noch die Fourier-Integrale der Kerne (5, 2) und (5, 3) angeben. 
Beachtet man, dass fiir alle ganzzahligen n 


ferrets ta] trai +R) 
0 


ist, falls der Realteil von A, das heisst RA + 0 ist und dem Ausdruck n/|1| 
fiir » = 0 der Wert 0 beigelegt wird. Setzt man weiter |f| < 2/2 voraus, so 
ast fur == 0 


Mie 6 
(Roepe leh 


und man erhalt fiir € + 0 und + 0 die Fourier-Integrale 


1 n g —nx er n 3 ni ©.) 
Ri”) (E) =| eP aeer) Fre oa Ae $ 
sowie 
WG Sra el Gees ap ae Pe Ape age ‘| ae 
we =— {eM (ae + aioe Mer at) om | 


ADE) = Page (1 Ei Tae. sin B; RM (é) = = (1 zig +H) cosP . (5,8) 


Diese Formeln gelten fiir € + 0. An der Stelle €= 0 ist das Geschwindigkeits- 
feld unstetig, da dort die Wirbelreihe liegt. 


278 HEINz SOHNGEN ZAMP 


6. Geschwindigkeitsfeld einer Wirbelreihe mit zeitlich veranderlicher 
Zirkulation 


In der (€, 7)-Ebene mégen in den Punkten E=0, Hy = m/3 (m= 0, +1, 
+2, ...) Wirbel der Zirkulation I'(7,,) e'®* liegen, wobei (yj) IN Hm die 
Periode 1 hat. Von diesen Wirbeln schwimmen in dem Masse, wie sich die 
Zirkulation andert, freie Wirbel mit der Geschwindigkeit W, ab. Bezeichnet 
man die Dichte dieser freien Wirbel im Punkte € > 0, 7,, mit e(€, 7m, S), So 
ist also 

e(E, Nm 8) = £0, Mm 8 — §) 


und 
 I'(tim) e'°§ 4+ € e(0, Mm, 8) =0. 


Aus beiden Gleichungen folgt 
e(é, Ym> s) ee t @ oN ~ I'm) e* 2s ‘) od (6,1) 
Nach (5,5) wird das Feld einer ruhenden Wirbelreihe, die aus Wirbeln der 


Zirkulation J"(7j) in den Punkten = 0, 7 = 0, +1/3, +2/3, ... aufgebaut 
ist, fiir den Fall kleiner Teilung durch 


1 1 

il , = , , 1 ‘ , , t 
w= = | Py) hal. 1) dy’, v= | T(r) ble. n— 7) ay 
0 


beschrieben. Aus diesem Elementarfeld kann das gesuchte Feld der instatio- 
naren Wirbelreihe aufgebaut werden, und man erhalt hierfiir 


1 co 
: Lee « ’ . , —iwé’ & , ’ , 
u= eit — I dn’ P(n') [halG. 9 — 1') i @ fe 1°* bal — 8, — 1) a8"), 
0 0 
1 [ae] 


iws 1 i ' , , . C ~ioak’ ’ ’ 1 
poe =/ dy! '(y’) lev, n—m!)— iw | eo byl — ',n — 0') dB], 
0 


0 

(6,2) 

wobei die beiden inneren Integrale den Einfluss der abschwimmenden freien 
Wirbel darstellen. 

Es sei ausdriicklich darauf hingewiesen, dass diese Formeln nur Mittelwerte 
sind. Dadurch, dass wir friiher den Summationsprozess durch eine Integration 
ersetzt haben, haben wir eine «Verschmierung» in der 7-Richtung vorgenom- 
men, durch die die Spriinge, welche die Tangentialgeschwindigkeiten an einer 
Wirbelschicht erleiden, herausfallen. Da wir aber zur Erfiillung der Randbe- 
dingung nur die Normalgeschwindigkeiten brauchen, so stért diese Mittelbil- 
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dung nicht. Die Tangentialgeschwindigkeiten bendtigen wir nur zur Berech- 
nung der an der Vorderkante angreifenden Kraft, und zwar in der Form eines 
Mittelwertes tiber eine volle Umgebung. 

Zur Abkiirzung fiihren wir noch in (6, 2) die Kerne 


K(E, 9) = RalE, 1) — i [ e¥°F hlG — £, 9) a8, 
0 


K (6,9) = byl, 9) — 6 f oF h(E — 8,9) a8 
ny 
ein, deren Fourier-Transformierte sich mit 
v 1 oa BF 1 
%n NuU= 10” Xn nNn#+1 ow 
fiir 1 + 0 in der Form 
7 t n € —NzH 
KG) =F { dtiww,) (+) er 
1 n é \ nH 
E ’ eos 
— 40 (L+ er) (CF xy + oH my) ent, 
- § . : t (6,4) 
n a He) Mele —nx& 
KE) = 7 — Ft ma) (Farce) 
e ss =, n & n we 
+ (l—itwx,) e Ab ae gels 
+10 (1 + ii, (ero ye See gees | | 
bzw. fiir n = 0 
: iwé i. g —iw 
K(é) : (1 # ) sin B ei; eae) = ; (1 a ie] ) cos f e€ é 
(6,5) 


schreiben lassen. Diese Funktionen sind wieder fiir € = 0 nicht definiert. 

Mit den so definierten Gréssen erhalt man dann nach dem Faltungssatz fir 
die Fourier-Transformierten des durch (6, 2) gegebenen Feldes der instationa- 
ren Wirbelreihe die Ausdriicke 


(mg s) a u” (é) eis em ets ee K™ () | 
a eaten) (6,6) 


: 4 J r 
yl(E, s) as v”) (E) eis —. gi@s — HES) ; | 


durch die das gesuchte Feld beschrieben wird. 
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7. Geschwindigkeitsfeld um den schwingenden Schaufelkranz, 
Randbedingung 


Dieses Feld setzt sich zusammen aus: 
1. dem Feld der Anstro6mung 
u,= W,cose, v,= Wysing; 
2. dem stationéren Wirbelfeld (ruhende Wirbel der Zirkulation J in den 
Vorderkanten) 


Why dh ; 

Mp = —> ROE) = — (14737 ) sin 8; 
ne yb 

pee is ROE) = — oe (1+ iE ye) 0086: 


3. dem instationéren Feld, das durch das Schwingen der Vorderkanten- 
wirbel der a ec I) entsteht. Hierfiir ist nach (5, 5) 


die 0A, (ts aS) , OA 0, fo Ss ’ , 
Dy 0A, (0, 0A, (0 / 3 
= 2 byl 9 — 9) — 2) bole — 9") | dn 


mit den ane E ese © 
ult ==. orn £9 [H0( &) A ” (0) aay. we ) As” (0)] giehe. 
Ne =—t2m%n — Fo THM (E) (3) A (0) i RW”) (€) Ay” (0) Fosse 
4. dem Feld der instationaren Vorderkantenwirbel der Zirkulation 
W, c G(n) geen | 
Die Fourier-Transformierten dieses Feldes ergeben sich nach (6, 6) zu 


c y : 
PY =W, 26% KME cle, vf = W, = 6M Kee 


uy . 
5. dem Feld der instationaren Schaufelbewegung mit der Dichte 
W, g(§, 9) e'°*. 


Die Fourier-Transformierten dieses Feldes sind 


1 
ue) = W, =f ge") KOE — 8) db’ fee, 
0 


1 
of) = WE f gery KING — £9 ae ef 
0 
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Damit kann das gesamte Geschwindigkeitsfeld um die Schaufeln in der Form 


U = Uy, + Ug + Uz + Uy + Us, V =v, + vg + 0g + 4 + 05 (7,1) 


geschrieben werden. Wir benutzen dabei die Zeichen U bzw. V, um anzudeu- 
ten, dass es sich nur um Mittelwerte handelt. Weiter sei hervorgehoben, dass 
die Felder 2, 3 und 4 an der Vorderkante (€ = 0) unstetig sind. 

Um die Randbedingung aufzustellen, denken wir uns das Gitter mit der 
Geschwindigkeit (—W, cosx, —W, sinw) bewegt und ausserdem schwingend. 
Die Randbedingung verlangt dann, dass die Wirbelbelegung des Gitters so be- 
schaffen ist, dass in jedem Schaufelpunkt die Normalkomponente der dort von 
den Wirbeln induzierten Geschwindigkeit gleich derjenigen Geschwindigkeit 
ist, mit der sich der betreffende Schaufelpunkt in Richtung seiner Schaufel- 
normalen bewegt. Bezeichnen wir also mit x den Vektor, der die jeweilige Lage 
des betrachteten Schaufelpunktes beschreibt, mit nm die Schaufelnormale in 
diesem Punkt und mit p die dort von den Wirbeln induzierte Geschwindig- 
keit, also 

D = {ty + Ug + Uy + Us, Vg + V3 + Vg + U5}, 


so lautet die Randbedingung 
= ON =D ON. (72) 


Hat der Schaufelpunkt die Gitterkoordinaten &, 7, so ist 
t= MES Viet) 
mit 
x(é, 4, t) = —Wytcosat+cé+aznsinB + a3 Ai(ésn;s),, 
y(é, 7, t) = —Wytsina +a3ncosp + a3 4,(€, 7, s). 
Weiter ist a, 
aA ey, 
1 \ ; oF? il ie 


Geht man damit in die Randbedingung (7, 2) ein und beriicksichtigt nur die 
in den Schwingungsamplituden linearen Glieder, so erhalt man die Gleichung 


; a 0A 
—W, sina + W, cosa, = - | 
(753) 
0A a 0A 
+W,+4. ae ? Uy Ai 02 1s Cae Ee | 


die fiir 0 < € < 1, alle 7 und s erfiillt sein muss. Dies ist die Bestimmungsglei- 
chung fiir die Wirbelbelegungen J’), W, c G(n), W, g(é, 7). Sie zerfallt in einen 
stationaren und einen instationadren Anteil, die fiir sich befriedigt sein miissen. 
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Der stationare Anteil von (7, 3) liefert 


; a W, sin« 
—W,sinna=v,, also [y= con ph 


Die stationare Abstromgeschwindigkeit ergibt sich damit fiir § > 0 zu 


W, cos(« + B) 


rem 
W, = Wy cosa + U2 = W, cosa — —~ sinB = aig 


ri, a 
und fiir den Mittelwert der stationaren Geschwindigkeit in der Umgebung der 
Vorderkante findet man 


Fares : ie ; 
Wao = W, cosa — ae sinB, W,sin« — i cos | = W,, {cosa,,, SiN 4} 


mit 
tga 1 sin 


W. Ww, 


18% = 2—tgatgp’ on Pil et ae 


Auf jede Schaufel wirkt senkrecht zu w,, die Kraft 


K =09W,I,=5 0 We Cn eset ee 


Co aCLeCOS 2 


Fir den Auftriebsbeiwert ergibt sich also das bekannte Ergebnis 


FR hE Ben 6S (7,5) 


‘c cos B 
Wir fihren hier noch den modifizierten Auftriebsbeiwert C, ein, der durch 


a sin « 
Com 2 oe (7,6) 


definiert wird. Damit ist dann 
1 
y= zWicC,. (7,7) 


Der instationare Anteil der Randbedingung (7, 3) liefert damit 


0A, (E : ; 
Us + pies W, a 3 [- a n) +a A,(é, 7) | gies _ Vs 


c 


fiir 0 < € < 1. Geht man darin zu den Fourier-Transformierten iiber und ver- 
steht fir 0 << € < 1 unter H™(€) den Ausdruck 


(n)(£ 
HO) = 2. 8 (ae +i Av\e)) | 


(7,8) 
+ iam Cy[HME) AMO) — aN) AMON], | 
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der durch die Schwingungsform und die stationare Anstrémung des Gitters 
gegeben ist, so liefert der instationare Anteil die Integralgleichung 


1 
GO KOE) + i, e(é") KM(E — &) dé’ = HE) fir O<XE<1 (7,9) 
0 


zur Bestimmung der noch fehlenden Wirbelbelegungen G” und g)(é). 

Durch die Gleichung (7, 8) ist die Funktion H™)(£) nur fiir € > 0 definiert. 
Unter H™” (0) wollen wir noch diejenige Grésse verstehen, die man erhalt, wenn 
man mit € von positiven Werten her gegen Null geht. Es ist dann 


H (0) = = - “4 [AY’(0) + iw AM (O)] + S Ca[]n| 40) + in AM (O)] c2*”, 
(7,10) 
wobei 
nN 
‘|n] B 7? Bn 


gesetzt ist. 


8. Lésung der Integralgleichung des schwingenden Gitters 


Die Integralgleichung (7,9) soll fiir gegebenes H™(€) gelést werden. Wir 
beschranken uns dabei zunachst einmal auf den Fall ~ > 0. Mit den Kon- 
stanten 

An = (L+t@ xp) e?; b,= 1 (6°? vy —e-*? Kp) 


nimmt die Integralgleichung die Form 


(S 
GM ane" + bye FP] — a, f ge) oO ae 
0 


; é 
ae a, f e(E) enHl(S—5') dé’ a8 b, fae oe t@(E-§') dé’ a 2H (é) 
é 0 


an. Multipliziert man diese Gleichung mit e~”**, differenziert und multipliziert 
mit e”**, so erhalt man 


Gln a(x fe x) ens a b, (i ie n #) e408] a ge (é) [— Phe Bn fe bn | 

& & 
+2 ay(% + 7”) | @€) e "*S—*) ge’ _ bi wm + 0H) f ee") eto) ges 
i 


vA 
0 


= 2[H" (6) — ax HE). 
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Durch entsprechende Operationen kann man auch noch die beiden restlichen 
Integrale wegschaffen und gelangt dann zu der Differentialgleichung 


eg” (&) cosB = —H (&) — i[w + 2n QsinB] H”” (é) 


+[(n 2)? + 2m@n Qsinp] HX" (&) +7 (n Q)? AME) . 


Setzt man : 


HOME) = -/ HOE) ds’ und HYG) =— [ HE) a’, — (8,2) 
so folgt 
g"(6) cos = —H™"() — io + 2m Qsing] HE) | 


+[(n Q)?+ 2@n Qsinp] HY) (8, 3) 
tio (nQ)? HME) + £,6+F,. | 


Um schliesslich noch die Konstanten G”), E,,, K, zu erhalten, gehe man mit 
(8, 3) in die Gleichung (8,1) ein. Damit die sich dann ergebende Beziehung 
befriedigt ist, muss 


io (n Q)? 67"? H™ (1) — w? (2 Q)* e*? HYY(0) + w? (n Q)8 HY"(0) 


En nQ+1%m(nQ + 2 cos B) 
Fe qo? (n Q) e?*8 (14 n Qe-*8) (HM (0 0) — nQe-*® HW 0)] 
ae n 2+ 1 (nQ+ 2 cos f) 
nQe-**B(io+nQe®) H rag 
nQ +i (nQ + 2 cos) 
sowie 


G™ cosB = —H™)(0) 


een e EAS HVC) + [n Qe? + ionQe? + io] HM(0) + iwanQ ef HM\0)} 
NQ+io(nQ+ 2 cos f) 


sein. Damit ist die Integralgleichung (7, 9) fiir positive n vollstandig geldést. 
Ist 2 < 0, so nimmt die Integralgleichung die Form 


[-a, en%é ot PB Psa G™ = ay fe) en%(E—s') dé’ 
a (8,4) 
eet fener (&’) e —nx( *) dé et b feree et ol (€-&) de" = = 2H" (é) 


an. Nun ist aber, wie man sich leicht uberzeugt, fiir n + 0 


a_|,)(B) = 4,,(—B), a_1,(B) = a »(—B) sowie b_\,\(B) ie bin(—B) , 
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das heisst, die Integralgleichung (8, 4) erhalt man dadurch, dass man in (8,1), 
abgesehen von G™, g) und H”), die Grosse durch || und £ durch — f er- 
setzt. Da wir aber bei der Lésung dieser Integralgleichung keine Voraussetzung 
tuber das Vorzeichen von 8 gemacht haben, so gilt diese Symmetrieeigenschaft 
auch fiir die Lésung, und man erhialt allgemein: 


Ist n + 0, so wird die Integralgleichung des schwingenden Schaufelkranzes 
(7, 9) durch 


|n|Qe-?tBn 
|n| 2+ 14 (|n| 2+ 2 cos) 


G™ cos 8 = —H™)(0) + ( 


{-H™(1) + [| Q etn + iw [n|Q etn + ¢ w] HY (0) + iw |n|QeiPn He(0)} 


BS |n|2Q+ 4a (|n|Q+ 2 cos f) 
(8,5) 
sowie 
gE) cosB = —H™” (€) —i[w + 2|n| Qsing,] H”(é) 
+[(n Q)2+ 2@ |n| QsinB,] H%(é) 
+i (n Q)? Hy (€)+£,6+F, 
mit 
iw (n Q)? e~ 18m HM (1) — w? (n Q)?2 e1'n HM (0) + w? (|x| Q)3 HY (0) 
i= [n|Q + 7 w (|n|Q + 2 cos) g 
ow? |n|.Q e2*8n (1 + |n|Q e~*6n) [HM (0) — |n| Q e~*8n HY (0)] 
F, =—( |2|2 + 74 (|n|2Q+ 2 cos B) 
|m|Q e~2*8n (6 wo + |n|Q e4?m) H™ (1) ) 
a In| Q + ia (|n[Q+ 2cosp) 
gelost. 


Zu erledigen bleibt schliesslich nur noch der Fall n = 0. Die Integralglei- 
chung nimmt dann die Form 


is 
2G cos fie" = [ cosB e~*@@-*) dé’ = HE) 


0 
an. Daraus folgt 
G cos B = —H")(0) 
sowie 


2 (E) cosB = —H'(E) — iw HME), 


das heisst, die oben angegebene Lésung (8, 5) ist auch fiir » = 0 gultig. 
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Damit ist die Integralgleichung des schwingenden Schaufelkranzes kleiner 
Teilung vollstandig gelést, falls die Auslenkungen der Schaufeln gegeben sind. 


9. Ca-Anteil der gebundenen Wirbel 


Die gebundenen instationaren Wirbel sind nach Abschnitt 8 durch die Funk- 
tion H™)(é) vollstindig bestimmt. Diese Funktion ist durch die Schwingungs- 
form und die Anstrémung des Schaufelkranzes gegeben; sie zerfallt in den von 
C, unabhangigen Anteil 

(n) 
Hye = 4. 23 [248O sw ame] 9,1) 


is 0g 


und in den C,-Anteil 
HE) = Fn CLAVE) AO) — BNE) APO), O<E<1) | 


a 


EG (n) ; (n) 2iBn—|n| QetBne (9,2) 
= 2 Cq[]n| AMO) + im AQ (O)] en | | 
Entsprechend der Zerlegung von H schreiben wir auch 
gE) = go) + ar(é) und G” = Gy) + Gm. (9,3) 


Die vom Anstellwinkel « abhangigen Anteile kénnen, da in diese die Schwin- 
gungsform der Schaufeln nur durch Konstante eingeht, nach dem vorhergehen- 
den Abschnitt berechnet werden, und man erhalt 

e'Pn + In| Q—|n|QE 


Ga (E) = 1 Ca[|n| AY (0) + ¢ m AL (O)] w? ef Fn ‘[n|Q+ ta (|n|Q+ 2c0sB)’ 


pg. (n) in AM(0)) effn 1212 (L+ to) +i@ethn — 
Go = —%Ca[|m] AMO) + im AM O)Je [n|Q+ ie ([n[Q+ 2cosf) * 


(9,4) 


10. Krafte am schwingenden Gitter 


Um aus den Wirbelverteilungen die Kriafte berechnen zu kénnen, brauchen 
wir noch die in dem betrachteten Punkt herrschende mittlere Geschwindigkeit. 
Wie schon frither dargelegt wurde, kann hierfiir im Innern des Gitters die 
Geschwindigkeit W, genommen werden. Bei der Bestimmung der auf die Vor- 
derkante wirkenden Kraft reicht aber die Kenntnis der mittleren stationdren 
Geschwindigkeit wW.. = (veo, Yoo) in der Umgebung der Vorderkante nicht aus, 
man benotigt vielmehr die dort herrschenden mittleren Geschwindigkeiten ein- 
schliesslich der in den Schwingungsamplituden linearen Glieder. Der Mitte- 
lungsprozess ist dabei iiber eine volle Umgebung des betrachteten Punktes zu 
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erstrecken. Durch den Ubergang von der Summendarstellung zum Integral 
haben wir den Mittelungsprozess in der 7-Richtung bereits durchgefiihrt und 
brauchen ihn jetzt nur noch in der é-Richtung zu vollziehen. 

Unstetig an der Vorderkante sind nur die Geschwindigkeiten wp, V9} U3, Us 
und 4, v,. Da diese Funktionen fiir € = 0 nicht definiert sind, so verstehen 
wir unter (0, 7), v2(0, 7), ... die Mittelwerte. Weiter verstehen wir unter 
Uz(+0, 7), Vo(+0, 4), ... diejenigen Werte, die sich ergeben, wenn man von 
positiven €-Werten her gegen Null geht. Dann ist fiir » + 0: 


u(0) = — 52 |n| a W, © {AQ(0) sin2 B + AY(0) cos? BY, 
2 Oa = et |n| 0 W, — {AY (0) cos2 6 — Al" (0) sin2 p}, 


w(+0) = v0) — 1 W, 7 Cy n{Al(0) sin2 B + AY(0) cos2 B} 


sowie 
Pe aM LG n ee 
ul”(0) = WG 22 | ray (b+ io mn) 6? 
t+ yup = Fw en) OFF — a (08? x + e#7,) |. 
v”)(0) =W, 6% 21, pL +400 mq) ef? 


v(+0) = o(0) — =— WH, G™ cosf. 


Die Geschwindigkeiten u”)(0) und v¥")(0), die an der Vorderkante stetig sind, 
kénnen durch Integration berechnet werden. Einfacher ist es, die Randbedin- 
gung zu benutzen. Nach Abschnitt 7 ist 


vf)(+0) + (+0) = = H(0), 


wobei H)(0) durch (7,10) gegeben ist. Aus dieser Beziehung kann dann vf) (+0) 


berechnet werden, und wegen der Stetigkeit an der Vorderkante ist diese Grésse 
auch gleich v¥)(0). Beachtet man weiter, dass fiir 7 + 0 


w0)(0) = 4 7 8") 


ist, so kann auch w")(0) berechnet werden. 
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Ist 2 = 0, so ist 


uy (0) = 0; vg(0) =0, 

1 
uy” (0) = ‘ W,— G sing; 10) =—-zW, — G6 cosB , 
u(0) = 0; vy(0) = 0. 


Schliesslich gilt noch fiir die stationaren Geschwindigkeiten 


u,(0) = Wy cosa; v,(0) = Wosina , 
; if 1 : 
u(0) = eo sin B = == WitgBsina; v,(0) = ocr cosB = —= Wysina, 


so dass also die mittlere stationare Anstrémgeschwindigkeit der Vorderkante 
durch 


Ung = (0) + U9(0) = W, (1 +7—C,sin 8), 
(10, 1) 
Voo = U;(0) + v,(0) = W, —— “s C, cosB 
gegeben wird. 
Beriicksichtigt man schliesslich noch, dass in der effektiven Anblasgeschwin- 


digkeit der Vorderkante auch die Eigenbewegung der Profile - mit dem um- 
gekehrten Vorzeichen — aufgenommen werden muss, so erhalt man 


U.(0, 1, 8) = Mey + Us(0, 7, 8) + U4(0, 7, S) + 45(0, 9, 5) 
—~W,iw “4 AX, n) e 


co 
Mee M8 De ul(0) e~127"" — w+ u,(0, 7, s); 
—0o 


| (10,2) 
Ve(0, 1, 8) = Uso + U9(0, 4, 8) + 24 (0, 4, 8) + 25(0, 9, “ 
Wit 24 A,(0, n) ¢ ec | 
= ty, + 28 Se (0) e227 = 0, + 04(0, 0,8) | 
mit ae 
ui(0) = —W, iw =3 Amo) 
+14 on W,— = {> G™ (cosB + e'?") + (0) (10,3) 


1 
Semlee oF < x |n|[sin2 B AY” (0) + cos2 B AS (0)], 
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i } g A n 1 n 
#0) = Tso 2. APO) «WE [HO + 4 eos] | 
(10,3) 


1 
WC, o x |n|[cos2 B AY” (0) — sin2 6 AS (0)] | 
sowie 
W)(0) = —W, io +3. AMO) — > W, 2 6 sing, 


Durch diese Formeln ist die mittlere effektive Anblasgeschwindigkeit der Vor- 
derkante vollstandig beschrieben. In den letzten vier Formeln gibt das erste 
Glied jeweils die negative Eigengeschwindigkeit der Vorderkante an, das zweite 
Glied riihrt von der instationaéren Wirbelbelegung her, das dritte Glied schliess- 
lich, falls es auftritt, wird von den schwingenden Vorderkantenwirbeln der 
Zirkulation J) erzeugt. Diese Gleichungen gelten auch dann noch, falls das 
Gitter nicht schwingt, sondern instationar angestro6mt wird. Nur ist dann 
natiirlich fiir H die dem Problem entsprechende Grdésse einzusetzen, weiter 
sind A,, A, gleich 0 zu setzen, und das erste Glied ist schliesslich durch die 
gegebene Storgeschwindigkeit am Orte der Vorderkante zu ersetzen. 

Damit sind wir nun in der Lage, die Krafte anzugeben, die auf die Schaufeln 
des Gitters wirken, und zwar greift an einem Element der Lange c dé und der 
Breite 1 aus dem Innern der m-ten Schaufel nach dem Kutta- Joukowskyschen 
Satz, da dort die Dichte der tragenden Wirbel gleich W g(&, 7, s) ist, die 
Kraft 

APP = OC WEEE, Himvs) as, dP e=0 


und an der Vorderkante entsprechend die Kraft 
Poy = @ U,(9, Mm, S) [Lo + Wy ¢ Gm, s)] , 
Pow ar wa’ V,(0, Nm s) [Vo + W, @ Ces s)] 


an. In diesen ist der stationare Anteil noch mit enthalten; um ihn von dem 
instationdren Anteil zu trennen, schreiben wir 


U,(0, 4, S) = Ug + u(0, 7, s): V0, 7, s) = 9, 0210, 7,3) 


und erhalten in einer linearen Theorie fiir die an der Vorderkante angreifende 


Kraft 
Pry = 0 Uo Lo + 0 Yoo Wr & GQ, S) + @ He(0, Im, 8) Lo» 


Ye = =0 Vgq 1 = 0 Yoo Wi € Gm, 8) — @ V2(0, Pe NR 
in der jeweils die ersten Glieder die stationare Kraft angeben. Ftir die gesamte 
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an der m-ten Schaufel angreifende Luftkraft ergibt sich somit 


P, = 0 U5 P9 + 0 Ue(0, Nm 8) 9 + © Yoo Wi € G(%m, S$) 
1 
+geW? [ eG, 7m.s) a, } (10,4) 
0 


Py = —0 Yoo Ly — 0 Ve(0, Ym; 8) P19 — @ Yoo Wy ¢ Gm, S) 


sowie fiir das an der Vorderkante angreifende und kopflastig positiv gezahlte 
Moment 


1 
M = oc? W, | € 66, Mm) dB. (10,5) 
0 


Durch die Formeln (10, 4) und (10, 5) ist die eingangs gestellte Aufgabe voll- 
standig gelést. Diese Formeln gelten fiir eine lineare Theorie in der Teilung a/c. 
Da die stationare Theorie — wie schon oben erwahnt — fiir Teilungswerte bis 
zu 0,8 eine sehr gute Naherung fiir den Auftrieb liefert, so wird man hoffen 
diirfen, dass fiir Schwingungsformen des Kranzes, bei denen die Anzahl | 7| 
der Harmonischen auf dem Umfang klein gegeniiber der Blattzahl 3 ist, die 
vorliegende Theorie brauchbare Werte ergibt. 

Bei der Herleitung der Formeln wurde vorausgesetzt, dass die Auslenkungen 
der Schaufeln dem Gesetz (3, 8) unterworfen sind, also harmonisch iiber dem 
Umfang verteilt sind. Hat man es allgemeiner mit einer Schwingung der Form 
(3, 6) bzw. (3, 7) zu tun, so sind die Krafte entsprechend zu iiberlagern. 

Fir die in die Gleichungen (10, 4) und (10, 5) eingehenden Grdssen gilt 
noch der folgende Wegweiser: Die stationaéren Gréssen t, Vo, W, und Ty 
werden bei gegebener Anstrémgeschwindigkeit W, und gegebenem stationarem 
Anstellwinkel « durch (7, 4), (7,6) und (7, 7) bzw. (10,1) geliefert. Die insta- 
tionéren Wirbelbelegungen G bzw. g sind den Formeln (8,5) und (8, 2) zu 
entnehmen, wobei die Hilfsgrésse H\(€) durch (7,8) mit (5,7) gegeben’ ist. 
Die mittleren instationaren Geschwindigkeiten an der Vorderkante ,(0) bzw. 
v(0) werden durch die Gleichungen (10, 3) mit (7, 10) gegeben. 


11. Zusammenstellung der Luftkrffte fiir die 
wichtigsten Schwingungsformen 


Es sollen hier explizit die Luftkrafte angegeben werden, die zu den drei ele- 
mentaren Schwingungsformen iiber die Tiefe, namlich 
1, Schlagschwingung normal zur Schaufel (3, 9), 
2. Schwingung tangential zur Schaufel (3,10), 
3. Drehschwingung um die Vorderkante (3,11), 
aus denen sich die allgemeinste Biege-Dreh-Schwingung eines starren Profils 
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linear kombinieren lasst, gehéren. Vorausgesetzt wird dabei, dass die Schwin- 
gungsform uber dem Umfang harmonisch ist, also in der Form (3, 8) dargestellt 
werden kann, wobei eine beliebige positive oder negative reelle Zahl und n eine 
beliebige Zahl ist. 

Durchweg verstehen wir dabei unter N die Grésse 


N = (n Q)? + w? (|m| 2 + 2 cos B)? 
und setzen noch 


Oey und 9572) 2 —or. 


Pett 
ala 
a) Schlagschwingung normal zur Schaufel 


Auslenkung der m-ten Schaufel: 
Al) Ay = ae ei es—t2aA Nm 


Krafte auf die m-te Schaufel: 


a ss 72 
ia 0 Ww? oie Cae 12mm nm 


; Cc / 2 3 
x A, +ioB, + G Ca) Ces Os) a: (e Ca) oy ae By). 
P al 0 Ww? A, a piws—i2annm 
Me c 


x { (= ¢.) (4, +4 b,) + (= OP ee ee b,,) 


ee 12 o% sin2 
b = 12 NcosB {n B 


+ w2 [2064+ 12 o? cosB + 24 o% cos? B + 24 o cos2 B cos A] 
— wi (o + 2 cos) [24 cosB + 24 ocos?B + 8 o? cosB + o°)}, 


1 : des 
ae IN {o3 sinB + w, «7 (o cosB — 2 sin? B) 


+ w% (o + 2 cos) [2 0cos2B + o* cos f}}, 


5 ees ae {o8 sin? B + w2 oa (o + 2 cos) [o sin? B — cosB cos2 B]} , 
at phe i 2 cos B) w 
ee a Nicos eh Coa 2 cos f) 8 ao sin cosB{o + 2cosf) o, 


+ w2[o4 + 8 o® cosB + 24 o* cos*B + 32 o cos? B + 16 cos? fj}, 


=A af = $ - ict 
ES + 2 cos) o? sin B 
co 2N | o, Sacoe tte B) 


+ w, 6 cos B[o? + 40 cosB + 4cos2 f])}, 


o*cosf, 


1 
Big MEN 
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a, = —_ {w, 6? sin B (o + 2 cos B) + w2 2 0cosB[o* + 4a cosB + 4 cos 2 B] 
+ w% 2osinB (o + 2 cosf) (6 + 4cosB)}, 


a, = se {—o3 sin B cos B + w, o? cos B — w3 a sin B cos B [o + 2 cos B}?}, 


by = ae o? [1+ cos2£ + ocosf] — o w, sinB [o? + 4 0 cosB + 8 cos? B] 


+ w2 (o + 2 cos) [4 o cos?B + o* cosB + 4 cosB}}, 


1 il 


Moment um die Vorderkante dey m-ten Schaufel (kopflastig positiv) : 
: c 4 iws—i2a 
M= 0 Wia A, | 6, + to D, + (= Ca) [C,, + iw Dy} fe’ ite 2: 
mit 
= RE Bh sma 3 4 
fea Nios B {@, 12 o? sin B — w2 [4 o* cosB + of] 
+ w3 osinB [6 o? + 32 ocosB + 48 cos? B] 
+ wf (o + 2 cos) [o3 + 8 o? cosB + 24 ao cos?B + 24 cosB]} , 


1 
Cy. = Gop (ene? sinB (6 cos B + o) + w§ o \3. cos2 B + ocosf) (o + 2 cosB)}, 
1 
a TS 2 4 3 y 3 aj 
D,, 12 N cosB {2 0*+ 60% cosB + ow, 2 0% sinB 


+ w? [ot + 7 o? cosB + 16 o? cos? B + 12 o cosB cos2 B] 


— w§ (o + 2cosf) osinB (4 cosB + a)}, 


1 
D, = TIN {@, 07 (3 cos2 B + o cos B) + w2 (3 sin2 B + osin§) («6 + 2 cos )}. 


b) Schwingung tangential zur Schaufel 
Auslenkung der m-ten Schaufel: 
Ax = A, e®-828m, Ay = 0. 
Kraft auf die m-te Schaufel: 
a c ; Cc 2 F ent 
Py = @ Wi A; 7 (= Ca) (4, +2 B,) + (= Ca) (4, +1 B,)h gOS t22N tes 


{ 
\ 
(< Ca) {a,, +14@ b,,\ et OS FanA tm 


1 
Palas aN {2 o® cos B + w?, o cos B [o? + 40 cosB + 4 cos2 B] 
—o(o+ 2cosf) (2sin2 6 + osinB) 8}, 
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23 = Tae {o* sin B cos B + w, 6? cosB + w2 sin B cos B (o + ZaCOS 3) 2 
2 ere i. ; ; 
Oe ee ea sinB + w, 6 [o sinB + 2osin2 8 + 4cosf sin2 B] 
+ w7, (6 + 2 cos 8B) [4 cos + 40 cos? B + a? cos fj}, 

ep Oy 

= 1 3 a 2 y 
LF oor eg cos? B + w2 o cos B (o + 2 cos #) [o cos B + cos2 B)}, 
re, 


Moment um die Vorderkante der m-ten Schaufel 
é : 
M = 0 Wiad, (= Ca) {C, +ioD,} 
mit 
1 
C= Amita {w2 o7 (3 cos2 B + ocosB) + w§ (3 sin2 B + asin) (o + 2cosf) o}, 


1 : 
= jaw {On o* (3 sin2 B+ osinB) — w? o (o + 2 cos) (3 cos28 + acosf)}. 


c) Drehschwingungen um die Vorderkante 


Auslenkung der m-ten Schaufel: 
AymOe Ayano Apel OS “ARN ime. 
Kraft auf die m-te Schaufel: 


Py = a) Wia Ag| Ag, + 1 o By + (= Ca) (Aq, a4 wo By) GiPS= 12 Ahn” 


P,= 0 Wiad, (<= Ca) fay, +io Cee te cee 


mit 
a7 aes o? (o + 2cosf)?+ 12 w, o sin (o* — 8 cos? B) 
+ w2 [904+ 80 o8 cosB + 192 o? + 48 ocosB (2 + 3 cos2 B) + 96 cos? B] 
+ w3 [2 o® sinB + 4 (0 + 2 cos) (3 o sin2 B + o? sin f)] 
+ w4 (o + 2 cos) [24 cosB + 24 a cos? B + 8 o* cosB + oF ]}, 
1 . 
ee abodes et o? sin2 B (« + 2cosf) — w, 12 a cosB (6 + 2 cos cos 2 B) 


— w2osin2 B [48 cos?B + 30 a cosB + 5 o*] 
+ w8 2 0cosB (o + 2 cos B) [3 cos2 8 + o cos B)}, 
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2 Pad o? sin B cos B — [4 o4 + 12 o® cos B — 48 o cos2 B cos B] 
» 24NcosB | o, 
+ ow, 40 sinB [o? + 12 6 cos? B + 24 cos? B] 
— w2 [6 o4 + 54 63 cosB + 192 o? cos*B + o cos f (192 + 120 cos2 B) 


+ 192 cos? f] — w§ o sin B (o + 2 cosB) (20+ 8 cos)}, 
1 


= — NPE cos2 Bf + o cos f] + 24 0sin2 B cos*B 
2s 24NcosB |p, NE: B Bi 
— w, 0 cos [10 6? cosB + o (12 + 30 cos2 B) + 48 cos B cos2 B] 
— w? o sin2 B (o + 2 cosB) [a + 6cosB}}, 
a4, = aq {6 8 (cos2.B + 0 cos f) — 12 ®, 6 Sin2 B cosB 


+ w2 0[24 cos B cos2 B + o (30 cos? B — 9) + 5 a? cos B] 
+ w§ osinB (o + 2 cos) [6cosB + o}}, 


= 1 : . 
by = Gan {or O78 B e+ 2 cos B) + 602+ 12 0cos2fcosB 


+ @, 0 sinf [48 cos? B + 30 o cosB + 5 a?) 

— w2 (o + 2cos) [o? cosB + 3 ocos2 B)}. 
Moment um die Vorderkante der m-ten Schaufel 

M = 9 Wiac Ag {Cy +i wD, beteseannm 
mit 
C, = eo te {30 
0 180 N cos B 

+ w,[21 of + 195 03 cos B + (30 + 240 cos? B) o? + 180 ocosB cos2 B] 

— wo sin B [30 o? + 180 o cos B + 240 cos? B] 

— wA (ao + 2 cos B) [4 63 + 33 o? cos B + (60 + 45 cos2 B) o + 120 cos f]}, 


=! { 90 ae Ae 
~ “180 N cos ma 6 
Di, 180 N cos B \ @, Oo sin B +156 lo +4o cos B an ] 


+ ow, 60 osin B [o* + 4a0cosB + 6 cos? B] 
+ wj7, [20 of + 180 o8 cos B + (45 + 600 cos? B) o® +360 o cos B (2 + cos2 B) 
+ 720 cos? B]}. 


o* + 90 o® cos B 


12. Biegeflattern eines Schaufelkranzes 


Um ein einfaches Beispiel fiir die Behandlung des Flatterproblems eines 
Schaufelkranzes zu geben, soll untersucht werden, ob die Schaufeln eines Kranzes 
auch bei Beriicksichtigung der Luftkrafte Biegeschwingungen ausfiihren kénnen. 
Dabei wollen wir annehmen, der Staffelungswinkel des Kranzes sei B = 0, die 
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Schaufeln im Vakuum haben die Eigenfrequenz », und seien iiber den Fuss hin- 
weg nicht gekoppelt. Weiter sei M die Masse dex Schaufel pro Spannweiten- 
einheit, und mit F,(¢) werde die auf die Spannweiteneinheit der k-ten Schaufel 
wirkende Normalkraft bezeichnet. Fiir die Bewegungsgleichungen erhalt man 
dann 


if 1 
AV EP)” AV aa Felt) o: (h = 1, 2,.+2,.3) (1254) 


dabei sind die F, selbst Funktionen der 4y,, Ays, .. , Ay. 
Es ist zu untersuchen, ob es eine Geschwindigkeit W,, eine Frequenz v sa 
42 nicht samtlich verschwindende Konstanten c,, so gibt, dass die 


i t2avt 
ANG, = Bren 


bei der Geschwindigkeit W, das Gleichungssystem (12,1) befriedigen. 
Um zunachst den Zusammenhang zwischen den Auslenkungen und den Kraf- 
ten herzustellen, berechnen wir nach (3, 2b) die Fourier-Transformierten 


Damit wird dann 
3 
a 1 lt) pi2nvt—-i2anne 
Ay ee e (A252) 
Zu dem Anteil cc” ¢1?7”!-#2"""k von Ay, gehért nach Abschnitt 11 die Kraft 
in der y-Richtung 
PY) = gac™ Wi efes—t2amme (FM) + 4 wo FM) (12, 3) 
mit 


a 2 

Fy) = Ay, + (5 Ce) ie ie Ca Aug 
re 2 

Fe) = B, + (= Ca) B, + (= Ca) Bae 


wobei die A und B reelle Funktionen der Gitterkonstanten 2, der reduzierten 
Frequenz w und von » sind. Gemass dem Superpositionsprinzip ist dann 


3 eee 
F(t) mega WH Sac) ets ene (Ea 2 kya) 
=g4g +1 


Geht man damit sowie mit der Zerlegung (12, 2) in die Bewegungsgleichungen 
ein, so erhalt man mit 


das Gleichungssystem 
» cl) e—t2mnmk fy — «2 — w (FY) + i FO) bar 0 6 (= eaten 8) 
—gtgtl 


+ 
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Wendet man darauf die Fourier-Transformation an, das heisst, multipliziert man 
die k-te Gleichung mit e'2™™"k und summiert iiber k = 1, 2,..., 3, so findet man, 
dass 


o) £2 — wo? — pp (FO) 4 4 oF™)$}=0 (m=—-3+3+4+1,...,3) (12,5) 


sein muss. Die Schaufeln des Kranzes kénnen also dann und nur dann Biege- 
schwingungen ausfiihren, wenn das Gleichungssystem (12, 5) eine nichttriviale 
Lésung zulasst. Dies ist aber dann und nur dann der Fall, wenn es wenigstens 
eine unter den ganzen Zahlen —3 +3 +1, ..., 3 so gibt, dass fiir diese Zahl n 


of — wo? — pw (FO +i FM) =0 
ist, das heisst aber 


FM =0 und o2%= wo? + pF 


ist. Die erste dieser beiden Gleichungen bestimmt die kritischen w-Werte, aus 
der zweiten kénnen dann die dazugehérigen kritischen Geschwindigkeiten be- 
rechnet werden. Dabei interessiert der Wert von w, der den gréssten Wert von 
@, liefert, denn @, ist ja der kritischen Geschwindigkeit umgekehrt proportional. 
Da nun yu im allgemeinen klein ist!), so lauft dies darauf hinaus, ein méglichst 
grosses w zu bestimmen, fiir das wenigstens eines der F\”) verschwindet. 

Im Fall der Biegeschwingungen ergibt sich mit 6 = 0 nach (12, 4) und Ab- 
schnitt 11 


c c 
(o + 2)4 iY goer rie o? 
= — ———— 1a, + 5 i} w2 + W, ont é 
4No, (o + 2) (2 + 0) (2+ 0)? 


Cc 
o—C, 
EO: eee ee o= |n| ee: 
n (o + 2). 5 = a 3 


Das Maximum wird angenommen, wenn o in der Nahe von 2 liegt. Da die Schau- 
felzahl 3 im allgemeinen gross ist, so diirfen wir voraussetzen, dass der Wert 2 
von @ angenommen werden kann und erhalten als kritische reduzierte Frequenz 


Le 
Dent ert Ca 


1) Es wird an ein Gas als str6mendes Medium gedacht. 
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und mithin (falls ~ = 0 gesetzt wird) als kritische Abstrémgeschwindigkeit 


16a wa 
GC 


oder mit dem iiblichen Auftriebsbeiwert c, fiir die kritische Anstrémgeschwin- 


digkeit 
Kae P(e) il c 2 
Wo, erit = : > /3 3 (5 her cal : 


Wie krit 
a 


a 


Damit ist also gezeigt, dass ein belasteter Schaufelkranz schon allein mit dem 
Freiheitsgrad Biegung flattern kann; wesentlich ist allerdings, dass die Schaufeln 
die Moglichkeit haben, mit Phasenverschiebungen zu schwingen, also nicht gegen- 
seitig abgestiitzt sind. Weiterhin darf auch die Eigenfrequenz der Schaufel nicht 
zu hoch liegen). Schliesslich ist noch zu beachten, dass in der vorliegenden 
Untersuchung ein inkompressibles Medium und anliegende StrOémung voraus- 
gesetzt wurde. Ist die Strémung abgerissen, so schwingt der Ablésepunkt; da- 
durch erfahren auch die Luftkrafte eine Phasenverschiebung, die im allgemeinen 
die Luftdampfung herabsetzt. Weiter ist noch zu vermuten, dass die Kompres- 
sibilitat der den Schaufelkranz umgebenden Gassdule von wesentlichem Einfluss 
ist. Auf diese Frage soll in einer spateren Arbeit eingegangen werden. 
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Ae, 


Summary 


For an investigation on the vibration behaviour of a ring of blades of an 
axial compressor, knowledge of the aerodynamic forces on the individual blades 
is required. Since, at small gap-chord-ratio and possibly at high aerodynamic 
loading of the blades, the aerodynamic interference between the blades is not 
negligible, a method including the interference is given for evaluating the aero- 
dynamic forces acting on an arbitrarily vibrating ring of blades. This method is 
valid for thin and slightly cambered blades, arranged with small gap-chord ratio, 
in nonviscous incompressible medium. 


(Eingegangen: 22. August 1952.) 


1) Zu einem ahnlichen Ergebnis gelangten auch Cu. BeLLenor und J. Lative p’Eprnay [6], 
wobei allerdings sehr stark vereinfachende Annahmen Uber die Luftkrafte an einem schwingenden 


Gitter gemacht wurden. 
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Wirbelsysteme in Schaufelgittern und Turbomaschinen 


Von WALTER TRAUPEL, Winterthur?) 


1. Einleitung 


Uber die Wirbelsysteme, die beim Durchtritt eines Strémungsmittels durch 
eine Schaufelreihe auftreten, besteht bis heute kein klares Bild. Betrachten wir 
etwa ein gerades Schaufelgitter, so fliesst diesem das Str6mungsmittel niemals ~ 
mit vollig gleichformiger Geschwindigkeitsverteilung zu, da sich an den Begren- 
zungswanden des Zustrémkanals Grenzschichten bilden. Die Endpartien der 
Schaufeln werden daher langsamer angestrémt als die mittleren, weshalb auch 
die Zirkulation um die Schaufel langs ihrer Erstreckung variert. Es wird daraus 
haufig gefolgert, dass, entsprechend dieser variablen Zirkulation, Wirbelschlep- 
pen von der Austrittskante stromabwarts gehen miissen — gemiss dem dritten 
Helmholtzschen Wirbelsatz — ahnlich wie beim einzelnen Tragfliigel im freien 
Luftraum. Anderseits entsteht bei jeder Ablenkung einer Str6mung der soge- 
nannte Sekundiarwirbel, da das wesentlich durch die Kernstr6mung bestimmte 
Druckfeld der Grenzschicht aufgepragt wird und dort entsprechende Querbe- 
wegungen herbeifiihrt. Nun stellt sich die Frage, ob es sich bei den beiden eben 
erwahnten Wirbelsystemen tatsachlich um zwei verschiedene Effekte handle 
oder ob nur zwei Aspekte oder Erklarungsversuche ein und desselben Effektes 
vorliegen oder ob wir gar zwei widersprechende Theorien vor uns haben, von 
denen mindestens eine falsch ware. Ziel der folgenden Ausfiihrungen ist es, diese 
und verwandte Fragen zu analysieren. — In den Figuren werden Wirbellinien 
entsprechend Stromlinien mit Pfeilen angegeben, wobei Pfeilsinn und Drehsinn 
des Wirbels durch die Rechtsschraubenregel einander zugeordnet sind. Figur 1 
stellt die Darstellungsweise klar. 


2. Der Einzelfliigel 


Beim Einzelfliigel im freien Luftraum gehen bekanntlich von den Fliigel- 
enden Wirbelschleppen stromabwarts, die den «nduzierten Widerstand» be- 
dingen. Herbeigefiihrt wird diese Erscheinung durch den Druckunterschied 
zwischen Fliigelober- und -unterseite, der eine Umstrémung der Fliigelenden 
im Sinne der Pfeile 7, Figur 2, bewirkt. Die Luft gleitet daher, wie die Pfeile 2 
zeigen, schief iiber die Fliigeloberseite und entsprechend im entgegengesetzten 
Sinne tiber die Unterseite, und von der Austrittskante gehen folglich Unstetig- 
keitsflachen — das heisst Wirbelflachen — ab, wie die Pfeile 3 darstellen. 


1) Gebriider Sulzer AG., Winterthur. 
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Man findet oft die Darstellung, wonach das Auftreten dieser Wirbelschlep- 
pen unmittelbar aus dem dritten Helmholtzschen Wirbelsatz folgt (Zirkulation 
um einen Wirbelfaden ist konstant, daher kann Wirbelfaden im Inneren der 
Fliissigkeit nicht beginnen oder enden, sondern ist in sich geschlossen oder er- 
streckt sich ins-Unendliche). Dies ist aber nicht korrekt, da der dritte Helm- 


Higa t Fig. 2 


Fig. 3 
Mathematisch denkbares Wirbelsystem am Einzelfligel. 


holtzsche Wirbelsatz ein rein mathematisches Theorem ist; er macht eine geo- 
metrische Aussage iiber jedes denkbare Vektorfeld, unabhangig von jeder phy- 
sikalischen Bedeutung. Die Entstehung der Wirbelschleppe ist aber durch die 
Gesetze der Dynamik bestimmt, kann also niemals durch mathematische Uber- 
legungen allein vorausgesagt werden. In der Tat sind mathematisch auch ganz 
andere Lésungen denkbar. So kénnte man etwa annehmen, der Fligel beein- 
flusse den Luftraum nur innerhalb des Bereiches seiner eigenen Spannweite, 
wahrend alles ausserhalb Liegende unbeeinflusst bleibt. Dann bilden sich, von 
den Fliigelenden ausgehend, senkrechte Trennflachen . vel. Figur 3 2 welche 
beeinflusstes und unbeeinflusstes Gebiet scheiden. Diese Losung erfiillt, wie 
noch unabsehbar viele andere, den dritten Helmholtzschen Wirbelsatz ebenso 
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wie die Lésung nach Figur 2. Es sind die Gesetze der Dynamik, die aus der 
unabsehbaren Fiille der an sich denkbaren Lésungen diejenige nach Figur 2 — 
auswahlen. 


3. Gerades Schaufelgitter, primares Wirbelsystem 


Figur 4 veranschaulicht das Wirbelsystem einer zwischen zwei parallelen 
Wanden verlaufenden Parallelstromung, deren Geschwindigkeitsprofil ebenfalls 
dargestellt ist. Es besteht aus im Bereiche der Grenzschichten quer zur Str6- 
mungsrichtung verlaufenden Wirbelfaden. 


N 
bate) 2? / 
[0) ® ra 
[o) 8 
oe S @ 
[o) 8 
on rs ® 
[o} @ 
°° 2 @ 
N SY 
Fig. 4 Fig. 5 
Geschwindigkeitsprofil und Wirbelsystem der Primares Wirbelsystem der Gitterstrémung. 


grenzschichtbehafteten Parallelstrémung. 


Wir betrachten nun eine solche Strémung, die in ein gerades Schaufelgitter 
eingeleitet und dort umgelenkt wird, um wieder als Parallelstromung fortzu- 
strémen (Figur 5). Die im Bild gezeigten Wirbellinien reprasentieren die Wirbel- 
faden der Grenzschicht der unteren Begrenzungsebene. Auf der Zustrémseite 
kommen sie aus dem Unendlichen und steigen dort, wo sie auf ein Schaufel- 
profil auftreffen, an diesem hoch — entsprechend der Grenzschicht am Schaufel- 
profil selbst — und gehen schliesslich der oberen Begrenzungsflache entlang ins 
Unendliche zuriick. Genau so verlaufen die Wirbelfaden auf der Austrittsseite. 
Innerhalb der Schaufelkanale sind die Wirbelfaden geschlossene Ringe, die 
an der Saugseite des Profils auf- und an der Druckseite absteigen. 

Bei endlicher Grenzschichtdicke erfiillen die Wirbelfaden den Raum stetig 
und schrumpfen in Unstetigkeitsflachen zusammen, wenn die Grenzschicht- 
dicke gegen Null strebt. In diesem Grenzfall springt also die Geschwindigkeit 
unmittelbar an den Wandflachen auf Null, was auch bei reibungsfreier Stré- 
mung so gedacht werden muss, soll die Kuttasche Abflussbedingung erhalten 
und das d’Alembertsche Paradoxon vermieden werden. — Bei endlicher Grenz- 
schichtdicke fallt die Zirkulation um das Profil gegen die Schaufelenden hin 


Vol. IV, 1953 Wirbelsysteme in Schaufelgittern und Turbomaschinen 301 


stetig auf Null ab. Daraus folgt nach dem dritten Helmholtzschen Wirbelsatz 
noch keineswegs, wie oft falschlicherweise geschlossen wird, dass von der Schau- 
felaustrittskante Wirbelschleppen abgehen miissen. Dies ware eine physika- 
lische Aussage, und eine solche kann aus dem genannten Theorem niemals 
gewonnen werden. Wie es im vorliegenden Falle erfiillt wird, zeigt eben Figur 5, 
und es muss dazu vom mathematischen Standpunkt aus nichts Weiteres bei- 
gefiigt werden. Mit der Annahme der Existenz der grenzschichtbehafteten 
Strémung ist auch die Existenz des dargestellten Wirbelsystems denknotwendig 
gegeben, denn dieses ist ja nichts anderes als eine besondere Art der mathema- 
tischen Beschreibung des Vorausgesetzten. 

Wir nennen das hier beschriebene Wirbelsystem das primdirve, da es ja mit 
dem Vorhandensein der Strémung zugleich schon besteht, und zwar — in beson- 
derer Form — selbst im Grenzfall der Potentialstr6mung. 


4. Sekund4res Wirbelsystem 


Betrachten wir wieder die Str6mung durch das gerade Schaufelgitter mit 
Grenzschichten an den ebenen Begrenzungswanden. Innerhalb des gekriimmten 
Schaufelkanals sind die in der Grenzschicht strémenden Fliissigkeitsteilchen 
den im Kanal herrschenden Druckgradienten ausgesetzt, doch sind ihre Ge- 
schwindigkeiten nicht so gross, dass beim Fortschreiten auf ungestérten Bahnen 
dynamisches Gleichgewicht entstaénde. Daher werden die wandnahen Teilchen 
in Richtung des Druckgefilles seitlich abgedrangt, das heisst also nach der 
Saugseite hin. 


Fig. 6 
Sekundarwirbelsystem der Gitterstromung. 


Das entsprechende Wirbelsystem zeigt Figur 6, in welcher wieder das Wir- 
belsystem langs der unteren Begrenzungsebene dargestellt ist. Uber den Quer- 
schnitt eines einen Schaufelkanal erfiillenden Stromfadens ist die Starke des 
Sekundarwirbels in bestimmter Weise stetig verteilt, wie in Figur 65 und ¢ 
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veranschaulicht. Nach Austritt aus dem Schaufelkanal erhalten die gegen- 
seitigen Begrenzungsflaichen dieser Stromfaden den Charakter von Diskonti- 
nuititsflachen (vgl. ebenfalls Figur 65 und c). 

Berechnen wir nun die Zirkulation der Sekundarstrémung langs einer Kon- 
trollkontur der Art ABCD (Figur 7). Sie ist 


B D 
ar = [ wdx + o(B) dy + | wdx—»(A) dy | 
A Cc (1) 


B B 
= [ ude + o(B) dy— [ (w+ $* ay) dx — v(A) dy. 
A A 


Da sich das Kontrollgebiet gerade iiber eine Teilung erstreckt, ist v(B) = v(A). 


Ferner ist gemass der Kontinuitatsgleichung 0u/0y = — 0v/0x. Folglich wird 
Gleichung (1) 
z “i 
Ae dy [ oY dx = dy | dv =[v(B) — v(A)] dy =0. (2) 
A A 


Ist aber dJ’ = 0, so muss die Wirbelstarke des in der Kontrollkontur hegenden 
Teiles der Trennfliche entgegengesetzt gleich der Wirbelstarke der gesamten 
Sekundarstr6mung innerhalb derselben Kontur sein. 

So gelangen wir zu der in Figur 6a dargestellten Gesamtstruktur des Wir- 
belsystems. Aus dem Unendlichen austrittsseitig laufen unendlich diinne 
Wirbelschichten an die Profilaustrittskanten heran, zerteilen sich hier und 
folgen den beiden Schaufeloberflachen. Langs diesen werden die Wirbelschich- 
ten immer schwacher, da von ihnen aus stetig verteilte Wirbel ins Innere der 
Stré6mung abzweigen und schliesslich in Abstrémrichtung wieder ins Unendliche 
verlaufen. Diese letzteren sind die eigentlichen Sekundirwirbel. Beim Durch- 
stromen des Schaufelkanals wird die Sekundarwirbelstarke immer grdésser, 
doch liegt keine Entstehung der Wirbelstarke im Inneren der Str6mung vor — 
das ware ja durch den dritten Helmholtzschen Wirbelsatz ausgeschlossen -, 
sondern die Wirbelfaden dringen, wie beschrieben, von den lings den Schaufel- 
flachen verlaufenden Wirbelschichten her in den Raum ein. Die Existenz dieser 
Wirbelschichten und ihre Veranderlichkeit langs der Schaufelflache ist gegeben 
durch das Gleiten der in Str6mungsrichtung sich verstarkenden Sekundarbewe- 
gung an dieser Flache. 

Die gesamte Zirkulation des Sekundarwirbelsystems nach Austritt aus der 
Schaufelung ist, wie gezeigt wurde, Null. Die Trennflachen, die itbrigens instabil 
sind, werden sich in Wirklichkeit schliesslich zufolge der Zahigkeit in der Stré- 
mung auflésen, womit aber gerade wieder die Sekundarbewegung «ausgeldscht » 
wird. So schliesst bei der zihigkeitsbehafteten Fliissigkeit das ganze Wirbel- 
system im Endlichen in sich zuriick. 
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Die Starke der an die Austrittskanten anschliessenden Wirbel hat keinen 
unmittelbaren Zusammenhang mit der Abnahme der Zirkulation gegen die 
Schaufelenden zu. Sgurre und Winter?) haben unter vereinfachenden An- 
nahmen (kleine Kanalbreite, keine Beschleunigung oder Verzégerung, Vernach- 
lassigung der Reibung) die Stirke des tiber den Querschnitt stetig verteilten 


Fig. 7 
Kontrollkontur zur Bestimmung der Wirbelstarke der Trennflache. 


Sekundarwirbels berechnet und fanden 


Ow 
2, (3) 


wo # der Umlenkungswinkel ist und 0w/0z der Ortliche Geschwindigkeitsgra- 
dient in der Grenzschicht an der ebenen Begrenzungswand. Die Theorie macht 
die Annahme, dass 0w/0z nur vom Wandabstand abhange. Berechnungen unter 
allgemeineren Voraussetzungen werden bei HAWTHORNE?) durchgefiihrt, doch 
k6nnen alle diese Rechnungen nur als erste Naherung betrachtet werden. Ein- 
fache, allgemein und streng giiltige Aussagen lassen sich jedenfalls nicht machen. 


5. Fir verinderliche Zirkulation entworfene Schaufelung 


Oft miissen Schaufelungen von Turbomaschinen so entworfen werden, dass 
selbst in reibungsfreier Strémung keine Konstanz der Zirkulation langs der 
Schaufel besteht. Im geraden Schaufelgitter, wo die Verhaltnisse leichter tiber- 
blickbar sind, lasst sich dieser Fall dadurch reproduzieren, dass man die Schau- 
feln hier, wo ohne Verdrehung eine konstante Zirkulation erreicht wird, ver- 
dreht. Figur 8 zeigt ein solches Gitter, bei dem, wie der Grundriss 0 zeigt, die 
Ablenkung und damit die Zirkulation gegen oben kleiner wird. Das Stromungs- 


1) H. B. Squire und K. G. Winter, The Secondary Flow in a Cascade of Airfoils in Non- 


uniform Stream, J. aeron. Sci. 18, Nr. 4, 271 (1951). , 
2) W. R. Hawrnorne, Secondary Circulation in Fluid Flow, Proc. Roy. Soc. London [A] 206, 


374 (1951). 
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mittel sei reibungsfrei, und es herrsche vor dem Gitter Potentialstrémung. 
Betrachten wir ein Fliissigkeitsteilchen wahrend seines Durchtrittes durch das 
Gitter, so folgt aus dem Thomsonschen Satz von der Erhaltung der Zirkulation, 
dass es das Gitter drehungsfrei verlassen muss, da es drehungsfrei eingetreten 
ist, dass also mithin die Str6mung auch nach dem Gitter eine Potentialstré- 
mung ist. 


o> 


ZY 


Fig. 8 
Gitter mit variabler Zirkulation und wirbelfreier Anstr6mung. 


Der Seitenriss (Figur 8c) ist ein Schnitt in der Ebene E (siehe Grundriss 6) 
durch den aus einem Schaufelkanal austretenden Stromfaden. Durch die ver- 
schiedene Lange der Pfeile 7 und 2 ist angedeutet, dass die entsprechende 
Geschwindigkeitskomponente gemass der unterschiedlichen Ablenkung oben 
kleiner ist als unten. Dieser Geschwindigkeitsunterschied, multipliziert mit der 
Schaufelteilung, ist gleich dem Unterschied der Zirkulation des untersten und 
obersten Schaufelprofils. Nun muss anderseits die Zirkulation um den ganzen 
Querschnitt des Stromfadens Null sein, da ja Potentialstrémung vorliegt. 
Daher miissen an den beiden anderen Begrenzungsflachen des Stromfadens 
Quergeschwindigkeiten auftreten, gekennzeichnet durch Pfeile 3 und 4, derart, 
dass ihre Wegintegrale entgegengesetzt gleich dem obgenannten Produkt aus 
Geschwindigkeitsdifferenz und Teilung sind. So werden die gegenseitigen Be- 
rihrungsflachen der aus den einzelnen Schaufelkandlen austretenden Strom- 
faden zu Unstetigkeitsflachen, und zwar ist nach der eben durchgefiihrten 
Uberlegung die Zirkulation einer solchen Unstetigkeitsfliche gleich dem ge- 
samten an einer Schaufel auftretenden Zirkulationsunterschied. 
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Hier haben wir also den Fall vor uns, wo in exakter Analogie zum Einzelfliigel 
von der Austrittskante eine Wirbelschicht abgeht, deren Starke durch den Zir- 
kulationsunterschied an der Schaufel gegeben ist. Figur 8a zeigt den Verlauf 
der Wirbellinien im Aufriss. Die Stromlinien 5 und 6 biegen in der gezeigten 
Weise ab und bewirken damit auch ein entsprechendes Abbiegen der lings der 
Schaufel verlaufenden Wirbellinien. Diese letzteren miinden damit zum Teil in 
die Austrittskante ein und laufen in der Wirbelschicht stromabwiarts fort. Die 
Existenz dieses Wirbelsystems wurde hergeleitet ausgehend vom Thomsonschen 
Satz, der eine unmittelbare Folge der dynamischen Bewegungsgleichung fiir 
zahigkeitsfreie Fliissigkeit ist. 

Man beachte, dass die Ausfiihrungen dieses Abschnittes zugleich fiir den 
Fall eines Spaltes zwischen Schaufelende und Begrenzungswand gelten. Ein 
solcher Spalt ist ja grundsatzlich nichts anderes als ein « Stiick Schaufel mit der 
Zirkulation Null». 


6. Schaufelung bei wirbeldurchsetzter Anstrémung 


Figur 9 zeigt das Gegenbeispiel der oben behandelten Str6mung. Einem 
Schaufelgitter stromt das reibungsfrei vorausgesetzte Str6mungsmittel mit 
langs der Schaufelhéhe variierender Richtung zu, zum Beispiel unten gemdss 
den gestrichelten Stromlinien 7, in der Mitte entsprechend den ausgezogenen 
Linien und oben gemass den gestrichelten Stromlinien 2, wobei die Grésse der 


Fig. 9 
Gitter mit variabler Zirkulation und wirbeliger Anstromung. 


Zustromgeschwindigkeit tiber die Schaufelhdhe konstant sei. Die Schaufelung 
sei so gestaltet, dass die Austrittsrichtung — abgesehen von schwachen, durch 
allfallige Wirbelbewegungen bedingten Quergeschwindigkeiten — konstant sel. 

Die Anderung des statischen Druckes beim Durchstrémen des Gitters ist fiir 
alle Fliissigkeitsteilchen gleich gross, da ja in den Raumen vor und nach dem 
Gitter ausgeglichener Druck herrscht. Nach BERNOULLI ist daher die Abstrom- 
geschwindigkeit fiir alle Teilchen gleich, weil dies auch fiir die Zustrd6mgeschwin- 


ZAMP Iv/20 
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digkeit gilt. Konnte das Gitter eine vollkommen unveranderliche Austrittsrich- 
tung erzwingen, so lage am Austritt eine Potentialstro6mung vor, und die Schau- 
felung wiirde folglich die zustr6menden Wirbel «aufschlucken ». Nach demThom- 
sonschen Satz kénnen aber die einzelnen Fliissigkeitsteilchen bei Reibungsfrei- 
heit ihre Drehung nicht verlieren, so dass die Wirbelfaden mit unveranderlicher 
Stirke das ganze Gitter durchsetzen (vgl. die ausgezogenen Linien 3). Die 
Stromfaden, die aus den einzelnen Schaufelkandlen austreten, sind also wirbel- 
durchsetzt und grenzen durch Trennflachen 4 aneinander an. Durch entspre- 
chende Schaufelgestaltung kann erreicht werden, dass die Zirkulation einer 
solchen Trennflache entgegengesetzt gleich derjenigen des aus einem Schaufel- 
kanal austretenden Stromfadens ist. Die resultierende Zirkulation der austreten- 
den Strémung ist dann Null, obwohl sie keine Potentialstromung ist. Mit Riick- 
sicht auf die Verluste in Schaufelungen ist es wichtig, aus dieser Uberlegung zu 
erkennen, dass eine Schaufelreihe einer wirbeldurchsetzten Strémung der hier 
angenommenen Art ihre Wirbel nicht einfach auf reversiblem Wege entziehen 
kann, sondern sie kann nur Gegenwirbel in die Strémung abgehen lassen, 
welche durch Reibungseffekte die urspriinglichen Wirbel zum Verschwinden 
bringen. Die Zirkulation der Trennflache ist hier wiederum gleich dem Zirku- 
lationsunterschied der beiden Endprofile der Schaufel, und die Gestalt des ent- 
sprechenden Wirbellinienbildes ist grundsatzlich die in Figur 8a angegebene. 


7. Wirbelverluste 


Trotz der Kompliziertheit der Vorginge kann man doch qualitative Aus- 
sagen tiber die Wirbelverluste machen. Betrachten wir zuerst den Sekundar- 
wirbel. Fiir ein nur ablenkendes (nicht die Geschwindigkeit 4nderndes) Gitter 
fanden SguirRE und WINTER Gleichung (3) als Naherungsausdruck fiir die Wir- 
belstarke ¢. Figur 10a veranschaulicht den senkrechten Querschnitt, durch den 
aus einem Schaufelkanal austretenden Stromfaden, dessen Héhe gleich 7 und 
dessen Breite gleich a ist. Ist die Gestalt des Geschwindigkeitsprofils gegeben, 
die Grenzschichtdicke also proportional /, so ist 0w/0z umgekehrt proportional / 
und proportional der Geschwindigkeit w,) in der Mitte des Stromfadens. Also 
wird 


- Wo 
C~ oS, (4) 


Die langs der Kontur A BCD gebildete Zirkulation [" ist das tiber die entspre- 
chende Flache erstreckte Integral von €, so dass die Proportionalitat 


Prall~buya (5) 


gilt. Anderseits ist J’ proportional den linearen Abmessungen und proportional 
den durch die Wirbelverteilung induzierten Geschwindigkeiten, fiir die ein in 
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irgendeinem festen Punkt vorhandener Wert w’ als Charakteristikum eingesetzt 
werden kann. Es ist daher ~ a w’ und somit aus (5) 


aw ~dwa. (6) 


Schliesslich ist der durch die Wirbelbewegung gegebene relative Verlust Z (kine- 
tische Energie der Wirbelbewegung) proportional (w’/w 9)”, also nach (6) 


Zo (7) 


Dieses Resultat gilt zwar nur unter den vereinfachenden Voraussetzungen 
der Theorie von SQuriRE und Wryrer, doch kann daraus der ganz allgemeine 


Fig. 10 
Kontrollkonturen zur Bestimmung des Energieverlustes durch Sekundarstrémung. 
a Fir extrem schlanke Schaufeln; 
b fiir extrem gedrungene Schaufeln. 


Schluss gezogen werden, dass der durch den Sekundarwirbel bedingte Verlust 
mit zunehmendem Ablenkungswinkel sehr stark ansteigt. Eine allgemeine Er- 
kenntnis lasst sich ferner gewinnen tiber den Einfluss des Erstreckungsverhalt- 
nisses der Schaufeln. Wir nehmen an, die Gestalt des Gitters im Schnitt parallel 
zu den Begrenzungsebenen sei gegeben, hingegen werde das Verhaltnis der 
Schaufelhéhe zur Sehnenlange (Erstreckungsverhaltnis) varitert. Die Gestalt 
des Geschwindigkeitsprofils sei ebenfalls vorgegeben. 

Figur 10a und 6 sind senkrechte Schnitte durch aus einem Schaufelkanal 
austretende Stromfaden, und zwar entspricht Figur 10a extrem schlanken, 
Figur 10 extrem gedrungenen Schaufeln. Betrachten wir in beiden Fallen die 
Zirkulation I’ langs der Kontrollkontur ABCD, so ist sie 


a B G D A 
Patal= [wf ds+ fwids+ | wjas+ | mds, (8) 
A PB C D 


wobei w’ die zur Kontur tangentiale Komponente der Sekundargeschwindig- 
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keit w’ ist, s der Integrationsweg und ¢ der iiber die Flache ABCD gebildete 
Mittelwert der Wirbelstarke. 

Beschranken wir die Betrachtung zuerst auf den Fall nach Figur 10a, so 
kénnen im Ausdruck (8) die Integrale von B nach C und von D nach A nahe- 
rungsweise vernachlassigt werden, da hier der Integrationsweg kurz und gleich- 
zeitig der Integrand klein ist. Die beiden restlichen Integrale ergeben zusammen 
1 @,, wo @, der Mittelwert von w iiber die betrachteten Strecken ist, also 


abla (9) 


Nun setzen wir 


aeSing ge’ my eo Aewee (10) 


wo W = w/w, z* = 2/l. Bei gegebener Gestalt des Geschwindigkeitsprofils ist 
dann 0W/dz* von den absoluten Abmessungen des Systems unabhangig. Der 
Ansatz (10) ist eine Verallgemeinerung der Formel von SQUIRE und WINTER 
und diirfte in dieser Form wohl unter sehr allgemeinen Voraussetzungen min- 
destens naherungsweise zutreffen. Mit (10) ist (9) in der Form 


aw OW ee 
ae err (11) 


darstellbar. Fiir eine gegebene Gestalt des Gitters und des Geschwindigkeits- 
profils wird dann der Verlust 


Wy 2 a\2 
ibaa post orl (12) 
oder, da anderseits die Sehnenlange s des Profils proportional a ist 
Z=K(-+). (13) 


Liegt umgekehrt der Extremfall des ganz gedrungenen Gitters nach Figur 10d 
vor, so verschwinden praktisch die Integrale von A bis B und C bis D, und die 
beiden anderen ergeben zusammen 2 a @;. Folglich tritt an die Stelle von (9) 


bo] sy 


—al=2aw,, (14) 


wahrend (11) zu ersetzen ist durch 


ow —- 
a Wo Oz* ~a Ws . (15) 
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Damit wird aber das fiir den Verlust massgebende Verhaltnis W,/W») unabhangig 
vom Erstreckungsverhiltnis, das heisst, wir haben ein konstantes Z: 


VANS (16) 


Uberblickt man (13) und (16), so erkennt man, dass der Verlust Z vom Er- 
streckungsverhdltnis s/l abhangt, gemass einer Funktion, die durch die Kurve 
der Figur 11 qualitativ dargestellt ist. Fiir sehr kleine s/J, also sehr schlanke 


Fig. 11 
Verlauf der Wirbelverluste in’ Funktion des Erstreckungsverhaltnisses. 


Schaufeln, hat die Kurve den Charakter einer quadratischen Parabel, geht aber 
fiir sehr grosse s//, also gedrungene Schaufeln, asymptotisch in den Wert Zax 
liber. 

Bekanntlich bestehen widersprechende Behauptungen dartiber, ob das Er- 
streckungsverh4ltnis den Wirkungsgrad einer Schaufelung massgebend beein- 
flusse, wobei sich diese Behauptungen alle auf Messergebnisse stiitzen. Beachtet 
man den Verlauf der Kurve der Figur 11 und den durch (7) wiedergegebenen 
Zusammenhang, so wird dies verstandlich. Es gibt Falle, wo der Einfluss von 
s/l gering oder der ganze durch Sekundarwirbel bedingte Verlust nach (7) sehr 
klein ist (zum Beispiel Axialverdichter, Kaplan-Turbine) und andere, wo das 
Gegenteil zutrifft. 

Die Uberlegungen, die zur Kurve der Figur 11 fiihren, gelten tibrigens auch 
fiir Verluste, wie sie durch die in Abschnitt 5 und 6 behandelten Strémungs- 
formen bedingt sind. In der Tat kénnen die Betrachtungen tiber die Zirkulation 
IT’ ebensowohl an Kontrollkonturen durchgefiihrt werden, welche den ganzen 
aus einem Kanal austretenden Stromfaden umspannen. Die Relation (10) be- 
steht allerdings nicht mehr, doch tritt an ihre Stelle eine der Form 


pes (17) 
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wo g eine Grdsse ist, die nur von der gesamten Anderung der Zustrémrichtung 
bzw. der Gittergestalt lings der Schaufelhéhe abhangt. Wichtig ist an (17), dass 
lim Nenner steht, was darauf beruht, dass die Wirbelstarke proportional der 
Anderung der Zirkulation (um die Schaufel) pro Langeneinheit ist. Von (17) 
aus gelangt man zu (13) und (16) genau gleich wie von (10) aus. 

Beachtlich ist, dass man mit dieser Uberlegung auch den in Abschnitt 5 
behandelten Fall mit umfasst, obwohl dort itiberhaupt keine iiber den Querschnitt 
stetig verteilte Wirbelstarke vorhanden ist. Wiirden namlich die durch die 
Pfeile 3 und 4 in Figur 8c angedeuteten Querbewegungen nicht auftreten, so 
lage eine wirbeldurchsetzte Strémung vor, und die genannten Querbewegungen 
bestimmen eine genau entgegengesetzte Wirbelverteilung ¢, welche gerade die 
Wirbelfreiheit wiederherstellt. Es ist diese Wirbelverteilung ¢, auf die obige 
Uberlegungen anzuwenden sind, denn diese Querbewegungen werden in nach- 
folgenden Schaufelungen nicht ausnutzbar sein, sondern ihre kinetische Energie 
wird durch Reibung aufgezehrt werden. 

Selbst die Beziehung (7) hat ein Analogon im Falle der in Abschnitt 5 und 6 
behandelten grenzschichtfreien Str6mungen mit langs der Schaufel variierender 
Zirkulation. In der Tat sind ja die Quergeschwindigkeiten proportional der 
Verdnderlichkeit der Zirkulation, wie aus dem dritten Helmholtzschen Wirbel- 


satz folgt. Damit wird aber der Energieverlust proportional dem Quadrat der — 


Veranderlichkeit der Zirkulation. Hier mag es zunachst scheinen, als sei eine 
physikalische Schlussfolgerung aus einem rein mathematischen Satz gezogen 
worden, doch trifft dies nicht zu, denn die Uberlegung wird nur dadurch még- 
lich, dass wir die Struktur des Wirbelsystems kennen, und diese wurde aus 
physikalischen Gesetzen erschlossen. Die quadratische Abhangigkeit der Ver- 
luste von der Variation der Zirkulation langs der Schaufel bedeutet, auf den 
Turbomaschinenbau iibertragen, dass fehlende oder falsche Verdrehung der 
Schaufeln den Wirkungsgrad nur ganz unmerklich beeintrachtigen, sofern der 
Fehler nicht allzu gross ist. Nach der Kurve der Figur 11 trifft dies in besonders 
hohem Masse fiir schlanke Schaufeln zu. Zusitzliche Verluste entstehen in 
solchen Fallen, wenn iiberhaupt, so in erster Linie durch ungiinstiges Anstr6men 
einzelner Schaufelprofile und damit verbundene Strémungsablésungen. 


8. Ubertragung auf die Turbomaschine 


Die Verhaltnisse in der Schaufelung einer Turbomaschine weichen in dreierlei 
Hinsicht von denen im geraden Schaufelgitter ab. 

Erstens ist die Str6mungsform, die beim Gitter als idealer Grundfall gilt — 
Richtung und Grésse der Geschwindigkeit sowie Druck sind tiberall langs der 
Schaufelhéhe konstant —, kein méglicher Gleichgewichtszustand mehr. An ihre 
Stelle tritt die im englischen Sprachgebrauch als «free vortex flow» bezeichnete 
Strémung : Axialkomponente tiber die Schaufelhéhe konstant, Umfangskompo- 


— 
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nente c, langs des Radius 7 variierend gemass + c,, = const. Abweichungen 
gegenttber den oben durchgefiihrten Uberlegungen entstehen dadurch nicht, 
denn man hat nur die Strémung 7 c,, = const als idealen Grundfall zu nehmen 
und Abweichungen hiervor in genau gleicher Weise zu betrachten, wie das oben 
fiir das gerade Gitter geschah. 

Zweitens tritt als weiterer Effekt hinzu, dass die Grenzschichten an den 
Laufschaufeln unter dem Einfluss der Fliehkraft nach aussen abgedrangt wer- 
den. Damit tiberlagert sich den bisherigen Wirbelstrémungen noch eine weitere, 
deren Mechanismus demjenigen der gewéhnlichen Sekundarstrémung analog 
ist. Die Bedeutung dieses Effektes wird iibrigens haufig tiberschatzt, denn die 
Spuren, welche die Str6mung oft an bestaubten Laufschaufeln hinterlisst, 
zeigen, dass die Grenzschichten nur verhaltnismassig wenig auszentrifugiert 
werden. 

Eine dritte Abweichung — die weitaus schwerwiegendste — besteht darin, 
dass die hier durchweg vorausgesetzte stationare Str6mung bei der Turbo- 
maschine héchstens im ersten Schaufelkranz verwirklicht ist. Schon die dem 
zweiten Schaufelkranz zufliessende Strémung hat (relativ zu diesem) einen 
komplizierten periodischen instationdaren Charakter. Was die Turbomaschinen- 
theorie bis heute tut — ohne es ausdriticklich zu sagen, meist wohl ohne sich tiber- 
haupt dariiber im klaren zu sein — ist folgendes: Man betrachtet diese ganzen 
instationaren Schwankungen als eine Art grobe Turbulenz, welche sich einer 
stationdaren Grundstrémung iiberlagert. Wenn der Charakter dieser stationaren 
Grundstroémung durch diese «geordnete Grobturbulenz» nicht massgebend be- 
einflusst wird — das hat man bisher angenommen -, dann lassen sich alle oben 
durchgefiihrten Uberlegungen ohne weiteres auf die Turbomaschine tibertragen. 
Ob aber diese Annahme zulassig ist, bleibt dahingestellt. Dafiir spricht lediglich 
die Tatsache, dass grobe, unerklarliche Abweichungen zwischen der Erfahrung 
und sonst gut fundierten Theorien im Turbomaschinenbau kaum aufgetreten 


sind. 


Summary 


When a fluid is flowing past rows of deviating blades, such as used in turbo 
machines, different vortex systems will arise. These vortex systems are largely 
uninvestigated and so far there has not even existed any survey of the pheno- 
mena to be expected. An attempt is made to give such a survey. 


(Eingegangen: 10. Februar 1953.) 
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Optimale Rechengenauigkeit bei Rechenanlagen mit gleitendem 
Komma 


Von Kiaus SAMELSON und FRIEDRICH L. BAUER, Miinchen’) 


§ 1. Problemstellung 


Fiir die Darstellung von Zahlen in programmgesteuerten Rechenautomaten 
sind zwei Varianten unter den Stichworten «festes Komma» und «gleitendes 
Komma» bekannt?). ' 

Beim Rechnen mit gleitendem Komma werden die Zahlen in der Maschine 
selbst halblogarithmisch, das heisst in der Form ¥ 2", x<1, dargestellt*). Die 
bei festem Komma notwendigen umstandlichen Vorbereitungen zur Feststellung 
des angemessenen Zahlenbereichs oder Programmierungsmassnahmen, die ein 
Uberlaufen des Rechenwerks verhindern, fallen dabei weg. Die Maschine errech net 
automatisch den jeweils giiltigen Skalenfaktor, der sich im Exponenten ausdriickt 
und iiblicherweise nach jeder arithmetischen Operation durch die Forderung 
1/2 < * <1 fiir die Mantisse * neu festgelegt wird. 

Bekanntlich werden dadurch haufig, etwa bei der Subtraktion nahezu gleicher 
Zahlen, die letzten, mit Ungenauigkeit (zum Beispiel Rundungsfehlern) behaf- 
teten Stellen unter Vortaéuschung einer nicht vorhandenen Genauigkeit weit nach 
links gertickt. Um dies zu verhindern, muss man die Forderung 1/2 S ~» fallen- 
lassen. Vielmehr hat die Maschine selbst (gegebenenfalls nach Anweisung durch 
das Programm) die Genauigkeit des Ergebnisses jeder arithmetischen Operation 
abzuschatzen und die Mantisse entsprechend zu verschieben. Dadurch wird hier 
der Exponent des Ergebnisses festgelegt. Die erstmalige Fixierung des Dualexpo- 
nenten bei der Zahleneingabe ist Sache der Umrechnungsroutine und wird in § 5 
dieser Arbeit behandelt. 


§ 2. Diskussion des Abbrechfehlers 


Ausdriicklich betont werden soll, dass die von der Maschine durchzufiihrende 
Abschatzung sich muy auf die einzelne Rechenoperation bezieht. Die Fortpflan- 
zung von Ungenauigkeiten in den Eingabewerten und von internen Rundungs- 
fehlern durch langere Serien von Rechenoperationen hindurch kann wohl kaum 
von der Maschine selbst kontrolliert werden und wird daher hier nicht diskutiert. 

Dementsprechend wird bei der Untersuchung der arithmetischen Operation 
die Annahme gemacht, dass die Mantissen der beiden der Maschine jeweils zur 
Verfiigung stehenden Operanden nach einem der in RSS. § 3.6, angegebenen Ver- 
fahren richtig gerundet sind und infolgedessen einen Fehler tragen, der etwa 
dem Verfahren b maximal die halbe Einheit der letzten angegebenen Stelle bei 


1) Mathematische Arbeitsgruppe ftir Rechenanlagen der Technischen Hochschule Miinchen. 
*) Vgl. zum Beispiel den Bericht von H. RurisHauser, A. SPEISER und E. STIEFEL, ZAMP 1, 
277, 339 (1950); 2, 1, 63 (1951). (Zitiert: RSS.) 
3) Wir beziehen uns der Einfachheit halber stets auf das reine Dualsystem. 
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betragt. Das formale Ergebnis der Rechenoperation ist dann mit einem gewissen 
Fehler / der Mantisse behaftet. Besitzt es mehr Stellen als die Operanden, so 
soll die Gréssenordnung des Fehlers von der Maschine zur Bestimmung der Zahl 
der brauchbaren Stellen beniitzt werden. 

Es ist selbstverstandlich, dass der zusdtzliche Rundungsfehler /, den das 
Abbrechen des formalen Ergebnisses mit sich bringt, von der gleichen Gréssen- 
ordnung ist wie der Fehler F,. Das Ergebnis tragt also einen Gesamtfehler, der 
stets grosser als F, ist und maximal etwa das Doppelte davon betragt. Es kann 
aber von per Maschine bei der nachsten Rechenoperation nur als richtig gerundet 
behandelt werden. Ein allmahliches Akkumulieren des Gesamtfehlers, wie es sich 
auch hier ergibt, kann demnach prinzipiell nicht verhindert werden, wenn man 
nicht stets so abbricht, dass F, gegen F, zu vernachlassigen ist. Dies wiirde jedoch 
den Verlust einer ganzen Reihe brauchbarer Stellen bei jeder Operation bedeuten, 
wie er zum Beispiel bei der Multiplikation mit festem Komma auftritt, wenn beide 
Faktoren klein gegen 1 sind. 


§ 3. Die arithmetischen Operationen 


Der obenerwahnte Fehler Ff, ist nun fiir die Grundrechenoperationen ‘with 
the machine’s eye-view”’ abzuschatzen. Dabei ist an eine Parallelmaschine mit 
N Mantissenstellen gedacht. Die Operanden seien mit 


Aa A) 2 und 25 = (et Ay) 2M 


bezeichnet. * und y sind die als Rechengréssen zur Verfiigung stehenden Man- 
tissen, Ax und Ay die Rundungsfehler, deren Betrage voraussetzungsgemass 
jeweils héchstens den Wert 2~(V+1) haben. 


a) Addition und Subtraktion 


Vor Beginn einer Addition oder Subtraktion sind nétigenfalls die Exponenten 
anzugleichen, selbstverstandlich durch Erhéhung des kleineven mit gleichzeitiger 
Rechtsverschiebung der zugehérigen Mantisse und Rundung auf N Stellen, zum 
Beispiel fiiyv n< m 

HOM 4 y 2M — (x ZM—M) 2M 4 y QM 
Die bei dem iiblichen Verfahren eventuell anschliessende Linksverschiebung der 
Summe unter Nachholung von Nullen ist zu unterlassen. Der Gesamtfehler 
betragt maximal eine Einheit der letzten Stelle, ist also héchstens gleich 2 Fy. 


b) Multiplikation 
Fiir den Fehler F, gilt bis auf Gréssen zweiter Ordnung 
R= |x Ay + y Ax| S (/rl + yl) 2°84) S Max(4]; |y|) 2-%. 
Sei etwa |x| < |y| und 2-(#+)) < |y|< 2-*, das heisst, |y| habe genau # Nullen 
“hinter dem Komma. Dann ist RK < 2~\V+*), Abzubrechen ist dementsprechend 
hinter der (N + k — 1)-ten Stelle; der Gesamtfehler ist kleiner als 2 F. 
Die (N + k — 1)-te Stelle befindet sich im allgemeinen in MR. Es sind also 


k — 1 Linksverschiebungen aus MR nach AC durchzuftihren. Dies ist stets még- 
lich, da voraussetzungsgemass das Produkt mindestens 2 k Nullen hinter dem 


Komma enthalt. 
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c) Division 


Fiir den Quotienten z,/z, muss |%|< |y| vorausgesetzt werden, damit die — 
Durchfiihrbarkeit der Division gewdhrleistet ist. Dies ist gegebenenfalls vorher 
durch Rechtsverschiebung der Zaéhleymantisse zu erzwingen. 

Fiir F, ergibt sich in erster Ordnung 


Ta Le <i ( ) wey = 2 ow 

fis ried #4n| S75 mail : = TA 7 

Sei wieder 2-@+1) < |y|< 2-*, dann ist 1/|y| S 2*+1!, also A S 2-(N-F-)), 
Abzubrechen ist also hinter der (N — k — 2)-ten Stelle, die in AC liegt. Demnach 
sind k + 2 Rechtsverschiebungen durchzufiihren. Der Gesamtfehler betragt wie- 
der héchstens 2 Fy. 

Die angegebenen Abschatzungen fiir F, sind unter Umstanden um einen 
Faktor 2 zu grob, soweit sie von der Zahl k der Nullen von || hinter dem Komma 
abhangen, fiir |¥| <|y| sogar nahezu um einen weiteren Faktor 2. Da tiber- 
dies der Mittelwert des Fehlerbetrags sicher merklich kleiner als die angegebenen 
Grenzen ist, k6nnte man daran denken, in der Regel bei der Multiplikation k 
Linksverschiebungen, bei der Division & + 1 oder gar nur k Rechtsverschiebungen 
des Zwischenergebnisses durchzufiihren und zusatzliche Rechtsverschiebungen 
nur bedarfsweise einzuprogrammieren. Dies ware auch deshalb zweckmiassig, weil 
der Fehler F auf die Halfte zuriickgeht, wenn einer der beiden Operanden (bei 
der Addition derjenige mit dem grésseren Exponenten) exakt ist. 

Sind beide Operanden exakt, so kann angenommen werden, dass sie unter 
Ausniitzung der vollen Stellenzahl, also ohne Nullen hinter dem Komma, vor- 
liegen. Das Ergebnis steht dann von selbst richtig. 


§ 4. Technische Durchfithrung 


Abschliessend soll diskutiert werden, durch welche Massnahmen die Maschine 
nach jeder Rechenoperation die richtige Lage der Mantisse und damit den Expo- 
nenten automatisch bestimmen kann. 

Bei der Addition und Subtraktion ist die technische Durchfiihrung nach 3a) 
im grundsatzlichen evident. Der Exponent wird Max(n; m). 

Die Multiplikation und Division ist abhangig von der Anzahl k der Nullen 
nach dem Komma im Betrag von y, das heisst der grésseren der beiden Faktor- 
mantissen bzw. der Mantisse.des Nenners. Eine einheitliche Behandlung ist még - 
lich, wenn man stets vor Beginn der Rechenoperation den Fall k = 0 erzwingt. 
Man hat dazu so lange Linksverschiebungen (bei der Division in MD, bei der Multi- 
plikation in MD oder MR) durchzufiihren, bis die 2°- und die 2-!-Stelle entgegen- 
gesetzte Werte annehmen?). Die dabei nachgeholten Nullen beeinflussen nur 
Stellen des Ergebnisses, die beim Abbrechen wieder wegfallen. Abschliessende 
Rechtsverschiebungen kénnten nach der im Schlussabsatz des § 3 gefiihrten Dis- 
kussion unter Umstanden wegfallen. 

Eine einheitliche Behandlung der Multiplikation und Division soll kurz skizziert 
werden: 

Beide Operanden werden gleichzeitig so lange nach links verschoben, bis sie 
am Komma auflaufen. Zur Erm6églichung der Division ist unter Umstinden eine 
Rechtsverschiebung in AC vorzunehmen. Die an dem kleineren Operanden durch- 


1) Nachdem bei der Division AC gegebenenfalls vorher auf die Stellenzahl von MD nach 
rechts verschoben worden ist, wie in § 8c) erwahnt. 
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geftihrten tiberzahligen Linksverschiebungen werden durch eine gleiche Anzahl 
Rechtsverschiebungen des Ergebnisses kompensiert. Der Mehraufwand an Schalt- 
elementen ist gering gegentiber dem gesamten Aufwand fiir eine Rechenmaschine 
mit gleitendem Komma in der iiblichen Bauweise. 


§5. Konvertierung halblogarithmisch dargestellter Zahlen vom 
Dezimal- ins Dualsystem 


Eine Zahl s laute in dezimaler, halblogarithmischer Darstellung 


Sdez = Paez LO: (Paez a 0,PiPops- : PL) 


Alle L Stellen sollen mit giiltigen Ziffern ausgefiillt sein, also 1/10 < p,,, < 1, 
falls sg,, exakt1) oder mit mindestens L Stellen bekannt ist. Im letzteren Fall 
betragt der Rundungsfehler héchstens eine Einheit der letzten Stelle. Falls jedoch 
Saez» Sei es etwa als Naherungsbruch einer Irrationalzahl oder als fehlerbehafteter 
Messwert einer physikalischen Grésse, mit weniger als L giiltigen Ziffern vorliegt, 
so soll, in Ubereinstimmung mit dem Grundsatz von § 1, die letzte gultige Ziffer 
an die L-te Stelle zu legen kommen?). 

Durch diese Forderung der «optimalen Genauigkeit » legen wir nun auch die 
Mantisse ug,,,, und den Exponenten & von s nach der Konvertierung 


eas pee t— 
Sdez = Paez 10* = Uaual Dees Sdual 
fest, das heisst, wir bestimmen & so, dass fiir 


10? 
Udual = (Pace 2k ), ~~, Pauat lauat 
Ua 


die Bedingung 1/2 S tg, = (10*/2") a4, < 1 erfiillt wird. 

Die Umrechnungsvorschrift lautet also: 

1. Konvertiere fg,, in Pgyqz, 1m wesentlichen nach bekannten Methoden (SSR. 

Se); 

2. konvertiere 10* in t4.4; 2° mit 1/2 = 74,9, < 1; 
3. bilde Pauat (tguat 2) = Uquat Ze 
Squat = “auar 2” ist die Dualdarstellung von s. 

Zu 1: Die Umrechnung einer in L Tetraden direkt verschliisselten Dezimal- 
mantisse pg,, in eine 4//-stellige Dualmantisse geschieht bekanntlich folgender- 
massen: 

M Loe 
Jo= 9, G=10g,, +16" p, (R= 1,2,..,£L), Pana = IL: Tem )* 
L Tetraden ergeben, dual umgerechnet, rund #log(10%) w 3,32 L Dualziffern. 
Ublicherweise hat man L = M und fiillt ~ 0,68 M Stellen der insgesamt N =4M 
Dualstellen nicht aus. Wenn man jedoch die Konvertierung direkt an die tetra- 
denweise Eingabe anschliesst, kann man eine fast vollstandige Ausniitzung der 
verfiigbaren Dualstellen erzielen, zum Beispiel mit L = 12 bei M= 10 (N = 40). 
Zu 2. Die Umrechung von 10* in ¢4,,4, 2" verlauft ebenfalls rekursiv, wobei jetzt 
das Exponentenrechenwerk herangezogen wird. Fiir 7 > 0 hat man 
1 5 : ; 
Yesrey 24, Va ee, 23 yaa,» Vi = (10) auat = banat 2* - 
1) Exakt (hinsichtlich der Stellenzahl L) soll heissen, dass sg,, eine ganze Zahl < 10L, divi- 


diert durch eine beliebige Zehnerpotenz, darstellt. i a 
2) Oder an eine andere, aber feste Stelle L’ < L, wenn in einem Rechengang grundsatzlich nur 


Zahlen mit weniger als L’ Stellen verarbeitet werden sollen. 
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Die Multiplikation mit 5/4 = 2° + 2-2>1, die nicht als Routinemultiplikation 
ausgefiihrt werden kann, lauft in einem einfachen Zyklus von Additionen und 
Rechtsverschiebungen ab. Dabei erfolgt immer wieder ein Uberlaufen des Man- 
tissenrechenwerks, das durchaus erwiinscht ist. Es fiihrt zur automatischen Kor- 
rektur durch Rechtsverschiebung unter Erhédhung des Exponenten um 1. Man 
hat somit fiir f4,,,, gerade immer die Forderung 1/2 S fgyq;< 1 erfiillt. 

Fiir 7 < 0 ist Payq, 10!*! zu berechnen. 

Verwendet man zweistellige Dezimalexponenten 7 = «10+ B, so bedeutet 
es eine Abkiirzung des Verfahrens, wenn man, mit Hilfe des gespeicherten Wertes 
(102°) 5.47 = (1020/234) 234 = v 234, 10 = 10% 10+ als r* 108. 24% berechnet. Man kann 
dann den Dezimalexponenten direkt tetradisch verschliisselt eingeben und behan- 
deln. Negative Exponenten werden zweckmassig in derselben Form «10 + 8, 
aber mit «< 0, 0< #6 <9 dargestellt, so dass hinsichtlich f keine Fallunter- 
scheidung mehr zu machen ist. Auf die verschiedenen méglichen Varianten des 
Verfahrens soll nicht weiter eingegangen werden. 


§ 6. Riickkonvertierung vom Dual- ins Dezimalsystem 


Die Riickkonvertierung von Sgyq_ = Uguay 2” stellt grundsatzlich eine Um- 

kehrung des Verfahrens von § 5 dar. Die Umrechnungsvorschrift lautet: 

1. Konvertiere 2” in 10*/tg,,4;, 1/2 S tauai << 1; speise 7 aus; 

2. bilde Uguat: fauar = Pauat’ 

3. konvertiere Pgyq_ iN Paez, im wesentlichen wieder nach bekannten Methoden 

(IRS 55 Sed) 

Paez 10* ist die Dezimaldarstellung von s = sgyq_- 

Zu 1: Fiir die Umrechnung von 2?” in 10*/tj,,,, wird die Rekursion von § 5, 1,. 
bentitzt. Sie ist abzubrechen bei der Zahl 2, fiir welche erstmals k = k’ wird"). An- 
schliessend sind noch k — k’ Rechtsverschiebungen von uq,,,, durchzufiihren, 

Zu 2: Die Division der reinen Mantissen, Wq,9;: ¢gyqyz, Kann nur noch als Ope- 
ration mit festem Komma ausgefiihrt werden, da der Exponent bereits abgebaut 
ist. Sie ist stets ausfiihrbar, wenn man (als einfachste Lésung) vor Ablauf des 
Schrittes') durch eine Rechtsverschiebung der Mantisse mit Exponentenanglei- 
chung |Uguaz|< 1/2 erzwingt. 

Auch die Ausgabe zweistelliger Dezimalexponenten ist in sinngemasser Um- 
kehr der Uberlegung am Schluss von § 5 direkt tetradenweise durchfiihrbar. 
Hierfiir ist zunachst die Zehnerstelle 


tire Sees ge ee (k* = k’ mod 34; 0S k” < 34) 
zu berechnen und wq,,, durch r* zu dividieren. Der Rest k” wird mit Hilfe der 
Zehnerrekursion abgebaut und liefert die Einerstelle 6 des Dezimalexponenten. 


Summary 


Digital automatic computing machines with floating binary point are investi- 
gated with respect to their accuracy. It is concluded that special devices can be 
designed without difficulty, which guarantee the highest accuracy obtainable 
with a fixed number of digits. For binary computers, a modification of the usual 
convertion routine is necessary, which is discussed in detail. 

(Eingegangen: 27. Februar 1953.) 


F 1) Praktisch diirfte es wohl am zweckmiassigsten sein, das Anfangsglied yy der Rekursion durch 
Yo = Vg 2~'*! zu ersetzen und nach dem Vorzeichenwechsel im Exponenten abzubrechen. 
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Frihjahrstagung der Schweizerischen Physikalischen Gesellschaft 
vom 2. Mai 1953 in Genf?) 


Berichte tiber angewandte Physik und Mathematik 


Biegeschwingungen verwundener Stabe, von A. TROscH?), New York, 
M. ANLIKER®), Ziirich, und H. Z1ieGrer, ETH., Ziirich. 

Mit Hilfe der programmgesteuerten Rechenmaschine des Instituts fiir ange- 
wandte Mathematik an der ETH. werden die natiirlichen Frequenzen der Biege- 
schwingungen von einseitig eingespannten, geraden Stiben, die im ungespannten 
Zustand wie ein Propellerblatt verdreht und deren Querschnitte ideal schmal 
sind, exakt berechnet. Die Masse pro Langeneinheit, der Verdrehungswinkel pro 
Langeneinheit sowie die Biegesteifigkeiten beziiglich der Hauptachsen werden 
dabei als konstant vorausgesetzt. Eine ausfiihrliche Veréffentlichung ist geplant. 


Projekt einer elektronischen Rechenmaschine am Institut fiir angewandte 
Mathematik der ETH., von A. P. SprEtser’). 

Eine der Aufgaben des Instituts fiir angewandte Mathematik ist die Lésung 
numerischer Aufgaben, welche aus Kreisen der Hochschule und der Industrie ge- 
stellt werden. Zu diesem Zweck wurde 1950 eine programmgesteuerte, mit Relais 
arbeitende Rechenmaschine in Betrieb genommen. Diese ist aber heute trotz 
einem fast durchgehenden Tag- und Nachtbetrieb dem Andrang an Problemen 
nicht mehr gewachsen, so dass viele Aufgaben zurtickgestellt werden miissen. 

Es ist daher ein Projekt einer elektronischen Rechenmaschine ausgearbeitet 
worden, welches sich vom bestehenden Geradt durch héhere Rechengeschwindig- 
keit, gréssere Speicherkapazitat und flexiblere Programmgebung unterscheidet. 
Die projektierte Maschine arbeitet im Dezimalsystem; die Zahlen werden in der 
Form + a-10+° dargestellt, wobei a 11 Dezimalstellen und 6 2 Dezimalstellen hat. 
Das Speicherwerk enthalt eine rotierende Trommel mit einem magnetischen 
Belag, auf dem die Ziffern als Dipole aufgezeichnet werden; insgesamt kénnen 
10000 Zahlen und Befehle gespeichert werden. Die Rechenzeit betragt 20 ms fir 
die Multiplikation und 5 ms fiir die Addition. Dazu kommen noch gewisse Warte- 
zeiten fiir die Ablesung aus dem Speicher. Die Inversion einer zehnreihigen Matrix 
dauert 2,5 min. 

Die elektronischen Teile der Maschine enthalten Vakuumroéhren und Kristall- 
dioden; ferner werden Relais iiberall dort eingesetzt, wo sie die Rechengeschwin- 
digkeit nicht reduzieren. Der Verkehr mit der Aussenwelt geschieht durch Loch- 
karten und automatische Schreibmaschinen. Auf die Entwicklung neuartiger 
Schaltelemente wird im Interesse einer schnellen Fertigstellung und einer hohen 
Betriebssicherheit verzichtet; dagegen wird die logische Planung im Sinne einer 
Reduktion des Aufwandes an Schaltelementen sehr sorgfaltig durchgefiihrt. Die 


1) Die Herbsttagung der Schweizerischen Physikalischen Gesellschaft findet am 6. September 
1953 in Lugano im Rahmen der Generalversammlung der Schweizerischen Naturforschenden 
Gesellschaft statt. ‘ 

2) Institute for Applied Mathematics and Mechanics, University of New York. 

3) Assistenz fiir technische Mechanik an der ETH., Zurich. 

4) Institut fiir angewandte Mathematik, ETH., Ztirich. 
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vollstandige Maschine wird ungefahr 1200 Réhren, 2000 Kristalldioden und 300 
Relais enthalten. 

Gegenwirtig sind Versuche im Gang, welche die Einzelheiten der Schaltele- 
mente, insbesondere der magnetischen Speicherung, abklaren werden. Erstes Ziel 
dieser Versuche ist der Bau eines Speicherwerkes in teduziertem Mafstab, welches 
im Dauerbetrieb auf seine Zuverlassigkeit hin erprobt werden kann. 


Warmeiibergang im Atomreaktor, von W. TRAUPEL}). 

Das Problem der Warmeiibertragung im Atomreaktor ist deshalb schwierig, 
weil in kleinem Raum sehr viel Warme entwickelt wird, zu deren Ubertragung 
weder eine sehr grosse Oberflache noch eine grosse Temperaturdifferenz zur Ver- 
fiigung steht. 

Von den drei grundsatzlich méglichen Arten der Warmeiibertragung — Strah- 
lung, Leitung und Konvektion (Warmetransport durch bewegtes Medium) — 
kommt fiir die Abgabe der im Reaktor entwickelten Warme an ein warmever- 
wertendes System wesentlich nur die dritte in Frage. Stromt der Warmetrager 
laminar, so liegt allerdings insofern noch ein Warmeleitungsvorgang vor, als fiir 
jedes Raumelement die algebraische Summe der eingeleiteten und durch Stré- 
mung hereintransportierten Warme verschwinden muss. Strémt das Mittel durch 
einen Kanal, von dessen Wandungen her Warme eintritt, so erfolgt die Warme- 
zufuhr zum Raumelement rein durch Leitung, die Warmeabfuhr durch die str6- 
mende Bewegung. Es kommt in diesem Falle also lediglich auf die weitméglichste 
Steigerung des Warmeleitvorganges an, das heisst Wahl eines Mediums mit hoher, 
Warmeleitzahl und kleine Kanalabmessungen (kleine Leitungswege, grosse Tem- 
peraturgradienten). Kleine Kanalabmessungen sind ohnehin notwendig, damit 
iiberhaupt laminare Str6mung entsteht. Obwohl gerade die beliebige Steigerung 
des Warmeiiberganges durch Verkleinerung der Kanalabmessungen theoretisch 
sehr verlockend ist, stehen dieser Lésung bis heute grosse technologische Schwie- 
rigkeiten im Wege. Bekanntlich muss das Uran von einer dichten Hille umgeben 
sein, die den Ubertritt radioaktiver Spaltprodukte in den Warmetrager verhin- 
dert. Die Erfillung dieser Bedingung wird schwierig, wenn man sich etwa vor- 
stellt, dass feine Kiithlkandle in grosser Zahl die Uranmasse durchsetzen. Ausser- 
dem stellt das Zufiihren und Wegfiihren des Warmetragers bei der Vielzahl und 
Kiirze der Kiihlkanale grosse konstruktive Probleme. 

Praktisch kommt daher nur Warmeiibertragung durch turbulente Str6mung 
in Betracht. Hierbei tritt die Warmeleitung gréssenordnungsmassig in den Hin- 
tergrund gegentiber dem Warmetransport durch die turbulente Querbewegung 
im Kanal. Die Warmebilanz fiir ein Raumelement lautet dann so, dass die durch 
die Hauptstr6mung wegtransportierte Warme gleich ist der durch die turbulente 
Querbewegung einfliessenden. Diese ist aber proportional der Warmekapazitat 
oc, des Mediums und dem statistischen Mittelwert der Quergeschwindigkeiten. 
Letzterer ist in grober Naherung proportional der Geschwindigkeit w der Haupt- 
bewegung, so dass das Produkt @ c, w wesentlich massgebend fiir die Warmeiiber- 
tragung wird. Einer beliebigen Steigerung dieser Grésse durch Wahl einer hohen 


1) Gebriider Sulzer AG., Winterthur. 
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dissipation, auszukommen ist. Fliissigkeiten sind deshalb ungleich giinstiger als 
Gase. Unter der engen Auswahl, die uns die kernphysikalischen Bedingungen noch 
lassen, findet sich jedoch keine, die von schweren Nachteilen frei ware. Wasser 
muss unter hohem Druck gehalten werden, soll der Reaktor mit einer Temperatur 
arbeiten, die eine auch nur einigermassen wirtschaftliche Umwandlung der Warme 
in Arbeit gestattet. Fliissige Metalle (Pb, Bi, Sn) bieten Aussicht auf Erfolg, doch 
sind die technischen Schwierigkeiten gross. Gase sind nur bei hohen Driicken (zum 
Beispiel 40 at und mehr) aussichtsreich, weil dann @ einigermassen giinstig ist; sie 
haben den Vorteil, dass bei innen Druck und Temperatur voneinander unabhangig 
sind. Am giinstigsten sind leichte Gase, vor allem He, weil bei diesem C» hoch 
liegt, was das Verhaltnis zwischen oc, w und (@/2) w? giinstig beeinflusst. Grund- 
satzlich sehr giinstig ist der Warmeiibergang an eine verdampfende Fliissigkeit, 
weil bei ihr der Verdampfungsvorgang selbst die Turbulenz anfacht. 


Kinematische Ahnlichkeit bei Hydrozyklonen, von R. GreGoric!). 

Nach einem Hinweis tiber den zu grossen Aufwand zur experimentellen Ermitt- 
lung des Wirkungsgrades eines Hydrozyklons nach den strengen Gesetzen der 
Modellahnlichkeit werden einige vereinfachende Annahmen zur Reduktion der 
Kennzahlen gemacht. Die wesentlichste darunter ist die Unterdriickung der Erd- 
beschleunigung gegeniiber der Zentripetalbeschleunigung. Als neue Kennzahl 
wird, neben der Reynolds-Zahl der Fliissigkeitsstr6mung relativ zu den Wanden 
des Hydrozyklons, die Reynolds-Zahl der Relativstromung um das Korn einge- 
fiihrt. Das bis heute leider noch etwas sparliche Versuchsmaterial zeigt, dass die 
Annahmen nicht falsch waren. 


Calculs et mesures relatifs au col de tuyéres supersoniques, par B. CHAIx?) 
et P. Henrici), (communiqué par B. CHarx). 

Le calcul et la mise en service de nouvelles tuyéres supersoniques dans la 
soufflerie de l'Institut d’Aérodynamique ont confirmé que pour construire une 
tuyére assurant un écoulement uniforme dans la section de mesure, la connaissance 
exacte de 1’écoulement transsonique au travers de la section minimum est déter- 
minante. En effet, les perturbations dans l’écoulement supersonique ne s’amor- 
tissent pas, malgré la distance séparant les sections considérées. Les mesures et 
les calculs de l’écoulement dans la zone critique du col de la tuyére nous ont 
permis de construire des veines de mesure dans lesquelles la vitesse est uniforme 
a+ 1,5% prés. Cet écart peut étre considérablement réduit si l’on tient compte, 

.en dessinant le col de la tuyére, de la variation de la couche limite: Pour le calcul 
de cette derniére, des mesures complémentaires en régime turbulent ont été faites. 

Dans les écoulements plans, que nous considérons ici, le calcul graphique de la 
détente supersonique par la méthode des caractéristiques ne pose pas de problémes 
particuliers. Mais il convient, eu égard a la limite de validité de cette méthode, 

-d@’amorcer le véseau non pas sur la ligne sonique, mais en aval de celle-ci, dans le 
domaine nettement supersonique. Comme, sous l’influence de la convexité de la 
paroi, le fluide atteint la vitesse du son en amont du col, il régne déja une vitesse 
supersonique — a laquelle correspond un angle caractéristique ,> 0 — au point 
ot la paroi a une tangente horizontale. La caractéristique partant de ce pownt et 
aboutissant plus en aval sur l’axe de symétrie, pourra donc servir de condition 


1) Escher Wyss AG., Ziirich. 
2) Institut d’Aérodynamique, E.P.F., Zurich. 
3) Actuellement: American University, Washington. 
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initiale au calcul du réseau de caractéristiques. Elle décrit 1’écoulement au col 
d’une maniére particuliérement simple: La direction de 1’écoulement est la méme 
aux deux extrémités de cette ligne — horizontale, 4 la paroi comme sur l’axe -, 
et les pressions sont par conséquent aussi égales, puisque ces deux points se 
trouvent sur la méme caractéristique. Entre les points extrémes de cette ligne 
de Mach, que nous appelons S a cause de sa sinuosité, la pression varie, passant 
par un maximum au point d’inflexion du S. En ce point l’écoulement s’éloigne de 
axe de la tuyére, formant avec celui-ci un angle ’. Comme notre caractéristique 
S est traversée en son point d’inflexion par une caractéristique de l'autre famille 
partant de la limite sonique sur l’axe (# = 0; y’ = 0), on trouve a l’intersection: 
9’ = a,/2. L’angle w, rend donc l’essentiel du S. 

A part les mesures exécutées pour déterminer les conditions initiales S, des 
calculs ont été faits, qui décrivent le S dans une approximation semblable a celle 
des formules de SAUER pour la ligne sonique. On obtient entre autres: 


Pee See aay Vale 
on 5 Ve (el 


9 3 R 


ot y est la demi-hauteur du col, R le rayon de courbure de la paroi, et y exposant 
adiabatique. 

Pour préciser encore les relations entre la forme de la tuyére en amont du col 
et les conditions initiales de l’écoulement supersonique, un développement en 
série contenant deux paramétres permet de tenir compte non seulement de la 
courbure de la paroi au col, mais encore de la tangente au point d’inflexion amont. 
Une tuyére dessinée en partant d’un col ainsi calculé a fait ses preuves tout aussi 
bien que celles qui furent dessinées a partir d’un col préalablement essayé. Ces 
calculs transsoniques peuvent donc remplacer les mesures préliminaires pour la 
détermination des conditions initiales S de l’écoulement supersonique. 


Untersuchungen tiber Strahlumlenkung, von K. ISERLAND}). 

Das Gegenstiick zum Bremspropeller bei Propellerflugzeugen wiirde bei Di- 
senflugzeugen eine «Strahlbremse» darstellen, das heisst eine Vorrichtung, die 
erlaubt, den Triebwerkstrahl umzukehren und dessen Schub zur Bremsung zu 
verwenden, 

Im Hinblick auf die Entwicklung einer solchen Strahlbremse wurden im In- 
stitut fiir Aerodynamik der ETH. Versuche, vorerst an einem inkompressiblen 
Luftstrahl, durchgefiihrt. Auf Vorschlag von Professor ACKERET kam folgendes | 
System zur Anwendung: Dem Strahl wird vor der Diise Drall gegeben. Beim Aus- | 
tritt aus der Diise spreizt er dann auseinander und wird von halbtorusférmigen 
Ringen aufgefangen, die ihn umlenken. Ohne Drall strémt der Strahl in der Mitte 
der Umlenkringe durch, deren Innendurchmesser grésser ist als der Strahldurch- 
messer. Ein Kranz von drehbaren Schaufeln im Schubrohr erlaubt es, dem Strahl 
beliebig viel oder gar keinen Drall zu erteilen, das heisst auf Brems- oder Antrieb- 
schub einzustellen. 

Messungen mit einer solchen Apparatur ergaben eine gute Umlenkung des 
ganzen Strahls von einem gewissen minimalen Drall an. Es entsteht jedoch vor 
der Diise, im Schubrohr, ein Druckanstieg, welcher am Triebwerk eine Erhéhung 
des Druckes nach der Turbine, das heisst Gefahr des Pumpens des Kompressors 
bedeutet. Um diesen zu beseitigen, wurden Diisen mit verstellbarem Querschnitt 
verwendet, sowohl mit Nadel- wie mit Klappenregulierung. Dadurch war es 


1) Institut fiir Aerodynamik, ETH., Ziirich. 
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méglich, den Druck vor der Diise bis auf 1,7—-6 % (ja nach Anordnung der Ringe) 
konstant zu halten. Es wurden Bremsschiibe von bis zu 60 % des normalen Diisen- 
schubes erreicht. 

Uber einen Teil dieser Versuche ist in der «Interavia» (Heft Nr. 3, 1953) be- 
richtet. 


Mittel und Methodik der Literaturbearbeitung, von JAN R. DE FRIES}). 

Am Beispiel der Einflussfaktoren, die bei der konstruktiven Gestaltung von 
Gleitlagerungen zu beriicksichtigen sind, lasst sich zeigen, welch wertvolle Dienste 
ein zentraler, systematisch arbeitender Literaturnachweis dem Forscher und dem 
Konstrukteur erweisen kann, wenn die Voraussetzungen einer harmonischen, 
offenen Zusammenarbeit beider Instanzen vorhanden sind. Es sind in der Praxis 
zu einem recht grossen Teil psychologische Hemmnisse, die einer Anpassung der 
Literaturnachweistechnik an den Entwicklungsstand anderer Hilfsmittel im Wege 
stehen. Diese Hemmunisse sind zum Teil vielleicht das Ergebnis vieler unfrucht- 
barer Versuche, durch eine abstrakte Titeldokumentation, welcher der so not- 
wendige technische Realismus des Kommentars fehlt, zum Ziele zu kommen, aber 
auch die problematische Bereitschaft der Bentitzer zu einem echten «Teamwork » 
soll nicht unerwahnt bleiben. 

Die Aufgaben der Literaturbearbeitung stellen sich nach zwei Seiten: nach 
der Antwort auf eine spezifische, vom Beniitzer gestellte Frage, die fast immer 
erst der analytischen Gegenfrage, Ergriindung und Prazisierung des Problems, 
bedarf, um den Nachweis mit jener dokumentarischen Breite fiihren zu k6nnen, 
die zur Erganzung spezialisierten Wissens so notwendig ist — und die im Prinzip 
ahnliche analytische Stellung der Vielzahl von Antworten potentiellen Fragen 
gegentiber, die der standig gesichteten Literatur zu entnehmen sind. Durch das 
Bilden einer «Modellvorstellung» des beschriebenen Vorganges, durch Einteilung 
und Kommentar k6nnen die sehr fruc&tbaren Querverbindungen der technischen 
Erfahrung geschaffen werden, die die dankbarste Aufgabe systematischer Literatur- 
bearbeitung darstellen. Aber iiber dieses mit der Breite des Arbeitsgebietes stets 
etwas laienhafte Urteilen und Einteilen hinaus gestattet das System des standigen 
Sichtens auch das Erfassen von Veranderungen, die fast unmerklich in einem Gebiet 
der Naturwissenschaft oder der angewandten Technik eintreten und dem Leser 
vereinzelter Artikel kaum sichtbar werden, aber fiir die Disposition zukiinftiger 
Aufgaben entscheidende Bedeutung haben kénnen. 

Die Flut technischer Publikationen stieg in die Zahl von 3 Millionen pro 1950, 
eine Zahl, die mit aller Deutlichkeit beweist, dass die technische Dokumentation 
analytischer Richtung eine elementare Forderung unserer Zeit ist. 


Influence des pertes de charge sur la stabilité de la vitesse d’un groupe 
hydroélectrique, par L. Borev’). 

Le but de cette étude est la recherche des conditions de stabilité du réglage d’un 
groupe hydroélectrique sous l’influence simultanée du coup de bélier et de la 
perte de charge dans la conduite d’amenée. Le calcul se fait selon les méthodes 
usuelles, en appliquant le critére de Hurwitz a une équation différentielle linéaire 
du troisiéme ordre. 

Désignant par k = H,/H, la valeur relative de la perte de charge en régime 
permanent, par m le temps caractéristique du dosage accélérométrique, par t’ le 


1) Eidgendssische Materialprifungs- und Versuchsanstalt, Zurich. 
2) Ateliers des Charmilles, Genéve. 
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temps caractéristique de la promptitude de réglage, par i i le temps caractéristique 
de l’inertie mécanique du groupe, par 6 le temps caractéristique de Vinertie hy- 
draulique du systéme d’alimentation et posant: 

qaeil 


a 


’ 


Zz 
a=1-—2k, B= >a(1+ h), y=|/ 
il est possible d’énoncer les conclusions suivantes: 
1° Le systéme est instable, si 
k = ae 
c’est-a-dire si la perte de charge est plus élevée que la moitié de la chute nette. 
2° Le systéme est instable, si 


m ie i 
7] Ba 
3° Le systéme est stable, si 
(3 6/2)? 27 By 


et si, simultanément, le rapport m/6 est suffisamment voisin de la valeur optimum: 


eerie re 


Si nous admettons que cette derniére circonstance est toujours réalisée, la princi- 
pale condition de stabilité est: 


ML RE nd (> 0) . 


Le calcul de stabilité dans lequel la perte de charge est négligée conduit 4 une 
condition de forme identique a la forme ci-dessus, mais ot le coefficient K a la 
valeur 1,66. 

Die Herbsttagung der Schweizerischen Physikalischen Gesellschaft findet am 
6. September 1953 in Lugano im Rahmen der Generalversammlung der Schwei- 
zerischen Naturforschenden Gesellschaft statt. 
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Einfiihrung in die freie Geometrie ebener Kurven, Reihe «Elemente 
der Mathematik vom héheren Standpunkt aus», Band 1. Von L. LocHrr-Ernst 
(Birkhauser, Basel 1952). 87 S., 169 Abb.; sFr. 12.50. 


Das Biichlein gibt eine Einfiihrung in die Geometrie der reellen, projektiven 
Ebene, wie sie vor allem von C, JurEL entwickelt wurde, und vereinigt dabei in 
gliicklicher Weise mathematische Strenge mit einer geometrischen Anschaulich- 
keit, die durch eine Fiille schéner Figuren noch gesteigert wird.-Der Verfasser 
stiitzt sich auf ein geschickt gewahltes Axiomensystem fiir eine projektive Geo- 
metrie mit Anordnung und Stetigkeit, definiert zundchst in einfacher Weise den 
elementaren Kurvenbogen und die geschlossene Elementarkurve so, dass patho- 
logische Falle von vornherein ausgeschlossen sind, und untersucht dann neben 
den Kurvenbogen ohne Wendepunkte und Spitzen (Spiralenbogen) in erster 
Linie die Elementarkurven dritter Ordnung und diejenigen dritter Klasse. Der 
gesamte Aufbau ist dabei in sich dual gestaltet. 
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Wenn auch diese Dinge den Anwendungen etwas fernerliegen, so verdient das 
Buch trotzdem in weitesten Kreisen Interesse wegen der Fiille an Anschauungs- 
material zur projektiven Geometrie und insbesondere zum Dualitatsprinzip, die 
es enthalt. Es sei deshalb allen Lehrern und Studenten der Mathematik wirm- 
stens empfohlen! Druck und Ausstattung sind vorziiglich und ansprechend. 

G. Bol, Freiburg i. Br. 


Statistical Theory with Engineering Applications. By A. Harp (John 
Wiley & Sons, New York, 1952). 782 pp.; $9.00. 


A. HALD, der als Professor fiir Statistik an der Universitat Kopenhagen wirkt, 
hat hier ein sehr ausfiihrliches und umfassendes Lehrbuch der mathematischen 
Statistik verfasst, wobei er sich in erster Linie an den Ingenieur richtet. Die ma- 
thematischen Grundlagen werden in méglichst elementarer Form geboten; an 
mathematischen Hilfsmitteln wird nichts vorausgesetzt, was tiber die tibliche 
Differential- und Integralrechnung hinausginge. Alle theoretischen Erérterungen 
werden durch Beispiele von Anwendungen anschaulich gemacht. 

Als Handbuch der modernen statistischen Methoden kann das Werk von 
HALpD, vor allem auch dank seiner klaren Darstellung, lebhaft empfohlen werden. 

A. Linder 


Statistical Tables and Formulas. By A. Harp (John Wiley & Sons, 
New York, 1952), 97 pp., $2.50. 

Zu dem oben angezeigten Werk hat HaLp eine besondere Sammlung von 
Tafeln zusammengestellt, die fiir die Durchfiihrung statistischer Arbeiten uner- 
lasslich sind. Vorausgeschickt sind die mathematischen Formeln und Angaben 


uber die verschiedenen Anwendungsméglichkeiten der verschiedenen Tafeln. 
A, Linder 


Schwerkraft und Weltall. Von P. Jorpan (Verlag Vieweg & Sohn, 
Braunschweig 1952), 207 S., 4 Abb.; DM 15.80. 

Kosmologie und Kosmogonie als Lehre von Struktur und Entstehung des 
Weltalls beschaftigen seit den Anfangen der Wissenschaft die spekulative Phan- 
tasie. Grossen Auftrieb erhielten sie nach Aufstellung der allgemeinen Relativi- 
tatstheorie, da der hier zutage tretende Zusammenhang zwischen Raumstruktur 
und Materieerfiillung bestimmte Modelle fiir den Kosmos als Ganzes lieferte. 
Doch blieben diese Theorien wesentlich unbefriedigend, weil ihre Aussagen in 
deutlichem Widerspruch zu Erfahrungstatsachen (zum Beispiel beziiglich des 
Alters der Welt) stehen, falls man nicht die Einsteinschen Feldgleichungen in 
etwas willkiirlicher Art durch Addition des «kosmologischen Glieds» abandert. 

Ein weiteres Problem, das erst von EDDINGTON als solches erkannt wurde 
und auf welches die allgemeine Relativitatstheorie keine Antwort gibt, ist die 
ausserordentliche Kleinheit der Gravitationskonstanten, das heisst die Tatsache, 
dass die Gravitationsanziehung der Teilchen etwa im Wasserstoffatom 2 - 103°- 
mal kleiner ist als die elektrische Anziehung. EDDINGTON selber versuchte, diese 
Zahl 2 - 1039 mit der auf 107° geschdtzten Anzahl Nukleonen im Weltall in Ver- 
bindung zu bringen und eine Deutung fiir diese riesigen Zahlen zu geben, die aber 
als misslungen zu betrachten ist. Ein neuer Gedanke wurde durch Dirac an das 
Problem herangetragen: Drrac vermutete (1937), dass jene beiden Zahlen nicht 
Naturkonstanten, sondern Zeitfunktionen seien, Funktionen des in Elementar- 
zeiten gemessen etwa 108° betragenden Weltalters. Dies bedeutet nichts anders, 
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als dass die Gravitationskonstante f selbst nicht als Konstante, sondern als 
Feldgrésse anzusehen ware. P. JoRDAN griff diesen Gedanken auf und entwik- 


kelte eine Verallgemeinerung der allgemeinen Relativitatstheorie, in der neben ~ 


den metrischen und elektrischen auch f als Feldgrésse erscheint. Einer der scho- 
nen Zufalle in der Geschichte der Wissenschaften — sind es wirklich Zufalle ? — 
wollte es, dass er das notwendige Werkzeug dazu fertig geschmiedet vorfand in 
der aus einem ganz andern Problemkreis entstandenen projektiven Relativitats- 
theorie. 

Bekanntlich wurde bald nach E1nste1ns Entdeckung des Zusammenhangs 
zwischen Metrik und Gravitation verschiedentlich versucht, auch das elektro- 
magnetische Feld zu «geometrisieren ». Einer dieser Versuche war der von KaLuza 
(1921), die Welt durch ein fiinfdimensionales Koordinatenkontinuum darzustellen. 
Es zeigte sich dann, dass die natiirliche Interpretation dieser fiinfdimensionalen 
Koordinaten war, sie als projektive Koordinaten fiir ein vierdimensionales Kon- 


tinuum aufzufassen, wobei aber eine quadratische Metrik in den projektiven ~ 


Koordinaten selber definiert ist. Diese Metrik spaltet beim Ubergang zu vier- 
dimensionalen (affinen) Koordinaten in eine vierdimensionale Metrik g;,, ein 


schiefes Tensorfeld F,,, und eine Invariante J; die den Einsteinschen analogen © 


fiinfdimensionalen Feldgleichungen spalten bei Hinzunahme einer Nebenbedin- 
gung J = 1 genau in das Einstein-Maxwellsche System von Gleichungen fiir die 
Gravitationspotentiale g;, und das elektromagnetische Feld F;,. 


Jorpans Entdeckung bestand nun darin, dass man bei Fallenlassen der Be- 


dingung J = 1 eine erweiterte Theorie erhalt, die gerade die Eigenschaft hat, die 
Gravitationskonstante — im Wesentlichen J — als neue Feldgrésse einzufiihren. 
Allerdings bleibt eine gewisse Willkiir in den Feldgleichungen bestehen, die sich 
jedoch weitgehend reduzieren lasst. 

Mit dieser erweiterten Theorie erhalt nun auch die kosmologische Frage eine 
neue Beleuchtung. Neue Modelle fiir den Kosmos werden méglich: insbesondere 


eines, das die Diracsche lineare Abnahme der Gravitationskonstanten mit der — 


Zeit liefert, sowie ein anderes, das vielleicht die Spontanentstehung von Sternen 


beschreibt, ein Prozess, der in solchen Theorien nétig ist, da der Energiesatz — 


kosmologische Modifikationen erleidet. (Gerade hier hat die Theorie allerdings 


eine wesentliche, von JORDAN selber betonte Liicke, indem die Neuentstehung ~ 
von Materie nach den Grundgleichungen eigentlich durch stetiges Wachstum der ~ 


vorhandenen Massen erfolgen sollte, was mit der Erfahrung unvereinbar zu sein 
scheint. JORDAN vermutet die Existenz einer in den Gleichungen nicht ausge- 
driickten Relaxation, welche diesen kontinuierlichen Prozess verhindert und 
spontane Entstehung von Sternen an seine Stelle setzt.) Dass diese Modelle im 
einzelnen noch ebenso unzulanglich sind wie die der bisherigen Theorie, ist nicht 
sehr schwerwiegend; nach der neuen Theorie existieren deren noch viele andere, 
die zu untersuchen eine Aufgabe der Zukunft ist. 

Diese ganze, bisher in Originalarbeiten verstreute, nur miihsam iiberblickbare 
Entwicklung wird im vorliegenden Werk ihres Hauptférderers in souverdner 
Weise zusammenhangend dargestellt. JoRDAN bringt das Meisterstiick zustande, 
in mathematisch einwandfreier und vollstandiger und dennoch leicht lesbarer 
Form die (gewéhnliche) allgemeine Relativitatstheorie (inklusive ihrer mathema- 
tischen Grundlegung, der Riemannschen Geometrie), die projektive Verallgemei- 
nerung samt Anwendungen und ausfiihrlicher Diskussion in einem Band von nur 
200 Seiten unterzubringen. Sicher sind die hier dargelegten Theorien noch unfer- 
tig; doch sind sie in voller Entwicklung begriffen und haben bereits einige theo- 
retische und empirische Erfolge zu buchen. Jedenfalls ist das Buch eine anre- 


~~ 
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gende, ja spannende und lohnende Lektiire fiir jeden Physiker und Astronomen; 
fur den, der sich naher mit der kosmologischen Forschung befasst, bildet es eine 
dauernde Quelle und unentbehrliche Hilfe. M. R. Schafroth 


Einfihrung in die theoretische Gasdynamik. Von R. SAUER (Springer- 
Verlag, Berlin 1951). 174 S., 107 Abb.; DM 16.50. 


Die 1943 erschienene ausgezeichnete Einfiihrung in die theoretische Gasdyna- 
mik von Prof. SAUER ist in demselben Verlag in neuer Auflage erschienen. Diese 
enthalt zahlreiche kleinere Erweiterungen und ist ganz allgemein auf den heu- 
tigen Stand der Technik gebracht worden. Wesentlich ist dabei die Einbeziehung 
des Charakteristikenverfahrens fiir nicht wirbelfreie Uberschallstr6mungen und 
ein ganz neuer Abschnitt tiber dreidimensionale Probleme, insbesondere nicht 
rotationssymmetrische Uberschallstrémungen um Drehk6rper mit Anstellwinkel 
und endliche Tragfliigel. (Siehe auch Rezension der franzésischen Ubersetzung: 
Ecoulements des fluides compressibles, ZAMP 3, Fasc. 4, 318 [1952].) 

Die neue Auflage von SAuERS Buch, die sich wiederum an die Ingenieure 
wendet, welche sich fiir die mathematische Behandlung der stationaren Stré- 
mungen bei hohen Geschwindigkeiten interessieren, ist durch seine klare und 
ubersichtliche Darstellung gekennzeichnet. Es ist ein vorziigliches Buch, dessen 
Studium jedermann, der sich in das Gebiet der modernen Gasdynamik ein- 
arbeiten will, zur Freude gereichen wird. R. Sdnger 


An Introduction to Acoustics.. By R. H. Randall (Addison Wesley Press, 
Inc., Cambridge, Mass., 1951). 340 pp., 140 figs; $6.00. 

Wenn der Verfasser im Vorwort feststellt, dass sich sein Buch in theoretischer 
Beziehung an Morses Vibration and Sound anlehnt und in praktischer Beziehung 
an BERANEKS Acoustic Measurements grenzt, so hat er es damit sehr zutreffend 
charakterisiert. Es ist ein Lehrbuch im besten Sinne des Wortes mit vielen mo- 
dernen Beispielen, die seine Lektiire sehr lebendig und anschaulich gestalten. 
Gleichzeitig wird auf diese Art ein Querschnitt durch den heutigen Stand der 
Technik vermittelt, wobei man natiirlich keine Vollstandigkeit erwarten darf. 

Das Buch ist in erster Linie fiir Studierende der Ingenieurwissenschaften und 
der Physik bestimmt. Es enthalt deshalb, wie das bei amerikanischen Biichern 
haufig iiblich ist, am Ende jedes Kapitels eine Anzahl typischer und gut ausge- 
wahlter Aufgaben. Die theoretischen Grundlagen sind einfach, klar und sauber 
dargestellt, ebenso die praktischen Beispiele, und tiberall ist groésster Wert auf 
physikalische Anschaulichkeit gelegt. Durchweg wird das cgs-MaBsystem ver- 
wendet, was erneut zeigt, dass die Einfiihrung des Giorgi-Systems in die Akustik 
einigen Hemmungen begegnet, die offenbar darauf zuriickzufiihren sind, dass 
hier das cgs-System direkt als technisches MaBsystem dient. Die trotzdem gele- 
gentlich anzutreffenden Angaben in Fuss und Zoll sind kleine Schonheitsfehler, 
die man anglo-sdchsischen Autoren tolerieren muss. ; 

Das Buch ist sehr sorgfaltig gedruckt und hervorragend ausgestattet und ist 
sicherlich eines der besten Lehrbiicher, die auf diesem Gebiete in den letzten 
Jahren erschienen sind. Aus dem reichhaltigen Inhalt seien die folgenden Punkte 
kurz erwahnt: Im ersten Kapitel wird ein sehr sch6ner Abriss tiber die allgemeine 
Schwingungslehre gegeben, wobei die Kenntnis der elementaren Differential- und 
Integralrechnung vorausgesetzt wird. Zwei weitere Abschnitte sind den ebenen 
und raumlichen Wellen sowie der Strahlung gewidmet, an die sich ein sehr schoner 
Teil iiber die so wichtigen Interferenzprobleme anschliesst; hier werden einfache 
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Beugungsprobleme, die in der Akustik so bedeutsam sind, behandelt. Weiter 
folgen Abschnitte tiber die akustische Impedanz, von welcher zu den H6érnern 


und Resonatoren und schliesslich zu der Wellenausbreitung iibergeleitet wird. — 


Besondere Erwahnung verdient dabei der Abschnitt iiber akustische Linsen. 
Diese schon an und fiir sich reichhaltige erste Halfte wird durch weitere 
Kapitel abgerundet, von deren Inhalt die folgenden Stichworte einen Begriff ver- 


mitteln: Saiten, Membranen, Musikinstrumente; Reflexion, Absorption; Sprache ~ 


und Gehér; Messmethoden; Schallwidergabe; Raumakustik; Ultraschall, Was- 
serschall, Materialpriifung. 
Dieser zweite Teil des Buches ist notwendigerweise sehr summarisch gehalten 


a 


und dient lediglich dazu, das Bild des Gebietes zu erganzen und abzurunden und © 


gegebenenfalls zum Studium der einschlagigen Spezialliteratur anzuregen. 

Das Buch kann nicht nur dem Studierenden, sondern auch dem praktisch 
tatigen Ingenieur und Physiker, der sich mit Schallfragen in irgendeiner Form 
zu befassen hat, warmstens empfohlen werden. W. Furrer 


Antennas, Theory and Practice. By S. A. SCHELKUNOFF and H. T. Frus 
(John Wiley & Sons Inc., New York 1952). 639 pp., $10.00. 

Dieses Werk iiber Antennen, verfasst von zwei berufenen Fachleuten der 
Bell Telephone Laboratories, zerfallt in drei Teile. Der erste Teil (Kapitel 1) gibt 
uns in Form eines Buches im Buche einen umfassenden Uberblick tiber das 


oe 


gesamte Gebiet der Antennen, der zweite Teil (Kapitel 2 bis 9) schafft ohne Ver- | 


wendung der Vektoranalysis und unter Beschrankung zeitraubender Ableitungen 
bekannter Grundgleichungen das Fundament fiir die im dritten Teile behandelten 
Antennen und Antennensysteme. Dieser zweite Teil umfasst, neben Definitionen, 


ausgehend von den Maxwellschen Gleichungen, die Probleme der ebenen und der 


spharischen Wellen, die gerichtete Strahlung und die Reflexionen an der Erde, 
die Stromverteilung in Antennen, die Impedanzen und Aquivalenzen von Anten- 
nen und die Berechnung wirksamer Antennenflachen und des Antennengewinnes. 
Der dritte Teil (Kapitel 10 bis 19) wendet nun die in den ersten beiden Teilen 
des Buches abgeleitete Theorie auf spezifische Antennen und Antennensysteme 
an. Kapitel 10 behandelt die sogenannten «kleinen Antennen» (Lineardimensionen 
der Antenne kleiner als ein Achtel der Wellenlange), Kapitel 11 die eigenresonan- 
ten Antennen (Dipolantennen), Kapitel 12 ist einer allgemeinen Theorie linearer 
Antennen gewidmet. Kapitel 13 befasst sich mit der Impedanz von Dipolanten- 
nen, wahrend das Kapitel 14 bzw. 15 die fiir Kurzwellenverbindungen wichtigen 
Rhombusantennen bzw. die linearen Antennensysteme erlautert. Die letzten 
vier Kapitel, 16 bis 19, behandeln zur Hauptsache Mikrowellenantennen, Horn- 
antennen, Schlitzantennen, Reflektoren und Mikrowellenlinsen. 

Neben einer ausgewahlten, Publikationen bis 1951 umfassenden Bibliographie 
enthalten eine Anzahl von Kapiteln das Selbststudium anregende Probleme nebst 
ihren Loésungen. Eine Zusammenstellung der Biicher, welche Antennen und ihre 
Probleme erértern, ferner eine fiir den Praktiker sehr wertvolle Zusammenfas- 
sung von Angaben und Formeln iiber Leitungen, Dipolantennen, Antennen- 
systeme, Hornantennen und Linsen sowie ein nach Namen und Sachgebieten 
geordneter Index runden das Buch ab. 

Das Werk von SCHELKUNOFF und Frits ist jedem Leser, welcher sich mehr 
als nur oberflachlich mit Antennenproblemen vertraut machen will, als Lehrbuch 
wie auch als Nachschlagewerk sehr zu empfehlen. Ich zweifle nicht daran, dass 
dieses Buch innert kiirzester Zeit in jeder Fachbibliothek den ihm gebiihrenden 
Ehrenplatz einnehmen wird. H. Hagger 
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